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Abstract. We study the asymptotic behavior of oscillatory Riemann-Hilbert 
problems arising in the AKNS hierarchy of intcgrable nonlinear PDE's. Our 
method is based on the Dcift-Zhou nonlinear steepest descent method in which 
the given Riemann-Hilbert problem localizes to small neighborhoods of station- 
ary phase points. In their original work, Deift and Zhou only considered ana- 
lytic phase functions. Subsequently Varzugin extended the Deift-Zhou method 
to a certain restricted class of non-analytic phase functions. In this paper, we 
extend Varzugin's method to a substantially more general class of non-analytic 
phase functions. In our work real variable methods play a key role. 



Contents 

1. Introduction 

1.1. Notational conventions 

2. Outline of the proof of Theorem ll.il 

3. Beals-Coifman operator formulation for RHPs 

4. Localization schemes and essence of the nonlinear stationary phase 

method 

4.1. The direct scheme 

4.2. The indirect scheme 

4.3. Combination of two schemes 

4.4. Essence of the nonlinear stationary phase method 

5. Hardy decomposition of oscillatory functions 

6. Factorizations of the oscillatory jump matrix J(A, i) 

6.1. Existence and boundedness 

6.2. Bounds on derivatives of 5 

6.3. Approximation of S 

7. Unique solvability of oscillatory RHPs 

7.1. Unique solvability of the original RHP 

7.2. Unique solvability of intermediate RHPs 

8. Reduction to model cases (I): Localization and phase reduction 

8.1. Localization to neighborhood of stationary points 

8.2. Reduction of phase 

8.3. Separation of contributions 

9. Reduction to model cases (II): Deift-Zhou's steepest descent argument 

9.1. Preparation for steepest descent 

9.2. Steepest descent reduction to the T-deformed model case 

10. Model Riemann-Hilbert problems 
10.1. Degenerate and defocusing cases 



2 



YEN DO 



10.2. Focusing cases 
References 



1. Introduction 

In many studies of asymptotical behaviors of nonlinear systems, the following 
Riemann-Hilbert problem (RHP) appears: given J an invertible matrix- valued func- 
tion on R, we look for a matrix-valued function M (A) which is analytic on C\t 
and satisfies the following jump condition: 

Af+(A) = M_(A)J(A), for a.e. A e 1 

M±(A) denote the upper and lower parts of M, with appropriate normalization 
at oo and non-tangential limits on R (denoted by the same symbols). The jump 
matrix J could be thought of as the multiplicative gain of M when moving from 
the lower half plane to the upper part. 

In this paper, we will be only interested in the 2x2 case since this appears in 
most of the applications. In the natural L 2 setting, the following normalization 
conditions are often imposed: 

(1) M+(X) -Ie H 2 (C+), M_(A) - I e ff a (C_), 

where i?2(C+) is the complex Hardy space H 2 of the upper half plane and /?a(C_) is 
the complex Hardy space H% of the lower half plane, and / is the identity matrix. It 
is not hard to see that these conditions are equivalent to the existence of h E L 2 (R) 
such that M± = I + C±(h). An overview of L p Riemann-Hilbert theory is available 
in [10]. 

For many equations in the AKNS hierarchy of integrable nonlinear PDEs, the 
jump matrix J has the following oscillatory structure pQ: 

l+p(\)q(\) p{X)e- lte ^ 
? ( A ) e it0(A) i 

Here p(A) and q(X) are the reflection coefficients of the initial data for the corre- 
sponding PDE (often q = ±p), and p, q will be sufficiently smooth with sufficient 
decay if the initial data is sufficiently nice (see [33] for further details). The phase 
is real-valued and depends on the algebraic structure of the current PDE. The 
precise technical assumptions on p, q, used in this paper will be specified later. 
The oscillatory factors e ±lte originate from the evolution of the reflection coeffi- 
cients when we let the initial data evolves according to the PDE. Several examples 
include 0(A) = (A — A ) 2 for NLS and 0(A) = 4(A 3 - 3AqA) for mKdV [J. 

In the setting of ([1]) and @, the question of interests is the long-time behavior 
of the potentials 

(3) u{t) = lim AA/i 2 (A,t) 

A— >OG 

(and similarly v(t) = lim^-i-oo AM2i(A, t)). Here in the limits A — > oo non-tangentially 
in C. In the AKNS setting, the function u(t) recovered from ^ is a constant mul- 
tiple of the solution u(x, t) to the respective nonlinear PDE, therefore the above 
question corresponds to the long-time behavior u(x, t). We note that the spatial 
parameter x in u(x,t) has been encoded in the phase and hence it is encoded in 
the stationary points {A : 0'(A) = 0}, and in this paper the stationary points will 



(2) J(A,<) := 
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be assumed constant (see also the technical assumptions on 9 and p, q below). This 
is the reason why we will simply write u(t) and v(t) instead of u{x,t) and v(x,t). 
For simplicity of notation, we'll often suppress t when writing M in this paper. 

Oscillatory RHPs also appear in other settings, in which the structure of the os- 
cillatory jump matrix J may be different from ([2]) and the Riemann-Hilbert solution 
may no longer have the L 2 normalization ((T|). For instance, the small-dispersion 
limit e —> + of the KdV equation 

u t - 6uu x + e 2 u xxx = 

corresponds to the large 1/e asymptotics of an oscillatory RHP with a different 
^-normalization [TSl [2B] - Other examples include the (non L 2 )Riemann-Hilbert 
formulation of Fokas, Its, Kitaev for orthogonal polynomials with varying weights 
e~ nV ( x 'dx [14] whose asymptotics plays an important role in random matrix theory 
(see 127] and the references therein, for related results involving random matrix 
theory and orthogonal polynomials see [241 E3] and the references therein) . 

The study of long-time behaviors of solutions to nonlinear integrable PDEs goes 
back at least to the work of Zakharov and Manakov [32], who were probably the 
first to write down a correct formula for the leading asymptotics of u(t) (and v(t)) 
in the NLS case. Using monodromy theory, Its [T^] was able to reduce the RHP 
formulation for NLS to a model case, which can then be solved explicitly, giving 
the desired asymptotics for u (and v). 

Dcift and Zhou [9] developed a rigorous nonlinear steepest descent method to 
study the oscillatory RHP associated with the mKdV equation. Their argument de- 
mands analyticity, at least for 9(X) (the reflection coefficients can be approximated 
by suitable analytic functions [10]). For the mKdV equation there is more than 
one stationary point and Deift and Zhou were able to separate their contributions 
using an operator formulation for RHPs that goes back to Beals and Coifman pQ . 

The first work in the non-analytic setting is due to Varzugin [31] in the set- 
ting (|1I2|) for 9 with stationary points of first order. Varzugin's asymptotics was 
previously ansatzed by Kitaev [18] using a method of isomonodromic deformation. 

The steepest descent method of Deift and Zhou has been extended in various 
ways. In their study of the small-dispersion limit of the KdV equation, Dcift, 
Venakides, and Zhou [7] discovered that the contribution to the asymptotics of the 
solution u(x, t, e) recovered from the solution of the corresponding RHP (via a limit 
similar to (|3]>) comes from a family of intervals instead of the isolated stationary 
points of the phase in the original RHP. The key idea is to introduce a ^-function 

such that after conjugating the given RHP with e l9 ^ a3 / t where 173 = is 

the third Pauli matrix, one arrives at an RHP whose jump matrix 

T — p -' icr 39-( A )/ e 7 , ,pi^39+W/e 
"new — ^ tJ ola^ 

is ready for a steepest descent argument (in Section l4~4l we shall discuss this further). 
Furthermore, a systematic procedure for determining the intervals contributing to 
the asymptotics of the Riemann-Hilbert problem was also developed. The methods 
of Deift, Venakides, Zhou have been successfully applied to the orthogonal polyno- 
mial setting [6] , which led to deep results in random matrix theory. The argument 
however still requires (real) analyticity of the phases of the new oscillatory jump 
matrices. In the dispersion KdV setting, this goes back to analyticity of the initial 
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data and in the orthogonal polynomial setting it goes back to analyticity of the 
underlying weights. 

The second extension of the steepest descent method of Deift and Zhou is the 
d-steepest descent method of McLaughlin and Miller, which first appeared in the 
orthogonal polynomial setting [26l [27] . The d method follows the general scheme 
of the steepest descent argument of Deift- Venakides-Zhou, however non-analytic 
data are now continued to the desired contours via the solution of a d equation. In 
particular, Green's theorem (which works for smooth functions) is used as a remedy 
for Cauchy's theorem (which requires analyticity), and the error terms are under 
control if appropriate constraints are imposed on the d problem. The 9-steepest 
descent method of McLaughlin and Miller allows for non-analytic phases with two 
Lipschitz derivatives near stationary points. In this method was adapted to 
the AKNS setting and was used to obtain long-time asymptotics for solutions to 
the defocusing NLS, with sharp error bound for initial data in H 1:1 (R) — {/ £ L 2 : 

/'. '•/ ( /.-!, 

The main goal of this paper is to further demonstrate that the asymptotics 
of oscillatory RHPs can be studied using analogue of real- variable tools from the 
linear theory of oscillatory integrals. We shall focus on the L 2 -normalized RHPs 
(TTJ) with oscillatory jump matrix of the form © (which corresponds to the study of 
long-time behavior of solutions to many nonlinear integrable PDEs) , however the 
author anticipates applications of the argument developed in this paper to other 
settings (especially those mentioned above where nonanalytic data are frequently 
encountered). Due to the length and the scope of the current paper, the extent 
of the applications of the argument is not fully explored, but we shall discuss the 
underlying ideas that can be used for these potential adaptations in Section 14.41 
The analysis used in this paper extends the work of Varzugin |31j . and for the 
RHP |l] [2]) our argument allows for non-analytic phases with a finite number of 
stationary points (of arbitrary orders), under fairly minimal regularity assumptions. 

To see a simple connection between the setting of (fTJ [5]) and the linear theory of 
oscillatory integrals, consider a degenerate situation when q(X) — 0. In this case, 
the oscillatory RHP (p] [2} becomes abelian and can be solved using the Hilbert 
transform. Recovering u(t) via the limit ([3]), we arrive at the linear problem: study 
the asymptotics behavior of J p(X)e~ lte ^dX, and this could be achieved by the 
classical stationary phase method [30] . which exploits cancellation resulted from 
rapid oscillation of e~ lte ^ away from the stationary points of 9. For this reason, 
the presentation of the proof of the main result of this paper will follow the spirit 
of the linear stationary phase method and is summarized in Section [2j 

In this paper, the following assumptions will be made on 9: 

(A) 9 is real valued and continuously differentiable, and on the complement of a 
finite set of points 9 has three locally integrable derivatives. 

(B) 9 has stationary points Ai, . . . , Xn of orders fci, . . . ,fcjv he- 9 is (kj + l)-time 
differentiable near Xj and 0'(Aj) = • • • = 6^(Xj) = ^ 9^ +1 \Xj). 

(C) If A is a stationary point of order k of 9 then we require 9^ k+1 \x) to be 
Holder continuous at A with some exponent (3 > 0, i.e. near A 



9^ k+1 \x) = 9^ k+1) (X) + 0(\x - Xf). 
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If k = 1 then we assume 9^ is integrable near A in some high L r , more precisely 
r > 4 (by Holder's inequality j3 > 1/r' , so this is automatic if r > 2). If fe = 2 and 

p(A)g(A) < then we assume 9^ is Lipschitz near A. 

For convenient, we say that a stationary point A is focusing if p(X)q(X) > and 
defocusing if p(X)q(X) < 0. Wc note that (C) is weaker than having two Lips- 
chitz derivatives near primary stationary points, which by Rademacher's theorem 
is equivalent to having a bounded 9^ there. The Lipschitz assumption of 9^ near 
defocusing secondary stationary points can also be weakened to a similar high L r 
integrability condition of 9^\ and is chosen here for simplicity of the argument. 

For k,j € Z + we recall the Sobolev space H k ' j := {/ e I? : /',..., f^ k) S 

/.-••'•/ <••''./••. /•-'!• 

The general assumption for p, q will be: 

(D) p,q £ H 1>0 (R) such that < 1 + p(x)q(x) = 0(1), and p, q have sufficiently 
decay (which depends on 9). 

(E) p, q have two I? derivatives near every stationary point A, and if its order 
k > 3 then p(X)q(X) < 1. If k > 3 or if k = 2 and p(X)q(X) < then we require a 
third L 2 derivative for p, q near A, and the assumption in (D) will be changed to 
p, q G i/ 2,0 overall. 

The regularity assumptions on p, q can be improved. For example, if N = 1 and 
Q(ki+i) j s cons tant near Ai (such as in the case of the NLS) then we only need one 
L 2 derivative for each p and q if Ai is a primary stationary point or a focusing 
secondary stationary point. For a stationary point A of order k, the requirement 
p(X)q(X) < 1 may be removed if a certain model RHP associated with stationary 
points of order k is better understood; the author plans to revisit this issue in a 
future manuscript. 

In this paper we will not try to optimize the decay requirements for p and q, this 
however could be done for explicit 9. 
For convenience of notation, denote 



1 1 0, kj even; 

-ln[l + KA i)9 (A,)], <iH^ (fl (M. l)(A))| fc . odd . 



Theorem 1.1. 3to such that the RHP ([7J H[) has unique solution for t > to- 
Furthermore, as t — > oo: 

(i) If N = then the recovered potentials u{t) and v(t) satisfies 



u(t),v{t) = O e {t- 1+e ), e>0 



This can be arbitrarily improved if stronger regularity assumption on p, q, 9 are 
given; 

(ii) If N > 1 then there are dj > such that 

TV / v 

M(<A -E[Q) u f) + o(t-^^) 



The contributions of Xj are of the form: 

(4) Uj(t) = UjpiX^r 1 ^ 1 exp(-i[td(Xj) + ajlnt-2u!j]j 

(5) Vj{t) = V j q{X j )t~ % i+ T exp (i [t9{X ) + ay ln< - 2uj]j 
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here Uj,Vj depend on Xj, kj, Qw +1 ) [Xj) and p(Xj)q(Xj), and Uj,Vj are anti- 
complex conjugates, i.e. 



Uj = -Vi 



while aj , ujj are real numbers and can be explicitly computed by 



2e j u j 



1 



— . In \Xj - y\dhx[l +pq](y) + e k v k In \Xj - \ k \ 

(the right-hand side for ujj makes sense if ln(l +pq) has sufficient decay, say H 1,1 , 
but ujj still exists under weaker assumption, see Section\$ and in particular h32\) ). 

Notes: 1. Theorcm ll.ll is related to a question in [5]. The unique solvability for large 
t of our oscillatory RHP can be proved for a larger class of p, q (see Theorem 1 7. 1[) . 
A similar form of the above asymptotics when 9 has one primary and one secondary 
phase point was previously given as an ansatz by Kitaev |18j . 

2. Compared to linear theory, there is an extra lni term in the exponent of the 
leading asymptotics of u(t) and v(t). There are also interactions between stationary 
points, which can be observed in the definition of uij. 

3. The error estimate 0(t ^ fe J + 1+dj - ) ) can be made explicit. For instance, if 
r = oo and fti = ■ ■ ■ = = 1 then we can take dj = gnr+I) — e ' furthermore if 
N = 1 then this can be improved to fc x +1 — e. In the NLS case, these estimates are 
well-known from the work of Deift and Zhou [10] . 

4. If kj = 1 then Uj (hence Vj) can be computed explicitly (see for instance 

my- 

ivj 1 /«7re 7 - ,\ 

Uj = — = — = exp — -i + i arg r(ze,i/, j 

(when p(Xj) = or q(Xj) = this should be understood in the limiting sense). 

1.1. Notational conventions. In this paper, we make the following conventions: 

1. Absolute value of a matrix: For any matrix M define \M\ — (trM* M) 1 / 2 
and for any matrix valued function A on R define ||A|| P := |||A||| . This makes 
the respective LP space a Banach space for 1 < p < oo. We'll abuse notation and 
refer to this space also as LP - it should be clear from the context what is being 
referred to. Notice that the above absolute value |.| satisfies the triangle inequality, 
furthermore \AB\ < \A\\B\ for any A, B. 

2. Action of operators: The action of any operator T on matrices is done entry- 
wise. If T is a bounded operator on LP then it is also bounded as an operator on 
the space of 2 x 2 matrices having LP entries, with comparable norm. 

3. Differentiability: / is said to be /c-time diffcrcntiable if / is (fc — l)-time 
continuously differentiable and f^' is locally integrable. 

4. Inequalities up to a constant: For two quantities A and B we say A < B 
if there exists an absolute constant C > such that \A\ < CB. If the constant 
C depends on the parameters p±, . . . ,p n we shall say that A < Pl ,..., Pn B. In some 
situations, when the dependence of C on certain parameters are not important for 
the proof or the discussion, we shall abuse notation and suppress those parameters 
in writing, in particular if the dependence of C on all the parameters are not im- 
portant we will simply write A < B. The implicit constants used in the manuscript 
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are not necessarily the same and may differ from line to line. The meaning and 
conventions used for > and ~£p lt ...,p n are understood similarly. 

5. oo— , p+, p—: For the sake of brevity of the argument, we shall often write oo— 
to denote a hnite number that is very large. For any finite number p we shall use 
p+ to denote a finite number q > p such that q — p is sufficiently small. Similarly, 
p— shall denote a finite number q < p such that p — q is sufficiently small. This 
notation will be useful when applying Holder inequality, for instance we can write 
Il/Slb < ll/lb+Hslloo- 

6. Sobolev spaces: For k, j E Z+ we shall denote H k >' := {/ 6 L 2 :/',... , f {k) £ 
L 2 ,x/,...,xV€L 2 }. 

2. Outline of the proof of Theorem 11.11 

In this section, we describe the main ideas used to prove Theorem ll.il 

One fundamental tool that will be used throughout is the operator formulation 

by Beals and Coifman for solutions to RHPs. Given a pair of L°° weights (w~, u> + ), 

Beals and Coifman introduced the following operator pQ 

C w f = C+(fw-) + C-(fw + ) 

here C + = H + , C_ = — H- and H + ,H^ are respectively the projections into the 
complex Hardy space H P (C+), H P (C-) of the upper and lower half planes. This 
operator can be thought of as a weighted version of the identity operator in the 
Hardy decomposition of L 2 

/ = C+(/)-C_(/) 

Usually, the pair of weights (w~ , w + ) is obtained from a factorization of the jump 
matrix 

J = {I -w-)-\l + w + ) 

Except for model cases that will be discussed below, the weights in this paper are in 
L 2 (in addition to being in L°°). In that case, if 1 — C w is invertible on L 2 then the 
respective RHP will be uniquely solvable and its solution can be formulated using 
C w pQ (see Section [3]) ■ Beals-Coifman's operator formulation is useful in studying 
perturbation and unique solvability of RHPs. 

Our proof of Theorem ll.il will be a sequence of reductions following the spirit of 
the classical stationary phase method. In each reduction, we reduce one RHP to 
another RHP such that, heuristically, if Theorem II .11 is true for the new RHP then 
it is true for the current RHP. More precisely, we will show that, for t large, unique 
solvability of the new RHP implies that of the current RHP, and the respective 
potentials u(t),v(t) differ by terms of highly decaying order of t. Note that for 
some of the intermediate RHPs that arise during our reductions, the defining limits 
for u, v may not exist (most likely the limits in the upper and lower half-planes 
don't agree). Denoting by AM the difference between the solutions of the old 
and new RHPs, we overcome this issue by showing that the nontangential limits 



lim sup 2 



decay strongly as t —> oo. 



zAMi 2 (z) and limsup^^ zAM 2 i(z) 

Since lim sup is sub- additive, these estimates are sufficient. 

We will start from the given oscillatory RHP whose jump matrix is defined in ([2]) 
and finally end up at N model cases. These model cases are similar and a general 
study of model cases will be done in Section [TUl In particular, there we'll show the 
conclusions of Theorem 1 1.1 1 for each of the model cases. 
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In many reductions, we will use a variant of Varzugin's argument |31j . Our 
variant will be summarized in two perturbation schemes described in Section |4] 
These schemes are useful when the weights are in L 2 (in addition to being in L°°), 
and they are based on the Beals-Coifman operator formulation. Several technical 
estimates will be needed to make these schemes possible; they will be discussed in 
Section [SJ 

For technical reasons, a certain phase- weight relation in (w~ ,w + ) is essential to 
our reductions, and a conjugation by a scalar RHP can be used to achieve this rela- 
tion. This type of RHP was introduced by Deift and Zhou [5] for a slightly different 
purpose, and was also used in [31]. If M solves our initial RHP, 03 is the third 
Pauli matrix, and 8 is the solution to the above scalar RHP then the corresponding 
RHP for M s := MS~ a3 will have the conjugated jump matrix J con j = 8 a }J8+ a3 . 
A good choice of 8 will ensure the existence of a factorization of J C onj that gives a 
pair of weights (w~ , w + ) with the desired phase-weight relation. 

The first reduction will be localization to small neighborhoods of stationary points. 
We say a RHP with a pair of weight (w~ , w + ) is localized to a set if these weights 
are supported on this set. After this reduction, the weights will be supported on a 
small neighborhood of Ai, . . . , Xn (if doesn't have any stationary points, we can 
directly reduce them to (0, 0)). In this reduction the respective potentials u(t), v[t) 
differ by terms of high decaying order. 

The next reduction will be phase reduction. The given RHP now has already been 
localized to a small neighborhood of the stationary points. In each neighborhood, 
the phase 8 can be approximated by an analytic function (so overall we have a 
locally analytic function). Near a stationary point A of order k, we'll approximate 
9 by the following Taylor approximation: 

The third reduction will be separation of contributions. To separate the contri- 
bution of a stationary point A , consider the following setting: Let Wq be the part 
of w that is supported near Ao, and w\ — w — Wq the rest. Each pair of weights 
Wi gives rise to a new normalized L 2 -RHP. We'll show that unique solvability of 
the new RHPs imply that of the current RHP, and the current potentials u(t), v(t) 
differ the sum of the respective potentials uo(t) + ui(t),vo(t) + v\(t) by terms of 
high decaying order. 

The separation of contributions is nontrivial in general, as matrix multiplication 
is noncommutative. The main idea is to handle commutators that arise when 
commuting terms during the calculation. To separate the contribution of Ao, we'll 
use an a priori estimate involving the behavior of the solution to RHPs localized 
to small neighborhood of Ao- The a priori estimate will be proved using the same 
reductions, i.e. starting from now we will have to show that our reductions will not 
harm the validity of these estimates. 

For each stationary point, we'll generally carry out phase reduction before sepa- 
rating its contribution. The reason is technical: regarding regularity, the reduction 
of above-mentioned a priori estimate is more expensive than the reduction of u, v. 
Because two different stationary points might require two different regularity as- 
sumptions, we will iterate the phase reduction and contribution separation steps 
through the list of stationary points. 
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The final reduction will be a reduction to model case. To reduce a localized RHP 
with a nice phase Q(x) = a + bx k+1 to an appropriate model RHP, we will use a 
steepest descent argument that goes back to Deift and Zhou [5]. The model RHP 
will not have L 2 weights (the weights however remain in L°°). This lack of L 2 
integrability prevents direct applications of Beals-Coifman operator formulation, 
and hence our perturbation schemes. 

To overcome the above difficulty, the main idea is to reduce the current RHP to 
an approximation of the model RHP, which we'll refer to as the pre-model RHP, 
which has L 2 D L°° analytic weights. The sub-reduction to the pre-model RHP can 
be done in the same way as before. To reduce the pre-model RHP to the model 
RHP, we will follow jTO] and exploit analytic continuation to deform the contour M. 
of our RHPs to the complex plane, such that on the new contour, the pre-model and 
model weights have strong decay. This means L 2 integrability of the weights will 
be available and hence we can use our favorite Beals-Coifman operator formulation. 

The second difficulty in our reduction to the model RHP is the fact that the 
model weights are no longer localized to a small neighborhood of the current sta- 
tionary point. This means the same argument that was used to reduce the a priori 
estimates needed in our separation argument will not be directly applicable. The 
main idea is to exploit analyticity of the model weights (and the pre-model weights) 
to make up for this lack of localization. 

Finally, Theorem 11.11 will be verified for the explicit model RHPs that arise as 
consequences of the above reductions. For each stationary point, we have one model 
RHP. Beals-Coifman techniques can be used in the classical settings, i.e. when 
— 1 < pq < 1 at the current stationary point. The quadratic non-classical case can 
be studied using confluent hypergeometric functions as in |31) . and for the cubic 
non-classical case we'll use a result from isomonodromy deformation theory [21]. 

In every reduction, it will be clear that if the new (i.e. reduced) RHP is uniquely 
solvable, then so is the old RHP. Thus, the unique solvability of the model cases 
will automatically imply the unique solvability of the given oscillatory RHP. Indeed, 
we'll prove a norm estimate which implies the unique solvability of RHPs. The proof 
of this estimate uses the same reductions in the proof of Theorem 11.11 

3. Beals-Coifman operator formulation for RHPs 

Given two weight functions w + , w~ g L°°, consider the following bounded Beals- 
Coifman operator (acting on L p , 1 < p < oo): 

C w f = C+(fw-) + C-(fw+) 

Here C+ = H + , C_ = — H- and H + ,H^ are respectively the projections into the 
complex Hardy spaces H P (C+), H P (C-) of the upper and lower half planes. C w 
can take L°° input if the weights are also in LP . 

Usually, we'll be interested in L 2 setting where the weights are obtained from a 
nice factorization of J: 

J(A) = (I - w-)- 1 ^ + w+), w ± g L 2 n L°° 

The L 2 fl L°° condition will always be assumed on the weights in this paper, unless 
otherwise specified. 

Theorem 3.1 (Beals-Coifman). Suppose that w g L 2 n L°° are obtained from a 
factorization of J. If 3p g I + L 2 satisfying p = I + C w fi then the L 2 -RHP has a 
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solution: 

M±(A) =7 + C±(n(w+ + uT))(A) 
Proof. Notice that M + = fi(I + w+) and M_ = /Lt(7 - w~). □ 

This theorem only gives us existence. Notice that the functional equation of \x 
can be rewritten as: 

(6) (I - C w )((i - I) = C W I 

Thus, if (1 — C w ) is invertible on I? then ([6]) always has a solution on L? (it 
is not hard to see that the invertibility of C w is independent of the choice of the 
factorization). 

Under mild assumption on the factorizations, the invertibility of (1 — C w ) indeed 
implies the uniqueness of the solution constructed in the theorem of Beals and Coif- 
man. This was observed by Deift and Zhou |10j under the additional assumption 
that (I + w+) ±1 , (I - to")* 1 G L°°. 

The following identity (or its variant) will be used frequently in future compu- 
tation: If h, g, -\h{Hg) - \g{Hh) G Ui< P <oo LP th en 

(7) C(h)C(g) = C (- l -h{Hg)~ l -g{Hh)) 

(H denotes the Hilbert transform). This identity is a consequence of Privalov's 
theorem (see for instance Garnett |15j). 

An immediate consequence of ([7]) is that (under the same assumption on g, h) 
C±(h)C±(g) respectively belong to H P (C±), so vanish under C T thanks to the 
orthogonality of C+ , C_ . Note that this orthogonality relies on the fact that R 
with the obvious orientation is a complete contour |34] , see Section |9] for more 
details. Essentially, this means this consequence of (0 remains true if M is replaced 
by any complete contour E. 

As an example, we'll show the following important lemma in [9] (which was 
attributed to Zhou [34 ). Let 1 < p < oo and let C$,Cj respectively denote L p 
Beals-Coifman operators with L°° weights (— $+), (— $+) such that 

(J + $±)(/ + $±) = / 
Let w, w<s> be two pairs of L°° weights such that: 

(I+w+) = (I + «;+)(/+*+) 

(/-«£) = (I-w-)(I + $-) 

Lemma 3.2. If for some 1 < pi,P2 < oo we respectively have <fr± G H Pl (<C±) and 
$± 6 H P2 (C±), then the following identities are true on L p : 

(8) i-a». = (i - c«) o (i - c™) 

(9) l-O, = (1 - C $ ) o (1 - C W J 

(10) (1 - C») o (1 - C 5 ) = 1 = (1 - C 5 ) o (1 - C») 

Below, we'll include a proof of Lemma 13.21 for future reference (which relies on 
orthogonality of C+, C_). As we'll see, L 2 integrability of w (or won't be 
needed. 
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Proof. The proof of ([8|) relies on the following identities 

(11) (C±h)(x)$ ± (x)=C±((C + h)$ + - (CL/0*_)(aO 

for every h G L' P (M.), any 1 < p < oo. Indeed, using boundedness of $± we can 
always write $± = C±(g) for some g S L Pl with g,Hg G L°°. Now, using ([7]) 
it is not hard to see that (C + h)$ + G H p (C + ), (C_/i)$_ G H p (C_), which easily 
implies (flTl) . 

The rest of the proof of © is simply algebraic. Notice that (JSJ can be rewritten 
as C$ o Cm = C<j. + C w — C WiS ,. Using the definition of w$, it is not hard to see 
that the right-hand side of this equation is a Beals-Coifman operator with the pair 
of weights 

(— io~<I>_, -w + $ + ) 

So essentially we need to understand C$ oC w . Fix any / G L P (R). Observe that 
any Beals-Coifman operator C w can be rewritten in two ways 

C w (f) = fw- + C_(f(w-+w+)) 

= -fw+ + C+(f(w- +«,+)), 

note that f(w + + w~) G L p . Consequently, (| 1 1 j) and orthogonality of C± imply 

c+(CM)(-*-)) = -c + (fw-$.) 
c.(cm*+) = -c-(fw+<i> + ) 

which completes the proof of (|8]l. 

The proof of ([S]) is exactly the same with suitable changes of symbols. 

Finally, the first equality in (fTO)) can be viewed as a special case of (HJ when the 
weights w ± are the same as the weights ±$± (i.e. = 0). The second equality 
in (fT0|) follows by symmetry. □ 

In some cases, 4>± and $± may not belong to the above type of Hardy spaces, 
but they are still analytic and uniformly bounded on the respective half planes C± . 
In that case, the following lemma will be useful: 

Lemma 3.3. The conclusion of Lemma \3.2\ remains true if we assume $± and Q± 
are analytic and uniformly bounded on the respective half planes C± . 

Proof. The main idea is to directly prove (1111) instead of using ([7J (which requires 
the assumptions in Lemma 13.21 that we want to avoid) . The rest of the argument 
is exactly the same. To show (fTTj) . essentially we want to show that (C±h)&± G 
H P (C±) for /ieL p . Fix any z G C+. We'll show by contour integration that 

(12) (C+h){z)*+(z) - c((C + h)3> + ) (z) 

Choose the contour 7 to be a semicircle of radius R centered at the origin lying in 
the upper half plane. Then for large R, 

.„ xai , , 1 f R (C+h)(x)$+(x) , . 

(C+h)(z)$ + (z) = / dx + contribution of j R 

2ni J_r x — z 

where 7^ := {Re l/3 : < (3 < it}. Since p < 00, to show that the contribution of 
7ii vanishes as R — > 00 it suffices to show that 

\\(C + h)(z)^ + (z)\\ LPhR) <\\h\\ LP(m 
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This is a consequence of the boundedness of $+ in the upper half plane and Car- 
leson's measure theorem (see for instance [17])) which gives, for every 1 < p < oo 

\\(C + h)(z)\\ LHlR) < \\h\\ LPm 

Consequently, (|T2|) is proved for z G C+. Sending z — > K nontangentially, we see 
that (C+/i)4>+ G -ff p (C+) (observe that the boundedness of $+ on R ensures that 
(C+/i)$+ G LP). 

Similarly, (C_/i)$_ G H P (C-), and (fTTj) can now be deduced easily. □ 

Remarks: Lemma l3.2l and r3.3l remain true when R is replaced by a complete contour, 
and this fact is needed for contour deformation of model RHPs in Section |9] Using 
(TTTfl) , we can show a partial converse to the above observation of Deift and Zhou. 
The following proposition is useful when it is possible to construct a nice solution 
to ([6]). Indeed, we are able to do that for model RHPs associated with stationary 
points of primary or secondary order. The proposition then implies the unique 
solvability of these RHPs. 

Proposition 3.1. If (I + w+) ±1 and {I - w~) ±:L are in L°°(R), and has a 
solution [i G I + L 2 (R) such that both [i and its matrix inverse /i -1 are in i°°(R), 
then (1 — C w ) is invertible on L P (R) for any 1 < p < oo. Furthermore, 

H(i -a,)- 1 ^-^ < P iiMiiooii^ 1 iioo(ii(/+u' + )" 1 iioo + ii(/-^)" 1 iio ) 

Proof. Let M± = I + C±((w + + w~)fi) = (i(I ± w ± ) G L°° as in the theorem of 
Beals and Coifman. For any matrix valued / G IP we have: 

(1-C w )f = C+(/(l-«r))+C_(-/(l + «;+)) 
= C+(fn- l M-) -C_(/ At - 1 M+) 

= (c M °»- x )(f) 

here fj^ 1 act by right multiplications on IP (these are bounded since /i ±: G L°°); 
Cm and Cm- 1 denote Beals-Coifman operators with L°° weights (M_, — M + ) and 
(MZ 1 ,-M+ 1 ). 

Using Liouville's theorem, it is not hard to see that det(M) = 1. By ([TO]) . 

CmCm- 1 — Cm-^Cm = 1 
Consequently, (1 — Cyf)^ 1 = p Cm- 1 , and so its norm is controlled by 

ll(i -^n < ||/*||oo||c M -i|| 

< iHioodlM^iu + 1^1100) 

< Niocii^iioc^a+^r'iu + iKi-^-r^ioo) 

□ 

4. Localization schemes and essence of the nonlinear stationary 

phase method 

In this section, we'll describe the two schemes used to localize our oscillatory 
RHP (jl!2p to small neighborhood of the stationary points. In Section l4~4l we shall 
discuss potential adaptations of these schemes to other oscillatory Riemann-Hilbert 
settings. 
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Recall that the oscillatory jump matrix is given by ([2]), with the following canon- 
ical factorization: 

(13) J = (j Pe T){qe M °^=(I-^i)- 1 (I + 

fO pe~ ite \ + ( 0\ 

i.e. w x = ( q q J an( l w i = ( ue qJ- This factorization is not always 

suitable for our reductions due to the lack of certain phase-weight relation (a suit- 
able scalar RHP is needed to get the right factorization; Section|6]will be devoted to 
this discussion), but here we'll use it to demonstrate the ideas behind our schemes, 
avoiding unnecessary technicalities. 

Now, a perturbation of p or q leads to a weight perturbation by a term e h. 
Such oscillatory term will be decomposed into two components, which are handled 
separately using two schemes. We assume throughout that the weights are in L°° n 
L 2 . 

4.1. The direct scheme. The Beals-Coifman formulation suggests that small 
changes in the weights won't have any strong effect on the solution of the RHP 
and hence u(t), v(t). Let (Aw",Aw + ) be added to (w^ , w x ) to get a new pair of 
weights , u>2~)' The effect oii«,d can be quantified by the non-tangential limit 
limsup A ^ DO | A (Mi (A) — M2(A))| and will be estimated below. 

For any pair (w~ , w + ), for convenience we denote w := w + + w~ and 



IFIIp : = IF \\p+ \\w \\ P 
H p (w) := \\C^w + )\\ p +\\C + (w-)\\ p 

Lemma 4.1. Suppose thatwi is a pair of strictly-triangular weights and ||iyi||oo 1 
as t — > oo. 

(i) Suppose that ||(1 — C Wi )~ 1 \\L2^ L 2 < 1 as t — >■ oo for each i = 1,2. Then the 
effect on u(t),v(t) can be controlled by 

\\Aw\\ 2 (H 2 (w 2 ) + H 2 ( Wl )^) + \\Aw\\ 00 H2{w2)H 2 {w 1 ) 
(14) + | / Aw + wf\ + \ / Aw~w±\ 



+ limsup |AC(Au>)(A)| 

A— > oo 

Moreover, ||^ti — M2II2 ^ l|Aw|| 2 + || Aw\\ooH 2 (w 2 ). 

(ii) Suppose that ||(1 - C Wl ) _1 ||l 2 ^l 2 ^ 1 and ||(1 - C,^) -1 \\lp^lp < p 1 for 
2<p<ooast-^oo. Then the above estimates remain true if we replace 

\\Aw\\ 0O H 2 (w2) by \\Aw\\^H p (w 2 ). 

p— 2 



Notes: The third term in (JT4J) is if the matrix structure of w 2 resembles that of 
w\. In practice, some noise may affect the structure of w 2 so it won't necessarily 
be 0, but will be very small for large t. 

Proof, (i) Let A £ K. Write M 2 (A) - Mi (A) as 

C(fi 2 w 2 ) - Ciniim) = C((M2 - A»i)toi) + c((/i 2 - I)Awj + C(Aw) 

Note that if /eL'fl L 2 then limA^oo A(C/)(A) = — ^7 / f( x )^ x assuming non- 
tangential limit, and ||/^2 — -f||a = II (1 ~C W2 ) ~ 1 C W2 I\\ 2 < H 2 {w 2 ). Thus, by Cauchy- 
Schwartz, to show (fT4")) it remains to estimate J (fi 2 — /zi)wi. 
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Since /U2 - Mi = Caw^2 + C Wl (^2 - we get 

1^2 — Hi = (1 — C Wl ) 1 Cam)M2 = Ca«jM2 + C Wl (f) 

where <f> = (1 — C Wi )~ 1 CawI^2- Observe that for any sign combination 
||C±( M 2A W ; ± )|| 2 < ||A«;|| 2 + ||AH|oo||M2-/||2 < \\Aw\\ 2 + || AwlU^K). 

Thus, ||0||2 can be controlled by the same estimate. Since ||tui||oc ^ lj the desired 
estimate for fj,± — ^2 follows immediately. 

Now, using L 2 orthogonality of the two Hardy spaces H 2 (C + ), H 2 (C-), we can 
rewrite / CawH2 W\ as 

(15) J C + (f, 2 Aw-)(c.(w+) +^r) + J C.(fi 2 Aw + )(w+ + C+(ti;r)) 

Using Cauchy-Schwarz, the contributions of C_(u>^) in the first integral and of 
C+(tuf ) in the second integral are controlled by 

^\\Aw\\ 2 + \\Aw\\ qo H2(w2))h 2 (w 1 ) 

Below we estimate the contribution of w± in the first integral of (|15p. The contri- 
bution of Wi in the second integral can be estimated similarly. Now, 

J C + ((J, 2 AlV~)Wi = J ^Aw^W^ + J C-{p 2 Aw~)w^ 

Aw~w^ + J (fi2 -I)Aw~Wi + J C-(fi 2 Aw-)C+(wi) 

The last term can be controlled as before. The second term can be controlled using 
Cauchy-Schwarz and ||wi||oo J$ 1- The eventual estimate is 

Aw-w±\ + \\Aw\\ 2 H 2 {w2) + (\\Aw\\ 2 + \\Aw\\ oc H 2 {w2) S )H 2 (w 1 ) 

Similarly, J C Wl cj> W\ can be estimated by 

(\\Aw\\ 2 + \\Aw\\ 00 H 2 {w 2 ) S jH 2 {w l ) + \ J <Pw^w^\ + \ J <j)W+W+\ 

Using strict-triangularity of w%, we have w^w^ = w^w± = 0. This completes the 
proof of (i). 

(ii) The only difference is in the estimate of ||C±(/i2Aw; ± )||2, where we'll use 
Holder's inequality to control it by ||Aw||2 + \\Aw\\ 2 P \\(j,<>, — I\\ p < ||Aw|| 2 + 

||Atu||_2p_/fp(u;2). □ 

p-2 

Remarks: 1. The above argument also proves that — n 2 \\ 2 ^ ||CAu;-f||2 + 
I [ ^Xit? 1 1 00 1 1 A*2 — I\\2- This is useful if we want to directly exploit oscillation in Aw. 

2. For H p (w) to decay as t — > 00, we'll see that certain phase-weight relation in 
the pair (w~,w + ) is desired. 

3. The last term limsup A ^ oc | \C(Aw)(X)\ can be crudely controlled by ||Aw||i, 
but in our applications Aw may consist of Hardy components (which might not 
even be in L 1 ). In that case, we will exploit cancellation of these terms under the 
Cauchy transform. 

4. If ||Aw;||oo decays as t — > 00 then it suffices to assume uniform boundedness 
(as t — > 00) of ||(1 — C u , 1 ) _1 ||. For sharp applications we will try to avoid estimates 
involving ||Au;||oo because it generally requires stronger regularity assumptions on 
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w (hence p, q). This means if applicable (ii) is better than (i). For this reason, we'll 
often use weaker assumptions on the resolvent operators. For instance, in part (ii) 
we might change the assumption on (1 — C^) -1 to 

• (1 — C W2 ) _1 is bounded from L p+ n L p to L p uniformly as t — > oo, for 
2 < p < oo. 

The above argument then shows that 

IIM1-M2H2 < ||A«;||a + \\Aw\\_^_(^H q (w 2 ) + H q+ (w 2 )) 
and the first two terms in (fT¥|) should be adjusted to 

\\Aw\\ 2 (h 2 (w 2 ) + H 2+ (w 2 ) + H^wi)) + \\Aw\\j^(H q+ (w 2 )+H q {w 2 j)H 2 ( Wl ) 

and the eventual effect is only an e in the order of our decay estimates. 

4.2. The indirect scheme. The goal of the indirect scheme is to perturb the 
weights by large but very well-structured terms. This scheme is intuitively a reverse 
process of the direct scheme: instead of modifying the weights, we'll modify the 
Riemann-Hilbert factorization J = MZ from outside so that 

• The leading asymptotics of u(t),v(t) are not destroyed. 

• The effect on J is essentially an addition of these large terms to the respec- 
tive weights. 

Below, we demonstrate this scheme by showing how to perturb wf by $ + a 
large but well-structured term. Typically, $+ is an analytic function on C+ which 
vanishes as z — > 00 non-tangentially, furthermore its boundary value on K is in 
L 2 n L°° and oscillates as t — > 00. It is important that $ + continutes analytically 
to the upper half-plane if we want to perturb it from the positive weight wf . Now, 
consider the following normalized L 2 -Riemann-Hilbert factorization 

J = MZ 1 M+ 

where M + := M(I + $+) and J = J(I + $ + ). Typically, the matrix structure 
of $+ will be similar to that of w + (which in our applications will be essentially 
strictly- upper or strictly-lower triangular); after all we are only working with the 
only non-zero entry of this weight. This implies 

J = (I- wV ) _1 (/ + w+ + $ + + w+<I> + ) 
w (I -w7y\l + w+ + $+) 

so the essential effect on is an addition by 3>+. In practice (when we are not 
using the canonical factorization (JT3J) ) it is possible that the triangular structure 
of wi,$+ are affected by some "noise" terms (that decay strongly as t — > 00), 
but these noise terms are indeed small and can be easily handled by the direct 
perturbation scheme. 

Now, the effect of this modification on u(t),v(t) can be controlled by 

limsup \z(m+(z) - M + (zj) \ = limsup \z$ + (z)\ 



Using Cauchy theorem, the last limit can be written as 

1 



lim sup 



.1 f ^ld X 



TTI 



® + (x)dx 
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which decays as t — > oo thanks to oscillation of $+ on R. The last equation is a 
heuristics (since <fr+ may not be in L l ) but it captures the essential idea. 

4.3. Combination of two schemes. To perturb the weights by an oscillatory 
function h{x)e ltQ ^ x > (here = ±6*), following [3T] we'll decompose it using the 
Hardy decomposition: 

h(x)e it@ ^ = C+{he lt0 )(x) - C-{he lte )(x) 

Under mild assumptions on h's regularity, it turns out that if 0' keeps the same sign 
on supp(ft-) then one Hardy component is small for large t, and can be perturbed 
by the first scheme. The remaining term, consequently, inherits the oscillation of 
^ e ite fo u ^ now en j y S analytic continuation to one half-plane. It turns out that to 
ensure that this remaining term has continuation to the correct half plane (so that 
it can be perturbed by the second scheme), we want to have a certain phase- weight 
relation in the pair of weights (w~ 1 w + ). 

In Section [5] estimates on Hardy components of oscillatory functions are proved, 
which then explains the reason why a correct phase-weight relation is desirable. 

We make an important observation that in the second scheme, the solution (if 
uniquely exists) /i to the basic equation \i = 7+C,„/z stays unchanged, thanks to the 
formula fi = M + (I + u> + ) _1 = A/_(J — w~). In the first perturbation scheme, we 
have an estimate controlling ||A^|| 2 . Consequently, during the reductions, ||A/i|| 2 
always remains small. We'll see later that the above two schemes are analogous to 
the two main steps of the steepest descent method of Deift and Zhou: the second 
scheme is analogous to a contour deformation and the first scheme is analogous to 
a steepest descent estimate on the deformed contour. 

4.4. Essence of the nonlinear stationary phase method. In this subsection, 
we shall discuss further the underlying ideas of the nonlinear stationary phase 
method developed in this paper and how they may be adapted to other oscillatory 
Riemann-Hilbert settings, in particular those that have been successfully studied 
by the steepest descent methods of Deift, Venakides, Zhou. Successful adaptations 
of these ideas will allow us to reduce the analyticity assumptions to fairly weak 
regularity assumptions on the relevant Riemann-Hilbert data in these applications. 

To understand asymptotics of a given oscillatory RHP, the general strategy is to 
exploit the oscillation of oscillatory terms to show that the given RHP converges to a 
model RHP as the large parameter goes to oo, and the solution of the limiting model 
RHP will provide the desired asymptotics of the quantity of interests. For instance, 
in the AKNS setting, if the phase 8 has only one stationary point then the limiting 
RHP is a local model RHP associated with that stationary point, and when there 
are more than one stationary point we have to work a little harder to separate their 
contributions. In the small-dispersion KdV and orthogonal polynomial settings, the 
limiting RHPs may have a multi-interval configuration, more precisely the jump 
matrix will be locally constant |6] . The finiteness of the number of intervals is 
implied by real analyticity |6] or certain convexity of relevant Riemann-Hilbert 
data [lj and we shall assume this finiteness in subsequent discussions. For these 
multi-interval model RHPs, often the contribution of different intervals are not 
separated and the usual strategy is to construct an explicit solution for the multi- 
interval configuration [6l [7] . 

The steepest descent methods of Deift, Venakides, Zhou achieve the above con- 
vergence through exploiting the strong decay of the oscillating terms in appropriate 
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deformation contours. As briefly mentioned in the Introduction section, often the 
jump matrix of the given oscillatory RHP may not be ready for the steepest descent 
argument. The sense in which the jump matrix is ready is better explained through 
an example. Suppose on a subinterval / of the given Riemann-Hilbert contour 
(which we assume R for simplicity) the jump matrix J has a factorization 

(16) J(A,e) =v-(\,e)v(\,e)v + (\,e) 

where 1 /e is the large parameter and v^ 1 consists of oscillatory terms e ie ^ e that we 
want to remove. As part of the Deift-Venakides-Zhou steepest descent scheme, the 
contour of the given RHP will be deformed, moving v + to the upper half plane and 
v~ to the lower half plane where the oscillating terms decays strongly as 1/e — > oo. 
Consequently v + and v~ converge to constant matrices as the large parameter goes 
to oo, therefore if lim 1 / > :_ > . 00 w(A, e) also exists then we can determine exactly the 
limiting jump matrix on the current interval / (here we ignore all the endpoint 
issues where strong decay is not available and the use of certain parametrix might 
be needed). 

In the AKNS setting, an function 5 solving a scalar RHP is used to conjugate 
the given jump matrix 

and the scalar RHP is chosen so that the new jump matrix is ready for the steepest 
descent argument. In the small-dispersion KdV and orthogonal polynomial settings, 
this is achieved by more sophisticated versions of S often referred to as the g- 
functions (which goes back at least to [8]). 

The ^-function for a given oscillatory RHP is often determined by a reverse- 
engineering process, as follows. Suppose that e lS ^ x ^ e is an oscillatory term in J new 
that one plans to move from an interval J C I to a contour in the upper half 
plane. In other words in the factorization (fT6|) this oscillatory term would be part 
of v + (\, e). It follows from the Cauchy-Riemann equation that if 9 is (real) analytic 
then having 9'(X) > on / will ensure strong decay of e l6 ( x ^ e on the new contour as 
the large parameter 1 /e goes to oo (this observation goes back at least to [7] , see also 
the beginning of Section [9] of this paper for a philosophical discussion). Similarly, 
if one plans to move e iS / £ to the lower half plane then it is desirable that 8'(X) < 
on I. These conditions on the oscillating phases translate back to conditions on g 
(there are also other conditions on <?, for instance in the factorization (|16|) of J new 
we also want v(X, e) to converge as 1/e — » oo, but this is a different issue). 

As discussed in the last two subsections, it is essential for the two perturbation 
schemes used in this paper that the weights in our Beals-Coifman factorization 
have correct phase- weight relation. In Section El it will be shown that the correct 
phase-weight relation is 

(i) The oscillating phase in w + is locally increasing. 

(ii) The oscillating phase in w~ is locally decreasing. 

(for details see the remarks after the proof of Lemma 15. lj) . Writing J = (I — 
?«~)~ 1 (7 + w + ) = v~v + , one sees that the above phase-weight relation is exactly 
the monotonicity required of the phases to ensure the success of the steepest-descent 
argument of Deift and Zhou. One can look at this paper from the following angle: 
in the AKNS setting, given the setup of the steepest descent argument of Deift 
and Zhou, one can prove the desired convergence to model RHPs using the above 
two perturbation schemes. In other oscillatory RH settings, given correct phase 
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monotonicity in the correct factors of J n ew (which in turn is achieved by the use 
of ^-functions), one should be able to establish the desired convergence (to the 
model RHP) by adapting these two schemes to those settings. In other words, the 
only place where the method developed in this paper differs significantly from the 
Deift-Venakides-Zhou methodology is in the way one achieves the convergence to 
the model RHP, where one exploits the cancellation resulted from rapid oscillation 
of the oscillatory terms away from their stationary points/intervals by the use of 
real- variable tools. 

To apply the two perturbation schemes described in the previous two sections to 
other settings, certain adaptations might be required. In settings where the quan- 
tity of interests is recovered from the Riemann-Hilbert solution by a limit similar to 
(J3j) (e.g. small-dispersion KdV), the second perturbation scheme should be immedi- 
ately applicable provided that certain triangularity are available (which is often the 
case) . In settings where one is more interested in the asymptotics of certain entries 
of the Riemann-Hilbert solution (e.g. orthogonal polynomials), the relative position 
of the oscillatory factors and the entries of interests should be taken into account. 
For adaptations of the first perturbation scheme, one observes that this scheme 
is based on the philosophy that Riemann-Hilbert solutions depend "continuously" 
on the corresponding jump matrices, and in the AKNS setting this philosophy is 
quantified by the use of the Beals-Coifman operator for normalized L 2 -RHP. For 
applications to other settings, we may have to adapt the Beals-Coifman formu- 
lation to RHPs with different normalizations or different factorization structures. 
Note that in the implementation of the first scheme in this paper as previously 
described, triangularity of the weights has been extensively exploited to minimize 
the regularity assumptions on the Riemann-Hilbert data, but it is not essential. 

The next section provides the estimates needed for our perturbation schemes. 

5. Hardy decomposition of oscillatory functions 

Let k > 1. In this section, 9{x) will denote a real valued function such that: 

(i) 9 is (k — l)-time continuously differentiable everywhere. On the complement 
of a finite set of points, 9 has (fc + 1) locally integrable derivatives. 

(ii) 9 has stationary points Ai, . . . , Ajv of order kx, . . . , fcjv i.e. 9 is (kj + l)-time 
differentiable near Xj with 9'{Xj) = ■■■ = 9 {ki) (Xj) = ^ 9 (k ^ +1 '>(X j ), and 9^+^ 
is continuous at Xj. 

The value of k will depend on the actual proposition, but it should be clear from 
the context. In this section, we'll refer to (i,ii) as the two phase conditions. The 
assumptions (A), (B) for Theorem II .11 is for k = 2. 

We'll denote by x\ < ■ ■ ■ < x$ the points mentioned in (i). 

We'll need several notions of decay relative to 9. For any fc-time differentiable 
function /, let fe,k be the following weighted sum of its derivatives: 

EjLol/^l 

where w k ^(9) := J2 a , |6>'|-( fc + Ql +-+ a <=)|6K 2 )| Q i . . . \9^ k+ ^ \ ak 
here the sum is taken over all nonnegative integers a± , . . . , a k such that 
(17) a\ + 2a 2 H Vka k = k- /3 

(to avoid the 0° situation, the convention is: if some otj is then the respective 
power |«j i s treated as 1). The terms in the above sum appear naturally when 
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we perform integration by parts on certain oscillatory integrals. If asymptotically 
9 and its derivatives behave like a polynomial of degree d then Wk.p(9) < |x|^~ fc<i 
and we can estimate fe,k more explicitly. On the other hand, if / is compactly 
supported and supported away from the stationary points then 

k 

(is) h,k<\f(x)e^ k+1 \x)\ + 

/3=0 

For any d S R, let D(9, d, k) be the set of fc-time differentiable functions / 
satisfying the following two conditions: 

• \im x ^ 00 (e'y 1 f ej (x) = for every < j < k - 1 

• \9'\ d fe,k is integrable at infinity 

In particular, every compactly supported function satisfies the decay requirements 
D(0, d, k) for any 0, d. 

The next lemma extends Varzugin's localization principle in |31) . here some 
modifications are made to allow for less stringent assumptions on 9 and for a wider 
class of /. 

Lemma 5.1. Let 2 < p < oo. Suppose f G L 2 is k-time differentiable and supported 
in {9' > 0} s.t. /,/',..., /(k- 1 ) vanish on the boundary of {9' > 0}. If f has 
sufficient decay (say, D(9, 1 — -, k)) and vanishes at every endpoint \j 's of {9' > 0} 
with high multiplicity, then 

l|C-(/e !te )|| P < i-( fc - 1+1 /p) 

If k > 2, this estimate is also true for the endpoint case p = oo. 

Remarks: If supp(/) C {9' < 0}, the estimate is true for C+(/e ) by symmetry. 

Proof. We'll largely follow the argument in [3T]. Indeed, we'll show, V2 < p < oo 

(19) \CMe U6 )\\ P < pAk t- (fe - 1+1/p) (||I^T- 1/p /,, fc ||i + EEl/ W Wl) 

n=l P=0 

as long as kp' > 1. Basically, this means for p = oo we need k > 2, while for p < oo 
we only require > 1, and these are exactly the conditions on k described above. 
First, using Hausdorff- Young inequality, for every 2 < p < oo we have: 

\\C-(fe M )\\ LP(R) < ||7?n-OIL ? ' ([ o, oo) ) 

Indeed, this is the only place the condition 2 < p < oo (i.e. 1 < p' < 2) is really 
needed. The rest of this argument works for 1 < p < oo. This remark will be useful 
for future computation, in particular the proof of Lemma 15.81 

For f > 0, write = {lix)- 1 e i(x ^ +te ^\f(x)dx, which looks like a 

standard oscillatory integral, except for the term xt;. This term is however harmless 
because for x G supp(/), £ and 9'(x) are of the same sign and so: 

(20) \4-(xt + t0(x))\ = \£ + t0'(x)\ > t\9'\ > 

dx 

Thus, if / behaves nicely on the boundary of {9' > 0} (which are essentially sta- 
tionary points Xj of 9), we expect this integral to decay as t — > oo. 
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Denote w(x) — £ + t9'(x) (we suppress i,£ in for simplicity). Let D w be the 
following operator that acts on differentiable functions: 

w ^ dx^w(x)^' 
Integration by parts k times gives us: 

/>oo 

fe ite {-0 = boundary terms + i k / {D k w f){x)e^ x+lte{x) dx 

J — oo 

The boundary terms are evaluations of £iJl e < x + ltS ( x ) Dj,/ at various endpoints for 
some absolute constant C(j). For < j < k — 1 they are evaluated at ±oo and 
the endpoints of {9' > 0}. For j = k — 1 there are also evaluations at x\, . . . , xs 
as left/right limits, since is not differentiable at these points. We'll see that 

the first two types of evaluations vanish, while the last is 0(i _ ( fc_1+1 / p ') in L| . 

Indeed, the evaluations at endpoints of {0' > 0} vanish since / vanishes at Xj 
with high multiplicity. To show that the evaluations at ±oo vanish, we will control 
these terms by showing 

(21) \DU\ < t-*f 0d . 

This estimate and the given assumption on decay of / will then give the desired 
claim. To prove (|21[) . first notice that every (D 3 w f)(x) is a linear combination of 

(22) fW (x) • — i— w'(x) ai . . . w U) (x) a * 

with o-i > 0, n > 0. Now, our first observation is: £ no longer appears in the 
derivatives of w, indeed = t6^ +1 ^ . The second observation is 

n = j + cti H h Q'j 

This can be easily seen by a scaling symmetry argument: replace w(x) by c ■ w(x) 
and notice that D" w — -^D w , then do a counting of c. In addition, using a dilation 

symmetry argument (i.e. replace the pair (w(x), f(x)) by (w(x) := w(cx), fix) := 
f{cx)) and notice that D^,f(x) — cD w f(cx)) we can also show that: 

j — (3 + ot\ + 2a 2 + • • • + jctj 

Using ([20]) . we then can dominate (|22l) by 



i/ w i-^i^ (2) r---i^' +1) r 

. 1 lfl( 2 )| Ql |flO'+ 1 )| Q i 

- 1 \J I |^/|i+aiH haj \° I '•■I' 7 I 

Thus, from the definition of f$j we get the desired estimate (|2"TT) . The last task is to 
estimate the boundary terms coming from x n , which modulo an absolute constant 
are of the form 

lim -e^+^^D^fix), and lim -e^ x+lt0{x) D^" 1 f{x) 



vi w x ^ x - w 
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These limits are nontrivial only if 0'(x n ) > 0, in which case their (K+) norms 
can be controlled by the summation over {a.\ + ■ • • + a^-i + /3 = k — 1} of 

\te^{x n )\ a ^ ...\te^ +1 \x n )\ a ^ 



l/ (/3) M 



< l/ (/3) (^)|<" 1 + - +Qfc - 1 



i? ([0,oo)) 



fe+aiH hQfc-i-r7 



= r (i - 1+ %(»K)|. 

Note that the condition hp' > 1 guarantees convergence of the norm. 
By Minkowski's inequality, we can control \\fe lte {— £)ll £ p'([ ^ by 



'if ([0,oo)) 
S fc-1 



II / (^/)(^)^ te(a!) ^IL r([ o )00 )) 



1 1 boundary terms 

S k — i 

< t - (fc -i + i/p) J- J- |/W)(a:„)| + / \\(Dif)(x)\\ , dx 



n=l /3=0 

Using the triangle inequality, \\(Dwf)( x )\\ l p' ^ ^ can be controlled by: 



/3=0 on 



Suppose a; G supp(/) is not a stationary points (so 8'(x) > 0). Since hp 1 > 1, the 
LP norms in the last sum converge and we can rewrite this sum as 

± |/(«(,)| £ • • • |^)(x)|- ^|±g^^g^ 

,3=0 a; I wi 

s E E i« (s °(*r • • • i* (fe+1) (*)r ^p^g^pw 

= r( fc - 1+ ^|fl'(x)| 1 /p , / e>fc (x) 

Consequently, || /e iie (-£) II ,-*>', rn is controlled by 



i if ao.oo)) 

S k-1 



/oo 
_., -oo 

i0f /p 7^ 



n=l ,3=0 
S fe-1 



= *- (fe - 1+ ^(EEi/ w wi 

n=l /3=0 

giving the desired estimates p^|) . 

Clearly if / 6 D{9,^,k) then l^'l 1 ^'/^ is integrable at oo, so the L in- 
tegrability of this term depends on its behavior at the stationary points, where 
0' = 0. Below, we'll show that if / vanishes at every Xj with high multiplicity, then 
|0'|!/p j k i s integrable at Xj. 
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By a simple application of L'Hospital's rule, at Xj, 6^- n ' contributes a zero of 
order kj + 1 — n for every 1 < n < kj + 1: 

*», e( y +1 =r +x)(Ai Uo 

£— >Aj (x - A J )^-™+ 1 (fcj + 1 - n)! 
If fc < fc,-, the order at Xj of |0'|-(k+«i+-+a*)+Vp' |6K 2 ) | Ql • • • |6>( fe+1 ) \ a - is: 

k I k 

-kj(k + J2 ai -^)+J2^-^ ai = /3-(l + ^)^ + -7 
1=1 p 1=1 p 

using (fl7|) . If fc > fcj, the order at Xj of this term is at least that much, as 

kj k 

y^(fc 3 - i)on > ^(fcj - i)oti 

i=l i=l 

Since is locally integrable, fg^ is integrable at Aj when / "vanishes" 

there with high multiplicity. □ 

Remarks: The above lemma indicates that: to ensure the success of the perturba- 
tion schemes in Section 0] we want the following phase- weight relation: 

• The oscillating phase of w + is (locally) increasing. 

• The oscillating phase of w~ is (locally) decreasing. 

To see why, consider for example a perturbation of w + by /e lte , where O is the 
oscillating phase of w + . Decompose according to our schemes 

/e" e = C + (...)-C_(...) 

If 9 is locally increasing then C+(. ..) is the large and oscillatory term, and it 
has analytic continuation to the upper half-plane and can be perturbed from w + 
using the indirect scheme. If was locally decreasing then the large term will be 
C_(. . . ), and we will not be able to perturb it away from w + using our second 
scheme because it has analytic continuation to the wrong half-plane. 

The canonical factorization ([T3|) may lack this phase-weight relation, as 9' is 
not always positive on M. In Section HI we'll use a scalar RHP to get the right 
factorization. This scalar RHP is the source of the interaction between stationary 
points and the In t term in the asymptotics in Theorem 11.11 The introduction of 
this scalar RHP leads to the following type of functions: 

Definition 5.2 (A^ + decomposition). A measurable function f2(x) has an 
Ak + Bk decomposition (with respect to 0) if it is k-time differentiate and has the 
following properties: 

(i) Boundedness: \\Q\\oo ^$ 1 

(ii) Q^(x) < 1 + maxi^jxjv [^zxq^ for 1 < n < fc — 1 while fl^(x) can 
be decomposed as A^{x) + B^{x), with ||Afc|j 2 <fc 1 and \B^{x)\ <^ 1 + 
max. x <j< N \ x _\\ k 

(we understand that maxKjxjv ^ x _\ -\ k = ^ when 8 has no stationary point). 

A trivial example is f2 = 1. Observe that if f2 has the above properties then so 
does h£l, for any h G H k, °. Furthermore, if / e H k '° is compactly supported and 
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supported away from the stationary points then using local integrability of 
it is not hard to see that 

IK/nkikllx < II/IIjj*.o. 

We'll sometimes say that a function has an + Bk decomposition on a set. This 
means that (i) and (ii) are only required to be true inside that set. 

Definition 5.3 (Vanish with multiplicity). A function f vanishes at X with multi- 
plicity m > up to the fc th derivative on a neighborhood P of X if: 

(i) f is k-time differ entiable on P, and 

(ii) f {l) (x) <\x- X\ m ~ l for anyO<i<k- 1, and 
(m) \\(x-\) k - m -ifW(x)\\ LHP) <l. 

We'll say that / vanishes at A with multiplicity m up to the k th derivative if such 
P exists. Note that m is not necessarily an integer and (ii) is trivial if i > m. As 
an example, / £ H 1,0 has multiplicity up to the first derivative. If furthermore 
/(A) = then the multiplicity can be improved to 5, and if /' is also bounded near 
A then the multiplicity can be improved to 1 — e. 

Lemma 5.4. Let f be k-time differ entiable and supported near a stationary point 
Xo of order fco, where it vanishes with multiplicity m up to the fc th derivative. Let 
f2 have an A^ + decomposition on this neighborhood. 

(i) Let k < k(j. Then for 2 < p < 00, 

\\C-(X{0'>o } me lte )\\ p + \\C + ( X{ e'<o}f^e ue )\\ p 

here \nt can be removed if m + - 7^ (fco + l)(fc — 1 + i). The estimate remains true 
for p — 00 provided that m > and (m + ^) k( ^ +1 = * ,s re pl ace d by fco " 1 — e. 

(ii) Letk>k + 1. If9^ k+ ^ is U integrable near Xq for some 1 < r < 00 then the 
above estimates remain true after the following adjustments: fco + 1 is replaced by 
fco + 1 + c(k,p), and if p = 00 then m + jj = m is replaced by m — c(k,p). Here 

( 1 

(23) c(k,p) -- 

Proof. By translation invariant, we can assume that our stationary point is 0. For 
convenience, denote by {\x\ < 1} a small neighborhood of where / is supported. 

(i) It suffices to show the estimate for the first term on the left hand side. 
Non-endpoint case: Assume that 2 < p < 00. 

Let ip be a C°° cutoff taking values in [0, 1] such that <p(x) = for \x\ > 2 and 
<p(x) = 1 for |x| < 1. Let e > be a small scale. 

Let <p e (x) := I — ip(j ) and g := X{fl'>o}'/ 9 e/^- Using the triangle inequality and 
applying estimate (fT9|) to g, we can control \\C '-{x{0'>o} f^ eltS )\\p by: 

M-)m\\ p +\\C-{ge M )\\ p 
< e 1 ^ sup \f(x)\ +t-^- 1+ ^\\\e'\ 1 - x ^gg jk \\ 1 

|x|<e 

We'll optimize the estimate over e (by carefully keeping track of their orders). It 
turns out that a good choice of e decays as t —> 00, so the assumption "e is small" 
is harmless. Note that below p is allowed to be in [1, 00]. 
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For < (3 < k, expand g^> (x) using the product rule and estimate it by 

g^{x) < m f ) m ^\+ E i/ (n3 ^ (Q) ^i 

(n,a, j)eSp 

where Sp := {(n, a, 7) 6 1? + : n + a + 7 = /3, 7 > 0}. 

Note that supp(<^ € ) C {e < \x\ < 1}, and for 7 > we have i^ 7 ) < e — ^ l{ e <| a .|<2 e } ■ 
Thus, we can control by 

l(Rf) (/3) |l{e<|*|<i} + E l/ (n) " (Q) l^ 7 l{|,|- } 

The last part of the proof of Lemma 15.11 shows that near the 0-dependcnt 
weights of \g^\ in the weighted sum \6'\ ~pgg t k is controlled by 

, ,^-(i+fc )fe+io 

notice that the assumption k < ko implies that 6^ k+1 ' is bounded near 0. 

For \x\ ~ e the weight of g^ is of size e ^~( 1+ ' Co ' fe+ ~ j so \9'\ 1 ^ 1 / p ge^k is controlled 
by the summation over < (3 < k of l((3) + II(/3), where 

W = E l/^n^lW- 114 *^^^!.^} 



II(/3) = £|/WoW|e-V- (1+fe ° )fe+ #l 



{|x|~e} 

By assumption, 

(24) /W (#) < \ x \ m ~ n for < n < k - 1; 

and Ill^^^—i/W^II^ <1 

Contribution of I ((3): If a < k, n < k then we have a pointwise bound for QS a ) and 
/("). Thus, the contribution in L 1 of the associated term in I ((3) is controlled by: 

\x\ m ~ n ~ a+l3 ~ {1+ko)k+ 7~ dx 

£<\x\<l 

Since n + a — /?, the above integral is controlled by 

0(1)+e m+l-(l+Mfc+^ |mg|! 

where | In e | is not needed if to + 1 (1 + ko)k — -3, which is equivalent to 

(25) m + - ^ (As - 1 + -)(*o + 1) 

P P 

and O(l) is needed only when m+1 > (fco + l)A:— fy, in which case it corresponds to 
a positive power of the size of our neighborhood of 0. This means that we can make 
this O(l) smaller by making our neighborhood smaller. These comments apply to 
similar estimates in the future. 

Using Cauchy-Schwarz, the contribution of (n, a) = (k,0) can be estimated by 



|i { „i< 1 p fc - m -¥ fc) (*)i| 2 ( / ^ mH - {1+ka)k+ ^) r 

K J{e<\x\<l} 



which is again controlled by 0(1) + e' n+1 ( 1+fe °)' c + P > | \ n e | \ _ 
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To estimate the contribution of (n,a) = (0, k), we write = Ak + Bk and 
notice first that for Bk (which has nice pointwise bound) we can proceed as before. 
For Ak (with bounded L 2 norm), by Cauchy-Schwarz 

Ixl^l A k (x)\dx < M 1 

x\<l 

Using m + 1 < (1 + fco)fc — we can similarly control the contribution of Ak by 

l ,m+l-(l+fc )fc+ ± 9- r>H\ , ™+ 1 -( 1 + fc o) fc +^ 
SUp \X\ < O(l) + €j 

e<|x|<l 

Consequently, the contribution of Yl@=a ^0^) m ^ can ^ e controlled by: 

o(1) + em+ i-(i +fc0 )fc+^| lne | 

Contribution of 1 1 (ft). For any term in 1 1 (ft) clearly a < k and n < fc, so using 
(l24l) and pointwise bounds on derivatives of f2, we can estimate this term by: 

e -7+ l 8-(i+fc )fc+^ /" bl™-™-"^ « e m+i-(i+fe )fe+-^ 



Choice of e. Now, we can estimate \\C-(x{B'>o}f^e lt0 )\\p by 

e h m + i-C^+D (O(l) + e »»+Mi+*o)*+£ | lne |) 
_ i 

The optimal choice e = t k <~> +1 gives the desired estimate. 

Endpoint case: Consider p — oo under the extra assumption m > (this forces 
/(0) = 0). Note that ||C_(</e )||oo can be estimated as before, so the main task is 
to estimate: 

l|C-(x {e '>0}^(-)/^ te )l|oo 

Let h := X{e'>o} ( / 5 (j)/^ el * 6 '- By standard Sobolev estimates, we have 

l|C-C0lloo <, INIg+ll^llg VI < q < 2 
Since /(0) = 0, it is not hard to see that h is diffcrentiable, with the following 
derivative: 

x{0'>o } ( V(-)/« + ¥>R/'n + + ¥>(-)/fii0')e itfl 

\e e e ej e / 

Notice that #'(a:) is about the size of |a;| fe ° for x near 0. Since both h and h' are 
supported on {|x| < e}, using (|24|) we can control /i'(a;) by: 

i|x| m + {x]™- 1 + l^na;]" 1 + |:r| m i|a;| fco 
< \x\ m ~ l + te m+ka 

Using m > and choose q sufficiently close to 1, after integrating we get \\h'\\ q < 
e m_1 +« + te m+ko+ *. On the other hand, \\h\\ q can be easily controlled by e m+ ? . 
Consequently 

\\h\\ q + \\h'\\ q <e m - 1+ ^ +te m+ka+ ^ 

_ 1 

The choice e — t k o+ 1 gives an overall estimate 

l|C_(x {e <>o } ^ e ite )l|oo < q r (ro - 1+ t)wint + t"(*- 1 ) 
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now choosing q sufficiently close to 1 we get the desired estimate. 

(ii) The proof is similar. The estimate for the weight of g = will involve 
0( k + 1 ) (which is not assumed bounded) and this has to be estimated more carefully. 
Using Holder's inequality, for that contribution we can get an estimate of 

1 + e ' p | In e 

instead of 1 + e m+Mfeo+1)fc+ ^| lne|. Eventually we have an overall estimate of 
e -+| + t -(fc-l+i) (j + € m+^-(ko+l)k + ^ i ln ^ 

then optimizing over e will give us the desired estimate. The fact that we have p- 
instead of 1 leads to the adjustment of our decaying order, as compared to (i). 
The proof for the endpoint case is entirely similar. □ 

In the above estimates, notice that (modulo a positive power) the implicit con- 
stant for i^'^i+p) j s proportional to the size of the given neighborhood. 

Recall the notation fcg := max{0, ki, . . . , fc/v} and H p (w) := ||C+(u; - )|| p + 
1 1 C- ( w + ) 1 1 p . The following corollaries are consequences of Lemma lSTTl and Lemma [5Tl 
for(m,fc) = (0,l),(i,l). 

Corollary 5.5 (Rough estimate). Suppose w E IT 1 ' and both Qi, ^2 have Ai+B± 
decompositions, and Q\ , 82 satisfy the two phase conditions for k = 1. If w have 
sufficient decay at 00 , then for 2 < p < 00 

for the pair of weights tofi :— (^{e' <o}W~ ^2e lt62 , X{e[>o} w+ ^i elt9l ^j ■ If w± are 
compactly supported, the implicit constant is O p (^\\w + \\ H i,o + ||w~||#i.o^ 

Corollary 5.6 (Vanish at stationary points). Let f E H ' be supported near a 
stationary point \q of order kg, where fl has anAi+Bi decomposition. Iff(Xo) = 
then for 2 < p < 00 

\\C.(x { e'>o}f^e M )\\ p + \\C + ( X{ e' m me lW )\\ p < p max(r I , t~ ^ ^ +e ) 

Taking linear combination of Corollary \5Jj\ we see that generally H p (w) decays 
as t — > 00 if there is a correct phase- weight relation in w — (w~ , w + ) (now the phase 
is allowed to be piece-wisely defined). In that case, if ||(1 — C w )~ 1 \\lp^lp ^ 1 then 
the solution [i to fi = I + C w fj, satisfies: 

Ha* - Hp = ll(i - c^c^Wp < h p (w) < rn^+ir. 

Corollary 5.7 (Almost orthogonality) . Assume that flf, ttf have A\-\-B\ decom- 
positions and 0i, 62 satisfy the two phase conditions for k = 1. For each j E {1, 2}, 
consider the restriction to {0'j > 0} of a pair of compactly supported H ' weights 

w 3 = {x{0' ] >o} w J e ~ Ue ^X{e' ] >Q}wje lt9 =). 

If supp(u>i) and supp(u>2) are disjoint then V2 < p < 00, the LP operators C^ni > C W2 q 1 
are almost orthogonal in the sense: 

\\C Wini C W2 n 2 \\Lp^Lp < P t~**fr™\\w 1 \\ H i.o\\w 2 \\w,<> 
\\C W2 n 2 C Wini \\Lp^Lp < P r p(fc »2 +1) |hi||jfi.o||«;2||j !f i,o 
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for large t. The implicit constants depend on 9 and the support of the weights. 

Proof. By symmetry, it suffices to show the first estimate. Let / € LP . Since u>i, w-i 
have disjoint supports, distonce(supp(w 1 ), supp(w2)) ^ 1 and (C U , 2 /)(A) is analytic 
at every A € w\. By Holder's inequality, for every n > and A G supp(wi): 



d\ n 



(C W2 a 2 f)(\) 



< 



/(l) X W >Q } <(^)^ d[ 



(x-X) 

,, ,X{6'>o}W^{x)Q.^{x)e- m - 
f (x) x 2 ~ ' , tt-ti dx 



(x - A)"+ : 



< 



< 



ll/UKI 



x-A|>l |z-A|(«+Df 



Now applying Corollary [531 we have 



|C l u 1 n 1 (C! u , 2 n 2 /)||p < p i ^sTro f ||«?j l "C 1J , 2 n a /|| i? .i t0 + C W2 Q 2 f\\ Hli o) 



< 



t ^^^ll^llooll^lll^i-oll/Hp 



□ 



Again, by taking linear combination, Corollary 15.71 remains true if there is a 
correct phase- weight relation in the pairs of weights Wi and u>2- 

The next lemma summarizes standard results in linear theory (see for instance 
Stein [30 ) and is included here for convenience. In this lemma, / have enough 
decay to ensure that relevant tail L 1 norms are finite. 

Lemma 5.8. Let f be k-time differentiable and supported near a stationary point 
Ao where it vanishes with multiplicity m up to the k th derivative. Let ft have an 
Ak + Bk decomposition on this neighborhood. Then 

r ... [ max(^ fc ,t _ ^T lni), k < k 

(26) / fSle tte dx < { _ m +i 

J lmax(£- fc ,i fc ° +1 +A mi), k > k + 1 

Corollary 5.9. With the same assumptions as in Lemma \5.8[ for any Ao £ C such 
that distance(Ao, supp(f)) > 1 we can control C(/Oe lte )(A ) by the same estimate 
as in 



(write the left-hand side in integral form and apply the result of the lemma. 



Proof of Lemma \5.8[ Take a smooth cutoff ip as before, we'll proceed as in the proof 
of Lemma 15.41 By translation invariant we can assume that the given stationary 
point is 0. Let g(x) = (l — (p(-)}f(x)Sl(x). Using integration by parts, we have: 



fVte m dx < 



< 



e sup |/(o;)n(a;)| + 
esup \f(x)\+t- k \ 



\D%,g\dx 



Ml 



Then rest of the argument is exactly the same as before, applied to p — 1 (recall that 
the condition p > 2 was used in that proof only for Hausdorff- Young inequality). 
We then have the desired bound. □ 
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Recall that in Lemma 14.11 we have two error terms of the form J Aw wi , which 
vanish if there is agreement in the triangularity structures of w\ and W2- The 
introduction of the scalar RHP in the next section might create some noise in their 
triangularity structure, and the next two lemmas will be used to handle these noise. 
These lemmas are complex variants of Lemma 15.11 and Lemma 15.41 Below, T is a 
ray originating from one stationary point that forms a nontrivial angle with the 
real line, and fl has an Ak + Bk decomposition. 

Lemma 5.10. Let 1 < p < oo. // / has k locally integrable derivatives and is 
supported away from the stationary points of 6, and f has sufficient decay then 

\\C(fe M )\\ LP{r) < p r k 

Proof. Without loss of generality, suppose that L is originated from 0. Let g(x, A) = 
. Using integration by parts we have 



J 9 e%tS ^ |boundary terms| + t k J g s _ k 



It is not hard to see that gg ^ is controlled by a weighted sum of the derivatives of 
/, where for each < f3 < k f^' has the following weight 

J2 \x-\\-( n+ Vw ktj (6). 

n+0=j<k 

Consequently, by Minkowski's inequality we have 

k 

J xK p=0 J n+0=j<k 

which is finite if / has sufficient decay. Note that the condition p > 1 ensures that 
the L^(T) norms of \x — A|~(' i+1 - ) are finite for x £ supp(J). Also, similar to the 
proof of Lemma l5.ll the boundary terms of the integration by parts are either or 
can be controlled by the evaluations at Xi, . . . , xs of 

Since x\, . . . , xs are not stationary points, after taking L P (T) in A these boundary 
terms contribute 0(t~ k ). □ 

Lemma 5.11. Let f has k locally integrable derivatives and is supported near a 
stationary point Ao of order k$. Let f vanish at Ao with multiplicity m up to the 
k th derivative. 

(i) If k < ko then for 1 < p < oo, 

\\C(fne M )\\ LP(r) < p max(t- fc ,r (m+ ^Wlnt) 

(ii) If k > ko + 1 and is assumed L r integrable near Ao then the estimate 
remains true if ko + 1 is replaced by ko + 1 + ^ • 

Proof. By symmetry can assume that Ao = 0. Take a normalized cutoff function 
<p which is supported on {|ar| < 2} and equals to 1 on {|x| < 1} as usual. Then 
decompose / = (p(-)f + (1 — <p(-))f and estimate 

\\C{ V (-)fne M )\\ LP{T) < ||^(-)/^ te || LP(R) < e^ +m 
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For g := (1 — ip(-))f we proceed as in Lemma 15.101 and get an estimate of 

t - k e / mf) i0) \\^- n - 

n+P=j<k J 

As before, if k < ko then for every /3, the weight of (fl/)W can be estimated by 

|x|?-"-Vl J " (fco+1)fe = la^+^-^+D* 
(note that the estimate is independent of the choice of (n, j) as long as they respect 
n + (3 — j < k). Thus, the same argument as before gives the following estimate 
for the contribution in L 1 of those terms 

t- k {\ + ev +m - {ka+1)k )\\ne\ 
and consequently we will have an overall estimate of 

ep +t [1 + ep y ' | In e | ) 

1 

and the choice e = t k »+ 1 gives the desired estimate. 

When k > fc , the estimate involving 6*( fe+1 ' has to be done more carefully. This 
appears when /3 = n = and the contribution of the respective term in the overall 
estimate can be estimated by 

r fc (l + ep +m -^ (fco+1)fc |lne|) 

so again we have an overall estimate of 

£ p +m + r k {\ + e ?+ m -7-(fco+i)fc| lne Q 

1 

and the optimal choice is now e = t fc o +1+: F¥ . □ 

Note: The endpoint estimate when p = oo in Lemma 15.111 can be formulated 
and proved similarly, and is left as an exercise. Lemma 15.101 and 15.111 remain true 
if r C M provided that there is a nontrivial distance from supp(f) to F. 

6. FACTORIZATIONS OF THE OSCILLATORY JUMP MATRIX J(A, t) 

The operator formulation of Beals and Coifman indicates that a good factoriza- 
tion of J is crucial to the study of our oscillatory RHP. Intuitively, a good factor- 
ization should separate the two terms e ±ite , since they decay in opposite regions. 
From Lemma I5TT1 we want in our factorization the following behavior: 

• The oscillating phase in w~ is locally decreasing; 

• The oscillating phase in w + is locally increasing. 

In our applications, the phase is either 6 or —9 and could be piecewise-defined. 
The jump matrix ^ has the following canonical factorization: 

This factorization has the correct weight-phase behavior on {#' > 0} only, thus 
extra work is needed when both {6' > 0} and {8' < 0} are nonempty. Below, we 
will discuss an important scalar RHP which will be used to fix this problem. 

Let D- = {9' < 0},D + = {8' > 0} and consider the scalar L 2 -normalized RHP 
with the following jump matrix: 

(27) {l+pq)l D _+l D+ 
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The existence, boundedness and several important properties of the solution 8 of 
this RHP will be proved in this section. In particular, its behavior near the sta- 
tionary points will be studied. This 8 was introduced by Deift and Zhou [9] in their 
studies of mKdV and was also used in [31] . 

We describe how 8 can be used to conjugate the jump matrix J(x,t). Indeed, 



with the following jump matrix (it will be shown that 8 is bounded away from 0). 



The factorization (f28|) has the desired phase-weight relation. It is not hard to 
see that the conjugation doesn't change the unique solvability or the recovered 
potentials u(t),v(t). Furthermore, the resolvent operators (1 — C^)" 1 associated 
with these RHPs have comparable norms on IP . The last claim is a consequence 
of Lemma 13.31 and the boundedness of 8± . As we'll see, the above scalar RHP is 
the source of the In t term and the interactions of stationary points in the leading 
asymptotics of u{t), v(t) stated in Theorem ll.il 

6.1. Existence and boundedness. 

Proposition 6.1. If ln(l + pq) G L 2 (R) n L°°(R) is real-valued then the L 2 - 
normalized scalar RHP associated with the jump matrix J£7| ] has unique solution, 



Proof. Let 8 be as in the conclusion of the proposition. The jump relation is 
automatic. We only need to show that 8± — 1 belongs to the respective Hardy 
spaces i?2 of the upper and lower half planes (by symmetry only need to show this 
for 8 + ). This will then easily imply the uniqueness of the solution. 
For convenience, denote g = ln(l +Pq)1d- ■ 

Let A be in the upper half plane, write A = x + iy with y > 0. Since g is 
real- valued, 




(28) 



J conj {% j t) 







oo 
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Now, using the elementary inequality \e z — 1| < z|el Rc ( z ) and the above estimate 
for Re(C+g)(A) , we easily have: 

||5+(a: + iy)-l||2 < \\C+g(. + iy)\\ 2 e c|lsll ~ 

for some absolute constant c > 0. Thus, S±(X) — 1 are bounded and belong to 
the respective Hardy spaces. Since g is real-valued, 8(z)8(z) = 1 for any z £ 
C\D-. Thus, \S + (x)5-(x)\ = 1 on R. This plus the jump relation show the desired 
equalities for (5±(x)|. 

Below we show uniqueness. Suppose 8\(z) is another normalized solution to the 
scalar RHP ([27)1 . Notice that if we replace g by —g in the above calculation, we 
would get <5 _1 instead of S. So 5± x — 1 are bounded and belong to the respective 
Hardy spaces. Now the function m(A) := Si(X)6~ 1 (X) satisfies the jump relation: 

m+(x) = m—(x), x £ R 

Using boundedness of S" 1 and Si, it is not hard to check that m±(A) belongs to 
the respective Hardy spaces. From there it is clear that m(A) = 1. □ 

6.2. Bounds on derivatives of S. We will prove some results about the behavior 
of <5± on R. Part (b) of the next lemma generalizes a result of Deift and Zhou [TU] , 
where it was proved for 9{x) — x 2 . 

From now on we will define Qq(x) := <5+(x)<5_ (x). 

Lemma 6.1 (A-B Decomposition Lemma). Assume that ln(l + pq) G H k '° . Then, 
(a) fly 1 are k-time differentiable inside R \ dD_, and for 1 < n < k — 1: 

(29) i^o 1 )^) < 1 



da;™ 1 ~ distance(a;,aD_) n 

dx k 

\\A k h S 1, \B k {x)\ < k 1 



(b) -^pr(Q^' 1 )(x) can be decomposed into Ak(x) + Bk(x), with 



distance^, dD_) k 

Remarks: 1. If 9 doesn't have any stationary points then D- = or R. In these 
cases, our convention in the lemma is 9£>_ = and -; — ; = 0. 

aistance(x,oD- ) 

2. If ln(l + pq) is only H k Q on a neighborhood of xo (still L 2 elsewhere), the 
above properties remain valid locally. More precisely, we can use a smooth cutoff 
centered at xo to decompose ln(l+pg) into two components, one with high regularity 
and supported near xq, and another with low regularity but supported away from 
.xo. Then notice that the contribution of the low regularity part is bounded if 
\ dD_ is near xq . 

Proof, (a) Let g = ln(l +pq)lr>_. Since Q^^x) = e^^^^^ are bounded, to show 
(|2T)1) it suffices to show that (H g) (x) is fc-time differentiable on R \ 3D- , with 

d n 1 

(30) ^(Hg^x) < 1+— Vl<n<fc-1 

dx n distance(x,dD^) n 

If doesn't have any stationary points then g € H k, ° (and hence Hg g H k, °), and 
(|30| follows from a standard Sobolev estimate: 

ll/ (fe) ||oc<||/||^ +l ,o 
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In the rest of the argument, we'll assume that 9 has N > 1 stationary points. 

Below, we'll assume that 1 < n < k unless otherwise specified. 

Let <fio(y) be a nice C°° function supported in [— i, i] and equal to 1 in [— j, |]. 
Take M > be a small number and consider x £ {J a ed~\- a — M,a + M}. 

Define <p>i(y), (p 2 (y) by: 

Pi(j/) = Vo(-j^-), <P2 = (1 - Vi)lr>_ 

Notice that <p 2 is zero near the endpoints of £>_, thus it is also C°° . Decompose 
£>{Hg){x) into: 

in jn . . 

(31)%) +n(x) = — H(In(l + j*j) ¥ >il 2 >_)0»0 + ^r^( ln (! + 
Estimation of II: Since ln(l + pq) E i? fc, °, for all n < fc we have 
H(z) = ff([ln(l+ M ) V2 ]W)w 

= H[ln(l+p q )^^ 2 )(x) + f2( n )H(Hl+pq)^- j ^ { 2 j) )(x) 

3=1 

= ni(a;)+n 2 (x) 

For 112(2;): Notice that for any j > 1, we have ip^ < jgj and supp((/4^) C {y : 
distance(y, dD ) < M/2}. This means a; is outside the support of tp^ 1 furthermore 

distance(x,swpp((p2^)) > M. 
Thus, we can write every H^\n(l + pq)( n ~^tp2^ (x) in the integral form. Since 
|| ln(l +pq) ( - n ~ j) \\ oc < || ]n(l+pq)\\ H n,o, II 2 (x) can be controlled by 

n 

Il 2 (x) < \\ln(l+pq)\\ H n,oJ2j7 / \^\y)\dy 

i=i 7 

< || M i +M )|| fl . l0 (^ + ... + _l_) 

For IIi(x): Notice that ||IIi(x)||2 < || ln(l + pq)\\H n -°- Thus, even when n = k this 
term exists a.e. When n < k we can estimate it pointwise as follows: 

Ih(x) < ||ln(l+ M )("V 2 || H1 ,o 

< ||ln(l+ OT )||^ +1 ,o(l + ^) 

Estimation of I: Notice that x supp(<^i), thus 

1 f (-l)"n!ln(l+ W )(y) yi (y) 

it, ^ 1 r [Mi+w)yi] (nl (y) J 

= boundary terms H : / dy 



7Ti y - a: 

using repeated integration by parts. Here, the boundary terms are of the forms 
[lnq+MVtfV) 

(a; - a)"~» ' " ' " " " ' " G 
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Since </?i(a) = 1 and ip^\a) = for any j > 1, a € t?-D_, these boundary terms 
can be controlled by 

\\Hi+ Pq )\\H^ + --- + ^). 

On the other hand, the remaining term J D [ln(l +pq)ipi\ ^ (y)/(y — x)dy is easily 

controlled in L 2 by || ln(l + pq)\\H™-°, even for n — k. When n < k, it can be 
estimated a. e. by: 

IMi+p?) W) (y)ll^ w "%)l 



< 



|y - x\ 

2imi +pg )ii w+ x.o- / i^d/jidy 



< || ln(l +pq) 11^+1,0(1 + — + ... + —) 
Consequently, -£-^(Hg)(x) exists for every n < k, and ifl<n<fc — 1 then: 

£-(H 9 )( x )<\\Hi + P<i)\UAi + ^) 

for x \J a edD^i a ~M,a + M], Taking M = distance{x, dUZ)/2 we get (|30|) . 
(b) We can use the chain rule to compute -£^{^ 1 ){x) as 

J^ki^O 1 )^) = T^oO) [ (Hg) (x) + remaining terms J 
= A k + B k 

By a symmetry trick x h- > cx, it is not hard to see that the remaining terms on 
the right-hand side are linear combinations of those {H g)^ ll \x) . . . (HgY^^(x) with 
ii + ■ ■ ■ + ij = k and 1 < ii, . . . ,ij < k — 1. Therefore using part (a) we see that 
these terms are bounded by 

distance (x, dD-) k 

This gives the desired bounds for Ak, Bk- □ 

6.3. Approximation of 6. Compared to the last section, the added assumption 
in this section is p, q G H 10 , which makes ln(l + pq) € H 1,0 (M.). The goal of 
this section is to find an approximation for S at a stationary point Xj of 9. Our 
approximation will satisfy the following model scalar RHP: 

5 j+ (x) = 5^{x)(l Dj+ {x) + l Dj _ 0)[1 + p(X j )q(X j )]), x e M > 

where D j± = {x e K : ±0j-(x) > 0} 

Here 0j(a;) := 0{Xj) H fF+iT! ( x — ^i) ^ s tne Taylor approximation of 6{x) at 



Notice that we don't impose any normalization condition for Sj , only analyticity 
on C \ R is required. Intuitively, the above jump matrix can be seen as a limiting 
approximation of the jump matrix of S, which explains why such Sj plays an im- 
portant role in modeling RHPs localized to very small neighborhood of Xj (where 
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p(x),q(x) are essentially p(Xj), q(\j) and 9 is essentially 9j). Note that Dj± can be 
seen as the local approximation near Xj of D± = {±8' > 0}. 

A consequence of the result in this section is the existence of such Sj so that 
along any ray 7 originated from Aj that forms non trivial angle with R, we'll have 

\6(z)-6j(z)\<\z-Xj\K ze 7 

The upper bound can be improved to \z — Aj| 1_e if ln(l + pq) <E H 2 ' a , and the 
argument is similar with minor adaptations. 

The approximation of d is based on the approximation of its logarithm, C(l d- hi(l4 
pq)). Essentially, if h approximates this logarithm and Re{h) is bounded then by 
the mean value theorem 

S-cxp(h) < \C(l D -ln(l+pq))-h\. 

For simplicity of notation, we'll denote p(Xj) and q{Xj) by pj and qj. Consider 
the following triangular function (supported on [—1,1]): 

if \x\ > 1; 
T(x) = { x + 1, if -1 < x < 0; 

if < x < 1. 

(note that T e if 1 ' , this is also the reason for choosing it triangular). Let c > 1 
be such that c < min^ |Ai — Xj \ and let T c (x) :— T((x — Xj)/c). 

To get a sense of the desired approximation, we'll compute C(1d_ \n(l+pjqj)T c )(z) 
for z G C \ R. Thanks to the constraint on c, we can write it as 

H 1 +PjQj) f T ^ x \ dx 




2ni 



Jd^ x - 



To compute this integral explicitly, we have to examine Dj_ 7 or equivalently the 
relative position of Xj within Let e be a function on R that assigns 1 to every 
right endpoint, —1 to every left endpoint, and to every interior /exterior point of 
D_ (note that endpoints of D_ are stationary points). Note that for every j, e(Xj) 
depends only on the parity of kj and the sign of 0( fe j +1 ) (Xj): 



0, if kj is even; 

sgn (6»(fc 3 +i)(A J )), iffcjisodd. 



Recall that Vj := — ^ln[l + Pjqj]. Consider the first case when e(Xj) = 1. Then 
the above integral can be rewritten as 



■ / X - X 3+ C i 

iv j I — ; — - — — ax 

JXj-c c(x-z) 

iv . c+(z-Xj+ c)[ln(z -Xj)- Hz - Xj + c)} 
3 c 
= ivj ln(z — Xj) + ivj g c (z) 

where g c (z) := 1 + (^{z — Xj) ln(z — Xj) — (z — Xj + c) ln(z — Xj + c)j/c. It i 
hard to see that 

11^(^11^(7) ^ 1 



is not 
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where the implicit constant depends on c and the angle between 7 and R (indeed 
near Xj, g' c is < | ln(z — Xj)\ which is L p -integrable there, while for large \z — Xj\ 
we can estimate g' c by ^— x~|)- Consequently by Cauchy-Schwarz's inequality, 

\g(zi) - g(z 2 )\ < \zi - z 2 \ 1/2 , Vz 1 ,z 2 £ 7 
so 3^.m. z ^\^ zei g(z), which indeed exists nontangentially (i.e. independent of 7): 

lim g(z) = 1 — lnc 

On the other hand, it is not hard to see that exp(ivj ln(z — Xj)) satisfies the de- 
sired model jump relation. Thus, the desired approximation for exp (c(1d- hi(l + 

Pjqj)T c )^j in this case is: 

exp {iuj ln(z — Xj) + icej 
for some real number a that can be computed by the following nontangential limit: 
a = - lim (c[1d_ hi(l + pjqj)T c ](z) — ivj ln(z — Xj)) 

When (.(Xj) = — 1, the above argument essentially works, except for one place: 
in the computation, we need to choose ln(x — z) as the anti-derivative of — 
instead of ln(z — x). This will ensure that the respective g c (z) (which is now 
1 + ^[z — Xj) ln(z — Xj) — (z — Xj + c) ln(z — Xj + c)^/c) has nontangential limit as 

z — > Xj independent of the approaching direction (more precisely, of Im(z — Xj)). 
Consequently we will have a slightly different approximation 

exp ^ — ivj ln(Aj — z) + iaj , a is defined by a similar limit. 

When e(Xj) = (i.e. kj is even), there are two possible scenarios: Dj- = 
(which corresponds to 9^ kj+1 ^(Xj) > or cquivalcntly Xj is an exterior point of 
£>_), or Dj- = M. (the opposite case). In the former case, clearly C(1d_ ln(l + 
Pjqj)T c )(z) = 0. In the latter case, we can "sum" the above approximations to get 
the desired approximation: 

exp {ivj ln(z — Xj) — ivj ln(Aj — z) + ic^j , a is defined by a similar limit. 

This can be further simplified using ln(^ — Xj) — ln^ — z) = 7risgn[Im(z)]. Con- 
sequently, the desired approximation of exp ((3[1d_ ln(l + PjQj)T c ]^ is of the form 
exp(ia + (3j(z)), where 

0, if Xj is an exterior point of Z?_; 

PjW = \ * e (Aj)^j In [e(Xj)(X — Xj)] , if Xj is an endpoint of D_; 

-7rz/jSgn[Im(A)]], if Xj is an interior point of _D_. 



= - lim (c[l D _]n(l+ P jqj)T c }(z)-pj(z)) 



In the next proposition, we'll show that the desired approximation of C(1d_ ln(l+ 
pq))(z) is of similar form, 

5j(z) := exp(iu!j + Pj{z)), 
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where uij is defined by the following nontangcntial limit: 

(32) Uj := - lim (c[l D _ln(l + pq)](z) - pj(z)) 

Remarks: Clearly, if the limit exists then uij G R. From there we can easily see that 
Sj is bounded on C \ R. Indeed, the only nontrivial case is when Xj is an endpoint 
of D-. In that case, VA G C, 

\Sj(X)\ = exp ( - e{Xj)vj arg[e(A j )(A - A,)]) < exp(|| In(l +pq)\\ O0 ) 

Proposition 6.2. Assume p, q, ln(l + pq) G H 10 and ln(l + pq) is real valued. 
Then: 

(i) The limit defining uij exists nontangentially. 

(ii) V ray 7 originated at Xj that forms an angle of measure > 1 with R: 

|5(A) -^(A)| < |A- A^l* 

Proof. Let h = ln(l +pq) — T c . It is good enough to show that \\m z ^ r \ j C(l _d_ h)(z) 
exists nontangentially, and on any such 7, C(1d_/i)(z) is Holder continuous of 
exponent 1/2. Indeed, we'll show that these claims are true for every h G H 1 -°(R) 
such that h(Xj) = 0. 

Decompose 1d_ h as hi + h 2 = lDj-h + l D _\ Dj _h and notice that hi is in H 1 ' 
and vanishes at Xj, while h 2 is in L 2 C\L°° and supported away from Xj. We'll show 
the claims for hi and h 2 using these properties. 

First, since hi G H 1 ' ( so is continuous), the nontangcntial limits lim z ^>,j {C±hi)(z) 
exist and equal to (C±hi)(Xj). Since (C + hi)(Xj) — (C_hi)(Xj) = hi(Xj) — 0, these 
limits are the same, so lim z ^\ j (Chi)(z) exists as desired. Furthermore, by a simple 
application of Carleson's measure theorem, we have 

||(C7n)'|| i2(7) = 11(0-^)11^(7) < WKWlhr) 
so by Cauchy-Schwarz (Chi) is Holder continuous of exponent 1/2 on 7. 

For h 2 , the fact that h 2 is supported away from Xj easily implies the existence 
of the limit, indeed (Ch 2 ) is now analytic at Xj. Again, for the estimate it is good 
enough to show that (Ch 2 )' G L 2 . Using the support condition it is not hard to see 
that for z G 7, 



(ck 2 y(z) = ^[j^- 2 d X 

2m J R (x - zf 



< \\h 



2||oo: 



l + \z-Xj\ 

notice the implicit constant depends on the angle between 7 and R. Consequently, 

||(C7, 2 )'(z)|| i2(7) < Ufclloo 
as desired. □ 
Remarks: When ln(l + pq) G H 1 ' 1 , we can compute ujj more explicitly by 

( 33 )_L f ]n|A J --y|dIn[l+pg](y)+ ^ e(A fc )i/(A fc )ln|A J - - A fc | 

^ •' D - l<fe<AT,A fc #A 3 

using suitable integration by parts. The condition ln(l + pq) G H > in particular 
ensures that the infinity boundary terms of the partial integration vanish. To see 
this it might be convenient to distinguish between values of e(Aj). 
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7. Unique solvability of oscillatory RHPs 

7.1. Unique solvability of the original RHP. In this section, we sketch the 
main ideas to show that for t large enough, the RHP ([TJ [2J is uniquely solvable. 
In light of the discussion in Section [3l it suffices to show that for large t there 
is a factorization of J(A,<) = (J — w~(\ ) t))~ 1 (I + w + (\,t)) such that I — C w is 
invertiblc. 

Theorem 7.1. Let 9 satisfies (A ) and (B) andp,q are continuous and vanish at oo, 
with < 1 + pq < 1. Assume that pq < 1 at stationary points of order > 3. TTiera 
i/iere exists a factorization of J(t) such that — C w rt)) \\L r -+L r ;$ 1 uniformly 
as t — > oo /or r > 2 sufficiently close to 2. I/fce < 2 and pq >0 at every secondary 
stationary point then we allow 2 < r < oo. 

Along the reductions used in the proof of Theorem 11.11 Theorem 17.11 can be 
proved for a nicer class of p, q, say p, q G C^; furthermore for large t the resolvent 
norms || (1 — C lu ) _1 1| are bounded by constants depending essentially on the values of 
\p\ and \q\ at stationary points. To extend this to more general p, q, just approximate 
them by Pi,qi G such that 

• < l+pxqi < 1; 

• (p,q) agrees with (pi,<Zi) at every stationary point; 

• p — pi and q — qi are uniformly small on M; 
and invoke an application of Newman series. 

Let p,q G Cg now. Conjugate J using the solution S to our scalar RHP as in the 
last section to get J con j; this does not affect the conclusion of Theorem 17.11 thanks 
to Lemma 13.31 After factorizing the jump matrix J con j{^,t) nicely as in Section [51 
we use Lemma 15.11 described in Section 2] to reduce the corresponding weights to 
a sufficiently small neighborhood of stationary points. In this reduction, we will 
show that for large t and for any 2 < p < oo: 

11(1 - Cv^Wl^l, < p 1 + 11(1 - C^J-Il^l, 

Near each stationary point Ao, 9 will be reduced to an analytic phase, after that 
the contribution of Ao is separated. Eventually we arrive at a number of simpler 
RHPs, one localized to a small neighborhood of each stationary point with a nice 
analytic phase. Each such RHP will be further reduced to a model RHP. The 
validity of Theorem 1 7. II for the model RHPs will be shown separately in Section [TOl 

For the separation of contribution, let w = Wq + Wi where wq = (wq,Wq') 
is supported near Ao and w± — (wY,Wi) is supported near the other stationary 
points. The main idea is to use the following paramctrix, introduced by Varzugin 



(i-c^ii + Y^c^ii-c^)- 1 ) = i- c Wj c Wn (i - c^y 1 

j=Q 0<j^n<l 

i 

(i+^c^i-c™,)" 1 )^-^) = i- i 1 - c w n y l c Wj c Wn 

j=0 0<j=jtn<l 

Thanks to the scalar RHP |2"T1 which places the oscillating terms e ±zt61 in the right 
places, we can apply Corollary 15.71 to see that \\C Wn C Wj \\lp^lp's decay as t — > oo 
(whenever n ^ j). This enables the success of the above parametrix. 
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7.2. Unique solvability of intermediate RHPs. In our reductions, we want to 
ensure that: 



Corollary 7.2. If p,q £ H ,Q then for large t the intermediate RHPs that appear 
as consequences of the reductions are uniquely solvable, with uniformly bounded 
resolvent norms (as t — >• oo) in relevant L p . 

Proof. Thanks to the above parametrix, it suffices to show the resolvent bound for 
an RHP associated with a pair of weights (w~ , w + ) localized to a neighborhood P 



of one stationary point. Let 
on P. Then we can assume 



be a normalized smooth cutoff function supported 



(w , w + ) = < 



'0 <j)8-5 + pe- lte \ I (P 

,0 J ' ^xT^+V 6 0, 

0\ (0 <j)5-5+pe~ lte 

,0 



o, 



if A £ 
if A £ £>_ 



Now, approximate p, q by p\ , gi £ such that — pi||iyi,o and ||q — 9i||/fi,o are 
small, furthermore p and p\ agree at the current stationary point, and so do q and 
q\. For 5 = e C<,lD - ln ( 1+P9 ^ ; we'll approximate ln(l + pq) by h £ real valued, 
such that they agree at every stationary point and || ln(l +pq) — h\\ H i,o is small. It 
is then not hard to see that Si := e c ^ lD - h ^ is close to 5 in L°°. Consequently, the 
approximated weights (W~~ , W + ) (where p,q,8 are replaced by Pi,Qi,Si) is close 
to (w~ , w + ) in L°° . 

Now, note that Sli = <5i_<5i + and its inverse have ^3 + B3 decompositions on P. 
The same reductions as in the proof of Theorem 17.11 can be used to show that the 
RHP associated with W satisfies the desired resolvent bound, and for large t the 
bound depends essentially on the value of [pi|, \q±\ at the current stationary point. 
Consequently if the above approximation errors are small enough then we get a 
comparable resolvent bound for (w~,w + ). □ 



8. Reduction to model cases (I): Localization and phase reduction 



In this section and the following section, we will use the localization schemes 
described in Section 2] to reduce the RHP (|1I2|1 to a number of model RHPs, at 
the expense of modifying the solution and the potentials u(t), v(t) by terms having 
sufficient decay as t — > 00. 

As discussed in Section [SI we first conjugate our RHP (|1I2[) by the scalar RHP 
(|27| and end up with the RHP (M s , J con j). For simplicity of notation, let P = 



n p = 8-5+p, Q 



n~ l q : 



J_ 1 (5 + 1 (7. From (p?gj) . the weights for J con j are: 



(w , w + ) = < 



'0 Pe- M \ ( tr 
v J ' \ Qe itB 
0\ fo Pe- iW 
Oe ite O/'lo 



if x £ D+: 
if x £ L>_ 
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It is clear that P,Q,]n(l + PQ) are still in£ 2 nL°°. Indeed, ln(l+PQ) = \n(l+pq). 
The reduction order is: 



The middle two reductions are iterated through the list of the stationary points. 
In general, the order of these two reductions can be interchanged; we choose the 
above order to avoid undesired regularity assumption on 9. Part I (this section) 
describes of the first three reductions. 

For simplicity, we'll assume that p,q £ H k '° with k > 2, but overall it will be 
clear from the argument that k > 2 is not always required. In particular, k = 1 is 
always allowed in the separation of contributions. 

We'll always assume that p and q and their relevant derivatives have sufficient 
decay at oo. 

8.1. Localization to neighborhood of stationary points. The argument used 
here will be used as a model for future reductions, and we will often refer to it 
whenever the needed proofs are essentially similar. 

Let 1 = 4> + + </)_+ 0o be a C°° partition of unity, where 

(a) < 0+,</>_,0o < 1, 

(b) <t>- are respectively supported in D +1 D_, and 

(c) supp(</>o) consists of small disjoint neighborhoods of Ai, . . . , \n, and 4>o(x) = 
1 for x sufficiently close to any stationary point. 

Our goal in this subsection is to reduce P, Q to P4*o,Q4>o- The new weights are 
denoted by wf and the new solution is denoted by Ml- 

In the schemes described in Section [H we decompose Aw := w — wl using 
the Hardy decomposition. Recall the notation H p (w) — \\C + (w~)\\ p + \\C-(w + )\\ p 
for a pair w = (w~,w + ). Thanks to correct phase- weight relation, the following 
estimates are true for 2 < p < oo: 



When k > 2, p — oo is allowed in these L°° decays will be useful for the 

reduction of unique solvability and resolvent bound (i.e. for Theorem 17. ip . On 
the other hand, if kg < 2 and there are no defocusing stationary point then the 
reduction of u(t), v(t) in this section can be done without requiring these L°° 
estimates, consequently k — 1 is allowed. 

For simplicity of notation, let weights w K be defined by 




(34) 



H p (l D _Aw),H p (l D+ Aw) < t^ k - 1+ ^ 



w 



w 



C-(Aw + ), w 



K 



w 



L 



C+(Aw-) 



(I3"4l implies that for 2 < p < oo 



\\W K -W L \\ P < t-( k - 1+ ^ 
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so in principle we can go from wk to wl using the first scheme. To go from w to 
wk, in the spirit of the second scheme we introduce M* = M s (/+$), where $± are 
respectively dehned on the upper half and lower half planes of C \ K. Heuristically, 

($)+ « -C+(Aw+) = -C + (l D+ Aw + )-C+(l d_A w +) = (*i)+ + (* 2 )+ 

(*)_ « -C_(Aio-) = -C , _(l D+ Aw;-)-C7_(lfl_Aw-) = ($i)_ + ($2)- 

The triangularity of (($i)_, ($i)+) and (Aw", Aw + ) agree on D + and disagree 
on Z3_ (opposite story for $2)- Thus, at a time only one term on the right-hand 
side has correct triangularity needed for the second scheme. It turns out that, as 
t — > 00, the other term is always small: ( < &i)± are small on £>_ and ( < &2)± are 
small on D + . This observation indicates that, modulo small "noise", the matrix 
structure of <&i and $2 are good enough for our applications. 

The simple choice $ = $1 + $2 turns out to be inconvenient since it is not 
guaranteed that I + $ is invertible on C \ R. Although when k > 2 this would not 
be a problem since it can be shown that $ is asymptotically strictly-triangular, it 
is better to avoid this restriction. A better "superposition" of <&! , <J> 2 is: 

(35) $ = $i+$2 + $2*l 

(this has the algebraic advantage that det(J + $) = det(J + $2) det(7 + 4>\) = 1 
since <I>i are strictly triangular.) This choice of <f> goes back to [31] but can be 
naturally interpreted as successive applications of the second scheme by $2 and $1. 
The jump matrix for M* is: 

J* = (I + <S>-)- 1 J(I + <P + ) = (I-w^)- 1 (I + w+) 

here the weights u>$ are defined via 

(I + w+) = (I + w+)(I + <£+) 

(I-wi) = (/ -«;-)(/ + $_) 
We can compute them explicitly below: 

w$ = I — (I — w~)(I + <&_) = w~ — $_ + 

w t = (I + + $+) - I = w + + $ + + 

Because of the complication arose from lack of uniform triangularity, u>$ differs wk 
by some noise, and the sequence of subreductions will be w — > w$ — > wk — > wl- 
For convenience we group the last two into Proposition 18.21 In the following two 
propositions, part (i) should always be viewed as the reduction of Theorem 17. II for 
nice p(x), q(x), while the other parts are for the reduction of u(t), v(t) assuming the 
validity of Theorem 17.11 and Corollary 17.21 

Proposition 8.1 (w — > w$). (i) If k > 2 then for any 1 < p < 00 the following 
estimates are equivalent: 

||(1 - C'uj^Wlp^lp = O p (l) as t -> 00 

||(1 - CUrii,-^ = O p (l) as t -> 00 
(ii) Assume that k > 1 and both resolvent operators exist on L 2 (although their 
norms may not be uniformly bounded as t — > 00). Then 

limsup \X(M 8 (\) - Af*(A))| < r k 

X— >OG 
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Proposition 8.2 (u>$ — > wjj). (i) If k > 2 then for any 1 < p < oo the following 
estimates are equivalent: 

||(1 - C WL y 1 \\ LP ^ LP = O p {l) as t -> oo 

||(1 - C W9 y l \\ LP ^ LP = 0„(1) os t -»• oo 

(»,) Suppose that as t — > oo oot/i (1 — Cu^)" 1 and (1 — C lu ) _1 are uniformly 
bounded on L? . Then for k > 2 we have 

limsup |A(M*(A) - Mi(A))| < e r ( ^ 1+ ^rr )+e 

fni) Suppose that as t — > oo both (1 — C,^) -1 and (1 — Cu,) -1 are uniformly 
bounded on L p for every 2 < p < oo. T/ien /or k > 1 we have 

limsup |A(A/* (A) - Mi(A))| < £ r (k '^ + ^P^) )+e 

A— >oo 

Recall that by the A-B decomposition lemma, fij^ 1 have Ak + B^ decompositions. 
Proof of Provosition \8.1\ (i) Let $ = -$i - $ 2 + < J > i'J > 2, then 
(7 + $)(/ + $) = (/+ $)(/ + $) = J 

The desired claim is now a consequence of and (JTUJ). Here, note that 

the L°°(R) norms of $± and $± can be controlled by ||($i)±|| 00 + || ($2)± ||oo + 
IK < J ) i)±l|oo||( < i > 2)±l|oo which are bounded as t — >• oo (for fc > 2). 
(ii) Let k > 1. Oscillation of Aw ensures that 

limsup|AC(l jD± Au) ± )(A)| = — | / Aw ± (x)dx| < i~ fe 

A^oo 27T J 

The above nontangential limits converge thanks to integrability at oo of p, q. This 
is the only place where this requirement is used and future implicit constants will 
not depend on these two L 1 norms. The above estimates imply 

Urn |A^-(A)| < r k 

A— >oc 

Consequently, 

limsup |A(M 5 (A) - M*(A))| < limsup |A*(A)| < r k 

A— >oo A— >oo 

here recall that lim^oo M 5 (X) = I, □ 
Proof of Provosition \8. c 2[ (i) We'll show that, for k > 1 and 2 < q < oo, 

(36) \\w*-w L \\ q <t- ik - 1+ ^ +e , 

and if k > 2 then q = oo is allowed. This will automatically imply (i). We'll show 
the estimate for [|1d + (ix;£ — w^)|| p and the argument for other sign combinations 
are entirely similar. Write 

1d+(w$ - u>t) = ln+(wj - w + ) + l D+ (w + - iu£) 

(37) = l D+ (3> + +w + $ + ) + l D+ Aw + 

Notice that as t — > oo, w + is uniformly bounded, ($1)+ is uniformly bounded 
in L°°~ , and lz? + ($ 2 )+ = -1 D+ C-(1 D _ Aw+) decays like Oft - ^ -1 " 1 "^) in IP. 
Consequently by Holder's inequality, ($i)+($2)+ decays like 0(< _ ^ _1+ p^ +e ) in 
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LP . Substitute $+ using ([35]) we see that the significant terms in 1d + ($+ + w + <&+) 
are: 

l J D + ($i)+,lr> + w + ($i)+ 
Now, since both w + , ($i)+ has the same strictly-upper-triangular matrix struc- 
ture on D + (this is the way the second perturbation scheme was designed), the 
term l£),w; + (<fri)+ is effectively 0. Consequently, the main contribution in (|37[) 
comes from 

1 d+ {<S>i)+ + 1d + Aw+ = 1 D+ + 1 D+ Aw + ^) 

= l D+ ( - C+(1d+ Aw+) + A W +) - C_(lr, + Aw+) 

= -l D+ C-(l D+ Aw + ) 

which is small L' p . This gives us the desired estimate. When k > 2 the argument 
can be repeated and the endpoint case q = oo is allowed. 

(ii) In this part we assume k > 2. By part (i), we also have uniform boundedness 
of (1 — Cto^,) -1 in L 2 as t — s- oo. Using Lemma T4. 11 for w\ = wl and t«2 = w®, we 
can control limsup A ^ co |A(M*(A) - Ml(A))| by: 

\\w<s> -w L \\ 2 \H 2 (w 9 ) + H 2 (w L )j + \\w$ - w L \\ 00 H2(w i )H 2 {w L ) 

(38) + | J(w+-w+)w+\ + \ J{w^-wl)wl\ 
+ limsup | AC(w$ — Wi)(A)| 

A— >-oo 

First two terms : By Corollary 15.51 

H 2 {w L ) < t~WPn 

thus using (|36j) this estimate is also true for H 2 {w§). Consequently, using Lemma l5.ll 
the first two terms in (|38|) are controlled by 

recall that kg :— max{0, hi,..., fcjv}). 

Two middle terms : We'll estimate J 1d + (u>J —w^)w^, the other cases can be done 
similarly. By triangularity, 

l D+ ( W i - W t) W t = 1 D+((®2)+Wl + ($l)+(*2) + w£ + W + ($ 2 ) + wt) 

The contribution of 1d + (^2)+^£ can be decomposed into 

Since ($2)+ is small on D + , using Cauchy-Schwarz the second integral is controlled 
by t~^ k ~^H 2 < fe e +1 > . To estimate the first integral, we'll use analytic continuation 
to C + and Cauchy theorem. Eventually it suffices to show that along any vertical 
ray V in C + originated from arbitrary stationary point, J r (^> 2 ) + C+(w^) is small. 
Using the complex variants of Lemma 15.11 and Lemma 15.41 (i.e. Lemma 15.101 and 
Lemma T5. lip , this follows easily: in L 2 (T), ($2)+ have strong decay and C+(w^) 
has some nontrivial decay as t —¥ 00. Thus, by Cauchy-Schwarz we get an estimate 
of size 

t -( k +2(kg+l)) 
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The contribution of other terms can be estimated similarly. It might be conve- 
nient to notice that 1d + w + ($2)+w^ has at most one nonzero entry and the above 
argument should be applied to that entry. 
Last term : We'll show that 

limsup \\C{w+ - w+){\)\ < f-^+Ww) 

A^oo 

the estimate for w^, — u>2 is similar. Now, 

— w\ = ($1)+ + ($2)+ + other terms 
= C+(Aw + ) + other terms 
where the other terms are in L 1 . Notice that C+C_ = and C_C_ = CL, so 
limsup |AC(C+(Au; + ))(A)| = limsup | AC + (Au; + )(A)| 

A— >oc A— >oo 

= 0{t~ k ), by linear theory 

The other terms are (<J>i) + ($ 2 ) + , w + ($i) + ($ 2 )+, w+($ 2 )+, which are 

in L 1 . Thus, 

limsup |AC(other terms) (A) | = — / other terms 

A->co 27T Jjg 

By orthogonality, the contribution of ($i)+(<f>2)+ in J R is 0, and the contribution 
of w + ($i) + (<f> 2 )+ becomes 



f c_( w +)($ 1 ) + ($ 2 )+ < ff 2 MII($i)+($ 2 



using Cauchy-Schwarz. Since ($1)+ is small on Z?_ and ($2)+ is small on D + , the 
above inequality gives an estimate of size t 2 W+ T 5 i~( fe ~ 3). 

Notice that by triangularity, J 1_d + w + ($i) + = J 1.d_iu + (<I>2)+ = 0. Now, the 
argument used to estimate the two middle terms of (|38| can be applied to estimate 
J 1_d_ui + ($i)+ and J l_D + ui + (<i>2)+- Consequently, we have an estimate of size 
£ _ 2(* e +i) i~(k-i) f or contribution of w + ($ 1 ) + and u> + ($ 2 )+- 

(hi) Suppose that k > 1 and ||(1 - C Wl )- 1 \\ L p^ L p,\\(1 - C^)" 1 \\ L p^ L p < p 1 
(as t — s> 00) for any 2 < p < 00. For p,q £ H 1 ' with sufficient decay, we'll have 
$1, $2 & L°° . This allows us to use ©,((9]) and (fT0|) . but these decay requirements 
of p and q will not contribute to the finiteness of the implicit constants. 

Under the assumption k > 1, ||$i||oo,||3>2||oo are however not known to be uni- 
formly bounded as t —¥ 00. Consequently, we won't have uniform boundedness 
of (1 — C^) - 1 in L p . On the other hand, we know that $i,$2 are uniformly 
bounded in L°° _ as t —> 00. Using ([8|),((9]) and (fT0|) . it is then not hard to show 
that (1 — Cw^)" 1 is uniformly bounded from L p+ n L p to L p for any 2 < p < 00. 
This allows us to use a variant of Lemma B~T1 (see the remarks after this lemma for 
details). 

Now, the same argument as in (ii) can be repeated, where to improve the estimate 
of (ii) to 

O e (£~ (fe ~^ + 2 < fc 8+ 1))+£ ) 

it suffices to show that for 2 < q < 00 

||w* - w L \\^_ (H q {w L ) + H q+ {w L )) < r (fc -'-^ 
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(after that choosing q sufficiently large will give us the desired estimate). We note 
that in other places an e in the decay order might also be lost, one reason is the 
lack of an uniform bound as t — > oo of ||$i||oo and H^Hoo (we need to use L°°~ 
instead of L°°). Now, the above inequality is a consequence of (|36|) . □ 



Remarks: Let A, B, C respectively be the decay estimates for H2(wl), \\^d±(w ( / ) — 

WXf) H2 and \\w,p — wl\\oo- Then typically (modulo an e) our argument gives an 

overall estimate of AB in part (hi), and an overall estimate of max(AB, A 2 C) in 

1 

part (ii). We'll always have A = t 2 <- k e+ 1 ) , thanks to Corollary 15.51 and the correct 
phase-weight relation in our weights. 

8.2. Reduction of phase. In this section, we study a RHP localized to a small 
neighborhood of the stationary points of 9. Let Ao be a stationary point of 0. 
Below, we'll show the phase reduction of 9 near Ao- 

By translation symmetry, assume that Ao = and is of order fco. Let the weights 
of our RHP be denoted by w ± . On the neighborhood of where the weights are 
supported, approximate 9 by 

0(01 + ^J- x ko+1 

m+ (fc + 1)! 

and keep the value of 9 near other stationary points. Let 9$ denote this approxima- 
tion. Without loss of generality we can assume that the portion of supp(w) near 
is small enough so on a neighborhood of this part 9 is (fco + l)-time differentiable, 
with 

9' (0) = ■ • ■ = 0< fe °)(O) = 0, 9 {ko+1 \0) ± 
0(*D+i)(a.) = 6K feo+1 )(0) + 0(\xf), > 

Below we show how to reduce 9 — t 9$. The following assumption is required only 
if 9 ^ 9q in the above neighborhood (since otherwise we don't need this reduction). 

(*) p, q have two L 2 derivatives on a neighborhood of if fco = 1. Three L 2 
derivatives are required if fco > 3, or fco = 2 and p(0)q(0) < 0. 

The regularity assumption on p, q for fco > 3 may be improved given better 
understanding of model RHPs associated with a stationary point of such orders. 
More precisely, we don't know if the corresponding resolvent operator (I — C w ) 
associated with such model RHP is invertible on L p for large p, therefore in our 
direct scheme we need some nontrivial decay for HAuiH^, which is the source of 
this extra requirement. For this reason, if \pq\ are very small at a stationary point 
of such order then this extra requirement can be removed. 

Notice that is also a stationary point of order fc for 9 : 9 ( Q ko+1) (0) = 6>( feo+1 ) (0) ^ 
0. Furthermore, by L'Hospital's rule and continuity of 9^ ko+1 ^ at we have: 

i im |M = ... =lim e;;;^) =1 

«ofl'(i) x^o 6»( fc « +1 )(a;) 

Thus, our neighborhood of can be chosen such that > 1 for every x ^ 
in this neighborhood. In particular, this means 6(x, s) := s9q(x) + (1 — s)9(x) has a 
stationary phase point at of order fco for every s € [0, 1]; furthermore, for x near 
0, 4-6{x, s) has the same signs as 6'{x). 

Let Wq be the weight obtained from after replacing 9 by 9q. For simplicity 
of notation, let Aw — wq — w, which is supported near 0. To reduce from w to 
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wo, we will proceed exactly as in the reduction from to wl- For convenient, in 
the following propositions we'll denote by k > 2 the number of L 2 derivatives that 
p and q are required to have in the assumption (*). The argument of Section |8.1 
works once we show the following estimates and its complex variant: 

Proposition 8.3. For any sign combinations, 

(39) limsup|Ac(lz 5± A w ± )| < r w+1)l ^ +c 



x 

Furthermore, if ko > 1 then for 2 < p < oo 

(40) H p (l D+ Aw), H p (l D _Aw) < n J ax(i _Cfc_a+ »\r (fl+ » ) W +e ) 

If ko — 1 and 9^ is assumed L r integrable near then the above estimates remains 
true after the following adjustments: ko + 1 is replaced by ko + 1 + c(2,p), and j3 is 
replaced by j5 — c(2,p) for 2 < p < oo, or f3 — 2c(2,p) for p = oo. 

Note: c{k,p) is defined in ([23| . In particular, c(2, 2) = ^ and c(2, oo) = Recall 
the assumption (B) on 9, which in particular says that 



if fco = 1 
0, if fc > 1- 

The reason for this assumption will be explained shortly. 

Proposition 8.4. Let 1 < p < oo. If T is a ray originating from a stationary point 
that forms a nontrivial angle with M and ko > 1, then 

\\C(l D± Aw)\\ LP{r) < max(H fe - 1 ),t-^ ) W+ e ) 

If ko = 1 then this estimate has to be adjusted as in Proposition [ 



Before proving these propositions, we make some comments. Recall that the 
phase reduction and separation of contribution will be iterated through the list of 
stationary points. Using the above estimates, the effect on u(t),v(t) can always be 

controlled by AB + A 2 C, where A = t~ 2(k e +1) , B is the estimate for H 2 (Aw), and C 
is the estimate for H ao (Aw), and all we need is I? boundedness of (1 — C™)^ 1 , which 
is always the case (posteriori). However, if (1 — C w )~ x is also bounded in LP for large 
p, we can improve this estimate to AB (modulo an e in the decay exponent). For 
this reason, if there are more than one stationary point we'll always start out with 
stationary points of order > 3 and defocusing secondary stationary points before the 
rest. The main reason is whenever these bad stationary points are involved, we will 
not have L p boundedness, but after separating them out L p boundedness suddenly 
becomes available, which enable us to avoid strong assumptions on regularity of 
9,p, q near the good stationary points. 

Below we show that our estimate of the effect on u(t), v(t) can be controlled by 
something decay better than t k e+ 1 . If we only use L 2 boundedness of (1 — C^,)" 1 
then we'll have a decay of t k e +i+ d o)+ € _ 

do = r4r > 

ko + l 

assuming three L? derivatives of p, q near and boundedness of 9^ 4 K These con- 
ditions are satisfied when ko > 3 or ko = 2 with p(0)q(0) < by our assumptions; 
we note that the boundedness of 9^ can be changed to some high L r integrability 
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condition similar to the assumption on primary stationary points; but we avoid 
doing that here for the sake of simplicity. 

On the other hand, if we assume high LP boundedness of (1 — C^) -1 then the 

effect on u(t),v(t) can be controlled by t ( 2 < fc » +1 > +d(l)+ \ where d > 2{k l +l) ■ More 
specifically, assuming only two L? derivatives of p, q near we can get 

d = min(-, (ft + -) ) > ' 



^2' v " 2'fco + r 2(fco + l) 
if ko > 1. When ko = 1, if #( 3 ) is assumed to be LP integrable near then 



3r 



In this case, using the assumption ft > - we have 



,„ 1 2 , 1 , . 1 1,1. 1 
(/3+ 2-37 ) 2Tf )>( 2 + 37 ) 2Tf ) -I 

Consequently, only two L 2 derivatives is needed for p, q near if is a good sta- 
tionary point (i.e. primary or focusing secondary). On the other hand, in that case 
the reduction regarding unique solvability and resolvent bound needs to be looked 
at more carefully since we won't have decay for ^(Aw) as t — > oo. In this case, 
the above assumption on ft saves us. Indeed, by (|40|) we have 



HooQ-d^w) = o(i) + o(r c ) 

as t — > oo for some c > 0. More specifically, 

r= /i£r- £ < iffco>i; 

iffco = l- 

so c > 0. Now, the implicit constant of O(l) is (modulo a positive power) propor- 
tional to the size of the current neighborhood. Thus, it can be made arbitrarily 
small. For large t, a posteriori the resolvent norm ||(1 — C wo )~ 1 \\ on LP can be 
controlled by an absolute p— constant depending essentially on the values of \p(x) | 
and \q(x) \ at the stationary points. For p in any given compact subset of [2, oo), if i 
is large then unique solvability follows by choosing our neighborhood small enough 
and using Newman series. Furthermore, on those L p we can control ||(1 — C^) -1 !! 
by a comparable p-constant, which is good enough for our applications. 

Below we'll prove Proposition [831 the proof for Proposition 18 .41 is similar. We'll 
use a lemma whose proof is a simple application of Fubini's theorem: 

Lemma 8.1. Let f(x,s) be in L^LL].' 1 ((Q , 1) x K), i.e. f Q \\f(x,s)\\ H i,ids < oo. 
Then C± and ds commute: 

(41) (C± J /(., s)ds) (y) = J (c±/(., »)) (y)ds 

Proof of lemma. Since LL^' 1 C L^, we know / S L\L\. Fubini's theorem gives us: 

LHS(0 = 1{±?>o} ^ I f f(x, s)e~ ix t dsdx = [ l {±? > 0} /(£, s)ds 
Z7T JrJo Jo 
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Similarly (use the above argument for s) := /(£, s)l{±£>o} S L^L^ instead 
of f(x, s) and the inverse Fourier transform instead of the Fourier transform), 



F~ 



l U g(.,s)ds)(y)= j F- 1 (g(.,s))(y)ds= f (C±/(., s))(y)d, 



□ 



Proof of Provosition \8.3[ To show (|40p and (|39[) . we want to estimate the respective 
Hardy projection of Aw. Below we'll estimate 

\C + (l D+ Aw-)\ = \C + {l D+ ^P{e- m - e**°))\, 

the other cases are similar. For simplicity we'll assume that ko > 1, the case kg = 1 
is entirely similar and the adjustments in this case followed from the corresponding 
adjustments in Lemma |5~21 Below we'll show the argument is entirely similar 
for 439]). 

Write l D+ cj)oP(e- lte - e- tte °) as: 

it [ f(x)e- it0 ^ s h D Jx)n o (x)ds 



ii 



where recall that 0(x, s) = s6o(x) + (1 — s)9(x) and 

/(re) :-0 o (x)(0 o ( 2 ;)-0( a ;))p(x) 

Note that /(x) is compactly supported, continuously differentiable, furthermore / 
vanishes at the endpoint of D + . Thus it is not hard to see that /1_d + ^o S L^H^' 1 . 
Apply the above lemma, we have 

C + (cb Q P(e- uf> - e- ite °)l D+ )(X) = it j C + (l D+ f(., S )n e- M ^)(X)ds 



Under assumption (*) and by the A-B decomposition lemma, Qq has an + 
decomposition in a neighborhood of the stationary point 0. Now, / has multiplicity 
ka + 1 + (3 at up to the k th derivative. The last claim is a consequence of the 
Holder condition of Q( k ° +1 '> at 0. By Lemma \5A\ we have 

\\C + (l D+ f(.,s)n )(y)\\ p < s max(t-( fe - 1 +^,r (fco+1+ ^ + ? ) ^) 

Below we'll show that the dependence on s in the above inequalities are mild, so 
that after integrating over s G [0, 1], the implicit constants remain finite. 

Since / is compactly supported, the s-dependence of the implicit constants is 
originated from the phase dependence weights in the weighted sum fg t k- 

d 2 d k+1 1 

\-^9(x, S )r ■ ■ ■ \-r^e(x, S )r 



>dx 2 v ' ; dx 2 y ' n |^0(x,s)|( fc +«i+-+«*)-VP 

Since s £ [0, 1], we are not concerned about those s's appearing in the numerators. 

For those s's that appear in the denominator, we note that: our choice of 4>q 
ensures that on its support, \4-0{x, s)| ~ uniformly over s G [0, 1]. Therefore 

we can remove the s-dependence in our estimates and complete the proof of this 
proposition. □ 
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8.3. Separation of contributions. In this section, let Ao is a stationary point of 
order k of 9 and w is a pair of weights supported near the stationary points. Let 
wq denote the part of w supported near Ao and w\ — w — wq be the remaining. 
We'll assume that both Wo,Wi have correct-phase weight relation and the RHPs 
associated with wq and with w\ are uniquely solvable. Furthermore, we assume 
quantitative bounds — C w . )~ 1 \\lp^lp < 1 as i — >■ oo for some 2 < p < oo. 
Regarding unique solvability, we'll prove that 

Proposition 8.5. The normalized RHP with weight w = "Wq + ui\ is uniquely 
solvable for large t, indeed 

l 

(42) ||(i -c^)- 1 !!^^ < i + Y. IK 1 - C^U^lp 

3=0 

Proof. To show that / — C w is invertible for large i, we'll use the parametrix / + 
E}=o Cw, {I ~ C Wj )" x = / + E^oi 1 - G Wi )- x C W3 , introduced in |3T]: 
l 

(/ - c w )(i +5^c^(j •- c Wj y i ) = i- J2 c Wi c Wj {i - c^y 1 

3=0 ' 0<i#j<l 

1 

3=0 0<ijij<l 

By Corollary 15. 71 for i ^ j, \\C Wi C Wj \\lp^lp decays as t —> oo: 
\\C Wi C Wj \\ LP ^ LP <r*t*TPT 

Thus, for large t we know that X] < t ^<i C Wi C Wj {I-C Wj )~ l and X)o<j^<i( /_ 
C Wj )~ 1 C Wj C Wi have small norm. Consequently, / — C w is invertible, and 



\(I — C w ) 1 \\lp^Lp ^ \\I + ^ Cwj (-f — C Wj ) 1 \\lp^Lp 

3=0 
1 

< C^Wlp^Lp 



3=0 



□ 



For < j < 1, denote by Mj,/j,j the respective ingredients of the normalized L 2 
RHP with weights (wj ,w^). To separate the contributions of different stationary 
points, we'll need the following a priori estimates: 

Lemma 8.2. Assuming that near Ao the phase 9 is of the form a+b(x — Xo) k+1 . As- 
sume that locally p, q have two L 2 derivatives. Then if P C K so that distance(P, supp(wo) 
1 cmrf distance(P, Ao) ^ 1 i/ien for a.e. X E P we have: 







up(t) ' 



(43) Mo(A,*)=J + ( „; (f) Y +° f ( rw+f ) 

, A — Ao 



.ffere uq and vq are the contribution of Xq to the leading asymptotics of u(t),v(t), 
defined as in Theorem \l.l\ 
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Wc note that the above assumption on 9 is automatic because at this moment 9 
has already been reduced to a nice analytic phase near Ao- Furthermore, if k = 1 
then the Lemma remains true for p, q with one L 2 derivative. On the other hand, 
for purely harmonic analysis interests, the above asymptotics remains true for more 
general class of 9, for instance those having three Lipschitz derivatives near Ao and 
satisfying 

(the error term has weaker decay if f3 < 1), however the main issue here is that the 
Lemma requires higher regularity assumptions on p, q to allow for this generality. 
Indeed, the phase reduction of this Lemma is fairly expensive, it may require p, q 
to have three L 2 derivatives near Ao- This is the reason why for every stationary 
point we always carry out phase reduction before separating its contribution. 

It would be interesting to see if the decay estimate O e {t ^k+iy-r ) can be im- 
proved to O e (t~ T + I+e ) assuming enough regularity of p, q. That improvement would 
allow us to obtain a stronger decay estimate for the error term in the leading asymp- 
totics of u(t) and v(t), when there are more than one stationary point. 

The proof of Lemma l8.2l is a sequence of reductions; at the end it will be proved 
for the corresponding model RHP associated with Xj. 

For technical reason, from now on we'll assume that 



(44) 



limsup^^ ( \X(M ) 12 (X,t)\ + |A(M ) 21 
limsup^^ (|A(M 1 ) 12 (A,t)| + |A(Mi)ai(A,t) 




Since the reduction of u(t),v(t) implies the reduction of these estimates, this as- 
sumption is harmless. 

Assuming these a priori estimates, we will show below that the leading asymp- 
totics of u(t) and v(t) are essentially the sum of the leading asymptotics of uo(t) + 
U\{t) and vo(t) + v\(t). As usual, we'll suppress t for brevity. 



Proposition 8.6. Assuming the a priori estimates we have 

(45) limsup |a(m 12 (A) - (M )i 2 (A) - (Mi)i 2 (A)) | < e t~^4+v +e 

(46) limsup |A(M 21 (A) - (M ) 2 i(A) - (M^^A)) | < e t~^^ +e 



Proof. We'll prove (|45l l46l) using a similar argument as in [31]. It is not hard to 
see that m(A) := Mi (A) Mo (A) solves the following normalized RHP: 

m + (A) = m_ (A) J m (X, t) 

J m (X,t) = ((7-«;r)Mo,_)~ 1 ((/ + OMo, + ) 
To show the proposition, we'll show that as A — > oo nontangentially 

flimsup A ^. DO |A(M(A)-m(A))| < e r^f+^ +c 
[limsup^^ | A Off(m(A) - Mi (A) - M (A)) | = 

Here, for any 2x2 matrix A, we denote by D(A) the diagonal part of A and Off (A) 
the off-diagonal part of A. Both D(A) and Off(A) are 2x2 matrices. 
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The first task is to find a suitable factorization for J m that will help us prove 
the first estimate. Recall that Mo,± = ^o(I ± ^o )> here 

MO = I + Cu, /Z.Q 

If matrix multiplication was commutative then we would have 

(7 + wt)M , + = (i (I + Wo)(I + wf) = f i (I + w+ + w+) 
thanks to disjoint support, and similarly (/— w^)A/o i _ = /j.q(I— — Wq). However, 
we can only estimate ||//o — I\\2 by t~ 2 ( k + r > which is far from the right-hand side 
of (|45l |46| if k ~ kg . Consequently, to effectively compare m with M using weight 
comparison (i.e. the first scheme), we don't want to use the immediate choice 
I±wt = (I±wf)M 0<± . 

Instead, we'll use the factorization J m = (I — u' rn ) _1 (/ + ui+), where: 

I + w+ n := ^(I + wf )M ,+, I - u;~ := Mo V - w^)M _ 
The idea is to somehow cancel out the effect of /xo Now, 

7 + = A*o 1 ( / + wt)fx (I + w+) 

= (I + w+)(I + w+) + ^[w+^o^I + w+) 

= (I + w + )+Aw + 

Aw + := ^ 1 [w;+,^ ](/ + w+) = ^[w^no] 

because wq — on the support of wi. 

Similarly, I + w+ n = (I — w~) — Aw" where Aw~ := /i 7 " 1 ^ 7 ", Mo]- We will see 

1 

that with this choice of weights for J m , we can achieve ||Aw|| p < t fc « +1 for any 
1 < p < 00. 

From the proof of (fT4| . limsup A ^ 00 |A(M(A) — m(A))| can be controlled by 
\\Aw\\ 2 H 2 (w) + \\n - pL m \\ 2 \\w m \\ % + limsup |AC(Ak;)(A)| 

A— foo 

EE X + Y+Z 

Here, ||m — Mm II 2 can be controlled by 

(47) HCWII2 + ||Au>|| 00 .H 2 (u;) 

using ||(1 — Cm;) _1 ||z, 2 -kl 2 < 1 (by the previous proposition). Now by Corollarv l5.51 

H 2 {w) < i~ 2 < fe 8+ 1 ) 
Thus it remains to estimate Aw. For simplicity of notation, let 

/ n Ug(t) \ 

\x-Ao / 

For A e supp(uii), we have distance(X, supp(wo) U {Ao}) > 1, and (|43|) gives 

(48) M o(A) =/ + M °(A) + O e (r5<AT7 +e ) 

Since det Mq(A) = 1 on R, we have det(/zo) = 1. Thus /i 7 / 1 is essentially a rearrange- 
ment of entries of /io with additional minus signs in appropriate places. Thus, (|48[) 
implies: 

(49) Mo '(A) = I - MoW + O e (r ™ +£ ) 
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Now, notice that Aw is supported on a small set, so using estimates (|48| I49|) we 
have: 

Aw ± = (j-/i°)[™± ) / + / i°] + o e (rOT+ £ ) 
= k±, ^] + o e (r ^W +£ ) 

Since u>i and Aw are supported on small sets, the 0(t~ 2 < fc + 1 ) +e ) term remains small 
in L p for any 1 < p < oo. In particular, ||Aw||2, HAwHoo < < _ ' E + T , and thus 

X < t _ 2(fc 8 3 +l) 



Y < ||Ca™/|| 2 + 0(t 2(fc « +1) ) 

Below we'll estimate ||Cato/||2- Notice that [^i ,/4)] are algebraic sums of terms 
like C '^^ 1 A [ ) :E ' t " > or Wl wnere C(t) denote matrix valued functions of t (inde- 

pendent of a;) and have order t^ 7 ^ 1 . In addition, the weights u^'s have the correct 
phase-weight relation. Since \x — Xq\ > 1 for x G supp(wj t ), Corollary 15.51 now 
gives: 

||CW||a ^ \C(t)\t~WV + ft(fTO +( ) = O e (t"^ii +£ ) 
Thus, Y < £ t _2 ( fc e+ 1 ' +e . Now to estimate Z notice that Aw e L 1 , thus 

„ ,. f°° Aw(x,t) , , i , r 00 A , , , , 

Z = limsup A / — dx\ = — / Aiu(a;, i)air 

A^oo J-oo ^ — A Z7T J_ 00 

so using linear theory of oscillatory integrals (see Lemma I5.8[) the contribution of 

i i i i 

[tuf , Hq] can be controlled by t k+1 k » +1 , here ^fpj- comes from constants like C(t) 

and fc * 1 is contributed by the respective oscillatory integral. Thus, Z and hence 

3 I c 

limsup A _^ oc |A(M(A) — m(A))| can be controlled by t 2 < fc e+ 1 > 

To finish the proof of our proposition, we need to show the second estimate 

(50) limsup | A Off(m(A) - Mi (A) - Mo (A)) | = 

A— >oo 

Observe that for any two 2x2 matrices A, B, 

|Off(A0)| = |Off(A)D(B) +D(A)OS(B)\ 

< |Off(A)||D(B)| + |D(A)||Off(B)| 

([50)) now follows by applying this inequality for A = M\ (A) — I and B = Mo (A) — I, 
noticing 

lim |D(M (A) - 7)| = lim |D(Mi(A) - J)| = 

A— >oo A— >oo 

limsup |AOff(M (A) - I)\ = limsup |AOff(M (A))| = O t (l) 

A— foo X^-oo 

limsup |AOff(Mi(A) - I)\ = limsup |AOff(Mi(A))| = O t (l) 

A— foo A— >oo 

□ 
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8.3.1. Reduction of the a priori estimate ^Bjj in Lemma \8.Sl For future reference, 
we'll set up a framework to reduce (j43|) along our reduction of u(t), v(t). 

Let \x and /j,q be associated with two pairs of differentiable weights w and wq 
that are localized to a given stationary point, which we'll assume of order kg. 
Assume that for each x the triangularity of w{x) and wo{x) agree, which in turn is 
invariant on each D + and £>_. For convenient, let A/j, = /i — (1q and Aw = w — u>o. 

Our setting will be as usual: the restriction on D + , D- of both Aw + and Aw~ 
are oscillatory functions of the form /e ite , where /(0) = and / has one L 2 
derivative. Since we have already done the phase reduction, in our setting O is of 
the form a + bx ka+1 , but the argument works for more general phases. 

Recall that P is away from {0} U supp(w). In this localized case, we'll also have 

distance (P, supp(u>o)) > 1 

In our final reduction (to model RHPs), the model weights will not be locally sup- 
ported. To overcome this, we'll exploit certain analytic continuation of the model 
weights to deform the Riemann-Hilbert contour M. Effectively, this deformation 
"moves" P away from the supports of the (deformed) model weights (so that the 
current argument can be reused). 

Proposition 8.7. If w and wq have the correct phase-weight relation and are sup- 
ported away from P then 

I|A/x||l-(p) < \\A(J,\\ 2 + t~*o£+VH 2 (Aw) + i~Wrr 
In applications of this proposition, observe that thanks to the correct phase- 

1 I 3 I 

weight relation in w and u>o, H2(Aw) can be estimated by t k »+ 1 or t 2 < fc o +1 > 
depending on availability of regularity. Also, as a consequence of the reduction 
of u(t), v(t), we'll have an a posteriori estimate for ||A/i|| 2 . To see that, note 
that A/j, = in the second scheme and small for large t in the first scheme, using 
Lemma [4.11 Often we'll have ||A/x||2 = 0(t~ 2 < fc o+ 1 ) +e ), which is enough for the 
desired reduction. 

Proof. Writing Afi = Caw^o + C w Afi and using Cauchy-Schwarz's inequality in 
the second term, we have: 

(51) I|A/x|U~ ( p) < \\C Awf i \\ L ~ {P) + HAHI2IMI2 

Thus, it comes down to estimating ||C , AiuAto|| i oo(p- ) - Intuitively, we would expect 
this term to be small if [Iq—I was sufficiently small (in some L p ), since by Lemma [5.8l 
it is not hard to see that Ca w I is small in L°°(P). However, from Lemma I5T41 (more 
precisely, Corollary 15.51) it is only known that 

llMo - ilia < i"W+u) 

which is not enough. To get the desired estimate, the main idea is to exploit the 
fact that P is away from {0} U supp(w) U supp(wo). 
Now, on P, write Caw/J-o as 

C Aw I + C Aw (c^/io) = 0(r^) + C+ (c u , mo Aur) + C_ (C wof j, Aw+) 

= 0(i~W)+X + Y 
Below we'll prove that X is small on P; the proof for Y is similar. 
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Write X = C+(E)+C+(F), where E = C+(^Wq)Avj- and F = C-(fi w£)Aw- . 
Notice that E and F are not symmetric in general because C+ is being applied to 
them. However, if C+(X) is evaluated on P then this symmetry is available as both 
E and F are supported away from P. 

Keeping in mind that ||C+(Aw; - ) H2 and ||C+(Aui~)|| £ =o(p) are small, we decom- 
pose F to establish the appearance of these terms in C+(F): 

F = C^ w+)C + (Aw-) - C^ w+)C^(Aw-) 

= C+(fi w+)C+{Aw-) - fi w+C+{Aw-) - C-(ij,ow£)C-(Aw-) 

Under C+, the last term vanishes, while the first term is unaffected. Using the 
distance condition, we can estimate ||C+(i 7 ')||£oo(p) by 

C + (fi Q w+)C+(Aw-) - C+(( M o - I)w+C + (Aw-)) - C + (w+C+(A w -)) 

= o(\\C+(Aw-)\\ L *° (P) + H 2 (w )H 2 (Aw)) - C+ (w+C+(A w -)) 

For A G P, let g(x) = Then 

C+(w+C + (Aw-))(\ ) = fgC + (Aw-) 

c-ia) c+(Aw-) 



< \\C^g)\\ 2 H 2 (Aw) 
Consequently, using the correct phase- weight relation in both w and wo we have 
\\C+(F)\\ L ^ iP) + rw^H 2 {Aw) 

□ 

9. Reduction to model cases (II): Deift-Zhou's steepest descent 

ARGUMENT 

In this section, we consider a RHP that is localized to a small neighborhood of a 
stationary point and has a nice analytic phase. The goal of this section is to reduce 
this RHP to a model RHP which will be explicitly studied in Section [TUl Without 
loss of generality we can assume that this stationary point is and is of order k, 
and the phase of this RHP is Q(x) = a + bx k+1 for a, b € R, 6^0. Let w^ 1 denote 
the weights of this RHP. 

Let 6q(A) be the local approximation of 5(A) near 0, and let £)_ := {9' < 0} 
and D + := {8' > 0}. Recall that Sq satisfies the following scalar RHP: 



5 a+ (x) = 5 -(x)(l D+ (x) + [l+p q ]l D _(x)y i 



x G 



here po = p(0) and go = <?(0). Also, D± can be made explicit by considering parity 
of k and sign of b; each of them is either 0, K_, K + , or R \ {0}. 
The current pair of weights w = (w~ , w + ) is of the form 



(w 
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here </> is a normalized cutoff supported near 0. Note that we implicitly used the 
jump relation of S to write w in the above form, w will be reduced to: 



(w M ,w+) = 



sr 2 




if x e J D^ 
if x e D- 



0- l+po<?o 

Intuitively, this means that as t —> oo the RHP (w~,w + ) localizes at 0: <fr(x) — > 1, 
p{x) -¥ po, q(x) -t q 0l and 8(x) -> 6 {x). 

The main difficulty in the above reduction is the lack of L 2 integrability of the 
model weights wm, despite the fact that they are still in L°°. This prevents Beals- 
Coifman's operator formulation which is essential to our perturbation schemes. To 
get around this issue, the idea introduced in [9] is to exploit analytic continuation 
of wm and deform the model RHP to a suitably chosen contour T, on which % 
has strong decay. We note that this approach was not taken by [31] . 

We'll use the following adaptation of a notation in [TU]: For any continuous /, 

let 

WM = 

1 + ix" 

here N is a large natural number that depends only on fc, so that [/] decays fast 
enough on K. Large choice of N in particular ensures the applicability to [/] of 
Lemma |5. II and related propositions. Intuitively, [/] is an analytic approximation 
at of / that stays in LP spaces. 

For simplicity of notation, let M.+ denote [0, oo) and M_ denote (— oo, 0]. Consider 
ro,ri, T2, T3, T±, ^5 six rays originating at such that Tq = M + , L3 = R_ and the 
others form small angles with R. More specifically, we will take 

ri = e ia R+, T 5 = e- ia R +1 

L 2 = e"' ia M_, T 4 = e ia R_ 

for a sufficiently small a > 0, say a = ™. Small value of a ensures that [/](sc) is 
analytic and decays strongly in the two angles formed by T^, T\ and T2, T4. 

For any (reasonable) oriented contour E, we can define (nontangential) ± bound- 
ary values on E of a function analytic on C \ E using the following standard con- 
vention: Along E following the local direction, the left side is + and the right side 
is — . Consequently, for any / G L P (E), 1 < p < 00, we can define the nontangential 
boundary values C^f of the Cauchy operator Csf. 

Let r = {J i= Q^i. On r, we'll orient Fo,r2,r4 outwards (i.e. — > 00), and 
ri,r3,r5 inwards (i.e. — > 0). In particular, i C T is oriented naturally. This 
orientation makes T a complete contour in the sense that it divides C into two 
(disconnected) regions: one stays entirely on the — side of T, and the other stays 
entirely on the + side of T (it is equivalent to saying that T is the common bound- 
ary of these regions, clockwise for one and counter-clockwise for the other) [T2] , 
For convenience, we'll use —Tj to refer to the reverse orientation on Tj, and |r\j| 
whenever orientation is not taken in account. 

The main advantage of T being complete is: the Cauchy operators Cp satisfy 
standard properties of (see for instance, Zhou [31]): 
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For any i ^ j, let Tij C M 2 denote the angle formed by rotating \Ti\ counter- 
clockwise to | r j | - With this notation, Toi U T23 U is the + side of T, and 
r i2 U r 34 U r 50 is the - side of T. 

The following observation is crucial in the analysis of this section: w M and w M 
have analytic continuation to small angles (pivoted at 0) respectively in the upper 
and lower half planes, where they decay strongly. More precisely, w M can be nicely 
continued to {z 6 C : Argz e (0, ppj) U (it — ^xx,7r)} and can be nicely 
continued to {z e C : Argz € (— 0) U (— tt, —it + jrpt)}- To verify these claims, 
it may be convenience to consider different cases based on the parity of k and the 
sign of Q( k+1 \0). 

Intuitively, the above strong decay of the deformed model weights is not a co- 
incidence, it is rather expected. If instead of real- variable methods, Lemma 15.41 
and Lemma 15.11 were proved using a steepest descent argument (assuming relevant 
things were analytic) then this type of decay would be exactly what we need. Given 
a complete contour E, in order to have C^T (w~) to be small it would be convenient 
that w~ has a decaying analytic continuation to the — side of E. Similarly, for 
(w + ) to be small we want w + to have a decaying analytic continuation to the 
+ side of the contour. Whenever a pair of weights has these properties, we'll say 
that they are natural. The naturality of our weights are indeed guaranteed by the 
correct phase- weight behavior in (w~,w + ), which is a consequence of the studies 
in Section [5] 

The reduction in this section consists of two major steps: 

Step 1: Preparation for steepest descent. In this step, we'll reduce w — > wm, &n 
analytic approximation of wm that stays in L? (R) : 



£|p]--" H 
,0 



5~ 2 < 




if x e D 0+ ; 
if x 6 D - 



'-l+pqi 

Step 2: Steepest descent. In this step, we'll reduce wm to the F-deformation of 
the model RHP associated with wm- To do this we'll deform the former to T. In 
this reduction, we largely follow the steepest descent method of Deift and Zhou 
[3 HO] ? with several adaptations to our settings. 

9.1. Preparation for steepest descent. Recall that the set P used in (|43|) is of 

distance > 1 from 0. Without loss of generality, we can assume that 

distance(P, supp(0)) > 1 
The reduction is broken up into two sub-steps 

w — > c/)Wm — > wm 

In the second sub-step, the localization argument as in Section 18.11 can be used 
to control the effect on u(t), v(t) by terms decaying highly as t — >• oo; hence the 
only task is to prove the reduction of (|4"5)l : observe that the same argument as in 
Section IB .3. II can not be repeated because wm is not localized. 
Sub-step 1: (w — I 4>wm)- Let Aw = w — cjywM- Using the two perturbation schemes 
described in Section |4] and implemented in Section [8j the following estimates are 
enough for our reduction: 
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Proposition 9.1. (i) If p and q has two L 2 derivatives near then for any sign 
combinations 

(52) H p (l D± Aw) < p r (1+ p } T^r\nt, 2 < p < oo 

(53) H 00 {l D± Aw) < e t-^r +e Ve>0 

(54) Kmsup|Acflr) ± At« ±N )| < 

(ii) If p and q has only one L 2 derivative near then \5J$ remains true, and 

(55) H p (l D± Aw) < p maxfr^r^+p^), 2 < p < oo 

(56) H 00 (l D± Aw) < 0(l)+0(fWU +£ ) 

in \56\) the 0(1) term is (modulo some positive power) proportional to the size of 
supp((p) . 

Proof. The above estimates are indeed direct consequences of Lemma 15.41 (in par- 
ticular Corollary I5.6j) and Lemma 15.81 For instance, we'll consider the setting in 
(i). Each lrj ± Aw ± has at most one nonzero entry \r> ± £lfe lW , with: 

(i) correct monotonicity for 9 on the support of this entry; 

(ii) / vanishes at with multiplicity 1 up to the second derivative, and 
(hi) fi has an A2 + B2 decomposition. 

For example, on D + the nontrivial entry of Aw^ is 

(^5- + 2 q - <M;W e = n„V(9 - 

The argument is similar for the setting in (ii), note that the corresponding multi- 
plicity will be i. □ 

Using the above proposition, the effect on u(t) , v{t) can be controlled by O e (t~ fc+r +£ ) 
if p and q have two L 2 derivatives near 0, or O e (t~ 2 < fc + 1 ) +e ) if they have only one 
L 2 derivative. In the latter case, choosing supp(4>) small enough will ensure that 
the reduction doesn't harm the unique solvability and the corresponding resolvent 
bound of the RHP associated with w (see Section l8T2|) : alternatively we can perform 
a suitable approximation as in Corollary 17.21 

Regarding the a priori estimate (|4"3")l . we'll use the argument of Section [8.3.11 
Note that this a priori estimate is not needed if N = 1 (thus one L 2 derivative for p, q 
near is enough). When TV > 1, for the reduction of (1431) we want || A^Hs = t~ 2 ( k + 1 ) , 
and this is satisfied if p and q have two L 2 derivatives near 0. 
Sub-step 2: (4>u>m to wm)- The only thing we have to show is the reduction of (|43|. 
As usual, P denotes a subset of K that is away from and supp(0), and Ao is a 
generic point in P. The main idea is to exploit analytic continuation, which was 
not available for us in Section 18.3. 1[ to move the weights away from R. Effectively, 
this re-establishes the distance condition between P and the support of the weights, 
which was essential in our reduction of ()43|) in the localized case. 

Proposition 9.2. If |^ 3| ) is true for (i (of wm) then it is true for fi (of 4>wm)- 

Proof. We'll use a variant of a steepest descent argument in |31j to show this propo- 
sition. 
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For simplicity of notation, in this proof the pair of weights wm will be simply 
denoted as w. Let /x be associated with w, and the solution of this RHP is denoted 
by M. We'll use A/u, Aw ± to denote the respective differences of two RHPs. 

Since distance(P, supp(0w)) > 1, Cauchy-Schwarz's inequality gives 

I|A/i||loo (p) < \\C Aw n\\ L oo {P) + ||A/i|| 2 ||<H|2 = \\ c AwV\\l°o(p) + 0(t~*&V) 

To estimate Ca-wH on P, we'll write it as C + (/iAw~) + C_(ttAw + ) = X + Y. Below 
we'll prove that ||X||^oo(p) is small; ||Y|| £ oo(p) can be estimated similarly. 
Using /j, = M + (I + u>+) _1 = M_(7 — ui - ) -1 , an explicit computation gives 

(57) l D+ fx u =1 D+ M U _ =l D+ ^M 11+ -Sf[q(l+pq]e ue M 12+ ^ 

(58) l D _n u =1 D _M U + =lD-(M u _+5l 2 [q/(l+ P q}e ue M 12 ^ 

Thus, has analytic continuation to a good sector (with pivot at 0) on which 

it is normalized to 1 as z — > oo nontangcntially. Similar observations can be made 
for //h1d_. 

Without loss of generality, we can assume that e( fe+1 ) (0) > 0. The picture (and 
our argument) is reflected across the imaginary axis if 6^ fc+1 ^(0) < 0. 

Under this assumption, the continuation sector can be made precise as follows: 

(i) k is odd: D + = M.+ ,D- = M_, and a good continuation sector for jUnlr> + 
is {—it < Argz < 2^}, and a good continuation sector for /^h1_d_ is {0 < Argz < 

(ii) k is even: D + = R\{0}, D_ = 0, and a good continuation sector for ^h1_d + 
is {-7r - ^ < Argz < 277}. 

Repeating this computation, we'll deduce similar conclusions for every entry of 
/1 (one small difference: the diagonal entries converge to f in their continuation 
sectors, while the off-diagonal entries vanish at 00). 

The computation also reveals that (the restrictions to D± of) different entries 
of ji do not always have the same continuation sectors. Since our estimation of X 
involves deformation of R by nontrivial angles, the knowledge of the continuation 
sector is important. This is the reason why we'll need to estimate X entry-wise, 
unlike our previous argument for the localized case. 

An explicit computation, with triangularity taken into account, gives 

Aw - = Avails- MiiAwr 2 lD + A 
\/X22Aio^ilD_ A 1 2iAwjT,lr>+/ 

Repeating the above computation (for /Jn) on other entries of /x, it is not hard 
to see that: in every entry of /iAw - , the respective restriction of the relevant entry 
of fi has a good continuation sector that contains the lower half plane. 

Intuitively, when estimating C+(/iAu;~), this observation allows us to reduce the 
respective Aw~j1d ± to their Hardy projections onto C + (i.e. C + (Au)r.l Di )) which 
are known to be small. Combining this observation with a contour deformation 
argument, we will be able to show the desired estimate for X. 

Below we'll estimate X12; detailed computation for other entries can be done 
similarly. Let / be the analytic continuation of /xh1d + mentioned above. Then / 
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Let g = (/ — 1)C+(Au> 12 1d + ). Consider two cases: 
Case 1: k is odd. Then D- = «_ and D+ = R+. 

Define 1^ = Ti, T° = T 4 , and T = T° + U r°_, orientated by increasing order of 
the real part. Notice now g has analytic continuation from D + to {0 < Argz < ^} 
and from £>_ to {— 7r < Argz < — ir + ^} (observe that C + (Aw^ 2 1d + ) is analytic 
outside D + ). We'll show that: 
Claim: For A in the upper half plane, 



Proof of claim. We'll use contour integration (twice). For R > large, consider 
a triangular contour 7 consists of three edges: 71 := [0, R], 72 := e la [0,R], and 
73 straight line segment connecting R and Re la . This contour is oriented counter 
clockwise, and 



We'll show that the contribution of 73 vanishes as R — >• 00. Indeed, a simple 
application of Carleson's measure theorem (see for instance [17]) shows that: for 
any < a < ir and h e L 2 {M), and for T R := {R + e ia x : x G [0, 00)} (which 
contains 73): 



Consequently, by (|57|) . by the boundedness of 6, and by the fact that the respective 
oscillating phase of w + is bounded (indeed it decays) on 73, we have, 





l|C+(fr)IU»(r*) < l|C + (fc)lk 2 (c +) ~ \\C+(h)\\ LHm 





Similar argument shows that, for A in the upper half plane, 




Adding these two equalities, we get the desired claim. 



□ 
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Taking limit of (|59[) from the upper half plane to any Ao G P, we have 

X 12 (A ) = 0(r**r) + 5 (Ao)l R+ + t^t / -^-<fe 

Now, |/i(Ao) — 1 1 < t^^+i since /i satisfies the a priori estimate (|43|) . Thus, using 
the definition of 3, we have 

5(A )1r+ < t-^\\C + (Aw w l D+ )\\ L ^ {P) 
< t^fctr 

On the other hand, by Cauchy-Schwarz's inequality, 
1 f 9{z) <r 1 



2-ni 



I -^dz < ^\\f-l\\ La(r) \\C + (Aw^lD + )h*(ro) 

J T o Z — Aq Aq 



/r° 

< \\C + (Aw^l D+ )\\ LHT o 



Applying Lemma l5.11| we have 

\\C + {Awi 2 l D+ )\\ LHr0) <t~*&V 
In fact, write Aw^l^ = h(x)5?_e~ lte( - x )l D+ with 

h = {i-4>)\p] 

Note that 5?_ has the A\ + B\ decomposition, 1 — <f> is smooth and supported away 
from 0, thus h vanishes at with multiplicity 1 up to the first derivative. Now 
Lemma 15.111 gives us the above estimate. 
Case 2: k is even. Then L>_ = and D + = K \ {0}. 

In this case, define r+ := Ti, T°_ := r 2 . The rest of the argument is essentially 
the same. Similar contour integration argument gives 

i r ( ) i r ( ) f 5(A) ' < ArgA < a; 

— / ^\dz = — / ^\dz + \ 5 (A), tt - a < ArgA < vr; 

ZTTl J & Z — A 27TI J r o Z — A 

[0, a < ArgA < tt — a. 

From there, for Ao £ P we have 

X 12 (A ) - O (i-^r)+ 5 (A ) + -L / -i^-dz 

=> X 12 (A ) = 0(t-^) + 0(||C + (A Wl2 l D+ )|U 3 (ro)) 

By Lemma TS-lli \\C+(Awi 2 ^D + )\\L 2 (r ) is small. Eventually Xi 2 (Ao) is small. □ 

9.2. Steepest descent reduction to the r-deformed model case. The reduc- 
tion consists of two steps: 

Deformation of the pre-model RHP to V 

I 

Reduction to the r-deformed model RHP 

However, regarding unique solvability and resolvent bounds, we can reduce a step 
further to the model RHP wm on R (the undeformed model RHP). The lack of 
L 2 (M) integrability of wm is the main reason that prevents the reduction to this 
actual model RHP of the leading asymptotics of u, v and the a priori estimate 
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On the other hand, we can still talk about invertibility of 1 — C WM because w M 
remains in L°°(R). 

The following analogue is very useful to keep in mind: 

contour deformation ~ second perturbation scheme, 

reduction on T ~ first perturbation scheme. 

We'll see many similarities which support this intuition, one of which was already 
mentioned above: in the first perturbation scheme and in the reduction on T, we 
estimate the effect on u(t),v(t) by directly estimating the weight differences. 
I. Contour deformation 

In this section, for simplicity the weights associated with the pre-model RHP 
will be denoted as (w~,w + ): 



(w , w + ) = < 



SZ 2 




if x e D + : 

if x e D- 



n+pql 

Here, S is the solution to the scalar RHP defined in Section [SJ 0(x) = a + bx k+1 is 
a real valued phase, and D + = {9 > 0}, D- = {9 < 0}. 

Our goal in this section is to deform this RHP to the following RHP: 

m + (z) = m,— (z)ji<(z), z 6 T 
m + -IeH 2 (T + ), m_-/etf 2 (r_) 

The jump matrix jy is defined on T and equals to 

(i) I on r Ur 3 ; and 

(ii) (I + on r*i U T 2 ; and 

(iii) / — w~ on r4 U T 5 . 

Here on T should be understood as the analytic continuation of the respective 
restriction of these weights. For instance, on Ti, w + is the continuation of w + |r + . 

Indeed, if n solves the pre-model RHP (w~,w + ) then m(z) := n(z)(I + < l>(z)) _1 
solves the above RHP, where: 

fo, ifzer 12 ur 45 ; 
$(z) = lw+, ifzer iur 23 ; 

if z e r 50 u r 34 . 

(as usual, t is suppressed for simplicity). To see that m satisfies the L 2 normalization 
condition, it might be convenient to exploit strong decay of the above $ (which goes 
back to strong decay of the respective oscillating terms e ±It0 ). However, as we'll 
see, uniform boundedness of $ is good enough for this purpose. 

Since T is complete, the same theory on R is applicable: For al 2 fl L°° factor- 
ization jr — (I — + Wp), as long as (1 — C^) -1 exists on L 2 (T), we can 
find a function fir on T such that fj,r = I + C Wr fir , and the above RHP has the 
following unique L 2 normalized solution: 

m±(A) := 1+ C^(fi r (w^ + ))(A), A g T 

The following relations also hold Vz 6 V: 

m + (z) = fir{I + Wp ), m_(z) = fJ,r(I — w r ) 
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Since w~, w + € L 2 n L , we can recover u, w from the limit as z — > oo of zn(z) 
along any non-tangential direction in C \ R. Our choice of factorization will ensure 
that Wp G L°° fl L 2 n L 1 , so we can recover ur(t),ur(i) by taking limit of zm(z) 
as z — > oo along any non-tangential direction in C \ T. Since m(z) and n{z) agree 
inside Ti2 U F45, we get: 

Corollary 9.1. Assuming unique solvability (of both RHPs), the deformation doesn't 
change the recovered potentials u,v of the RHP associated with (w~ ,w + ). 

Remarks: 1. The above deformation can be carried out for the model RHP ()52[) . 

2. Our second scheme is based on a conjugation of the jump matrix, which 
corresponds to a right multiplication on each Riemann-Hilbert factor. This is alge- 
braically similar to the above contour deformation; the main difference is that the 
deformed and original RHPs arc on different contours. Thus, contour deformation 
can be seen as a contour extension followed by the second scheme. 

To enable contour extension, we'll show that the invcrtibility of the respective 
1 — C w are not affected when extending the weights to trivially to a larger contour; 
furthermore the old and new resolvent norms are comparable. 

9.2.1. Contour extension. Our setting is general: £ = £i U £2 is a union of two 
oriented contours (intersecting at finitely many points); w — (w + ,w~) is a pair of 
L°° weights defined on £ such that u;|e 2 = 0; Cs, w and C^ lyW are Beals-Coifman 
operators on £ and £1 respectively. Note that we do not require the weights to be 
in L 2 ; this allows us to apply the following proposition to both model and pre- model 
weights. 

Proposition 9.3 (Contour extension). Let 1 < p < 00. The invertibility of 1 — 
Ct,x,w and 1 — Cs jt0 respectively on L p (£i) and L p (£) are equivalent, furthermore 
the resolvent norms are "comparable" in the sense that the following s are equivalent: 

1(1 _ Ce,«>)~ ||lp(s)-kl?>(e) ^ 1 

IK 1 - Cs 1: «)) _1 ||lp(E 1 )^Lp(S 1 ) ^ 1 

Proof. (=>): Suppose that 1 — Ce,™ is invertible on L p (£). Let / G L p (£i) and 
F e L p (£) is its trivial extension to £. Then 3G £ L p (£) such that (1 - Ct,, w )G = 
F. Since w is supported on £1, we can write 

G = F + Cy,, w G = F + Ct, 1iW (G\^ 1 ) 

Letting g — G|si € £ p (£i), we have 

(l-Cz uW )g = f 

so 1 — w is invertible on L p (Y<i). Furthermore, we can easily show |j(l — 
Cs^r 1 !! < ||(1 - Ce,™)" 1 !!, by noticing ||<?||lp( Si ) < ||G|| iP(s) ', \\F\\ LP(S) = 
ll/l|jy>(Ex)) and 

l|G|| LP(s) < ||(1 -Cs,™)- 1 1| -P^i:) 
(<=) Now, suppose that 1 — Cs ltW is invertible on L p (£i). Let F £ L p (£) and 
/ = .F|ei- Let g := (1 — Cy, 1 . w )~ 1 f . If (1 — Cs iW )G — F then similar computations 
as before imply that 

G ( Z \ = U( Z )> ifze£i; 
{F(z) + (Cx uw g)(z), ifze£ 2 . 
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Thus G exists and is unique, furthermore using ||io||oo < 1 we have: 

\\G\\ L , m < + ||F|U, (S2 ) < (l + ||(1 - Cs,^)- 1 !!) ||F|| iP(s) 

as desired. □ 

Remarks: 1. The above proof reveals that if both 1 — C^ llW and 1 — Cs^m are 
invertible and if two functions / S L p (Si) and F E L P (S) agree on Si, then 

((l-C Sl , 10 )- 1 /)(2)= ((I-Ce,™) -1 ^)^) for any zG Si 

(note that the left-hand side is in L p (Si) while the right-hand side is in L P (S), 
however the above equality is about values on Si). Consequently, 

((1 - Ce!,™,)- 1 ^,™/) (z) = ((1 - Ci,™)- 1 ^,™/) (z) for any z e Sj 

Thus, if the weights are in L 2 , the new [i will agree with the old fj, on the old 
contour, provided that the new contour inherits the orientation of the old contour. 
As a corollary, we have 

Corollary 9.2. The trivial extension of w from R to T doesn't affect the a priori 
estimate (which is an estimate onM). 

We'll take our common contour to be T. 

9.2.2. Factorization of jr- Below, we'll determine a natural factorization for jr. 
Let j = (7 — w _ ) _1 (7 + w + ) be the jump matrix associated with w. Inspired by 
the second scheme, our factorization of jr = (I + $_)J(7 + $+)~ 1 will be chosen 
such that 

/ + $+ = {I + w+)- 1 (I + w+) 

/ + $_ = {I -w-)- 1 ^ -vj-). 

This will allow us to show the T analogue of (|8] — ^ [T0|) (using completeness of T) , 
and hence prove the norm equivalence claim. Detailed computation give: 

f(o,o), onr ur 3 ; 

(w^,w+) = l(0,(I + w+)- 1 -I) = (0,-w+), ouriUTa; 

[ (i -(j -to-)- 1 , o) = (-w-, o), onr 4 ur 5 . 

This factorization of jr is natural, it is consistent with our intuition that wjt should 
have a decaying analytic continuation to the + side of T, while should have a 
decaying analytic continuation to the — side of T (see the discussion near the 
beginning of Section [9}. 

Using the standard relation \i = M + (I + w + )~ 1 = M_(7 — w - ) -1 , the above 
choice of factorization will ensure that conjugation from (the trivial extension to 
r) j to jr leaves \i unchanged. 

Corollary 9.3. The contour deformation of the pre-model RHP doesn't affect the 
a priori estimate \4'ty - 
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9.2.3. Contour deformation. Now, we are ready to show that our deformation of the 
pre-model/model RHPs doesn't affect the invertibility of the respective resolvent 
operator: 

Proposition 9.4 (Contour deformation). Let 1 < p < oo. 

(i) For large t, if one of 1 — C Wr and 1 — C w is invertible (on the respective IP 
space) then so is the other. In that case, their inverse norms are "comparable" in 
the sense that the following bounds are equivalent: 

||(1 - C„, r ) _1 ||ip(r)^Lp(r) i$ 1 

||(1 — C w )~~ ||lp(M)->Lp(K) ^ 1 

(ii) The same conclusions are true for the deformation of the model RHP. 

Proof. Note that by Proposition 19.31 it suffices to show Proposition 19.41 for the 
extension of the weights to T. The following argument will be essentially the same 
as in the proof of Lemma 13.31 

(i) We want to show that on L P (T), the invertibility of 1 — Cr, Wr and 1 — Cr.w are 
equivalent and the resolvent norms are comparable. For simplicity, we'll suppress T 
when writing these operators. The desired claim is now an immediate consequence 
of the following analogue of dU — > [10|) 

(l-C*)o(l-C lor ) = l-C™ 

(i-a s ,-i)o(i-c w ) = i-c wr 

(1 - C*) o (1 - C # -i) = (1 - C s -i) o (1 - C») = 1 

where Cj> is the Beals-Coifman operators on T with weights (— $_ , $+), and C$-i is 
the Beals-Coifman operators on T with weights (/—(/ + <&_) _1 , (I + C E>+) _1 — I) = 
(<f>_, — <ji_l_) (thanks to triangularity). The proof of these identities can be done 
similarly, using completeness of T and the following analogue of ([TT]) : For any 
heL p (T), 

(59) (C r h)(z)^(z) = Cr((C+h)^ + -(C^h)^^(z), 

Using (ISTR) and orthogonality of Cp , the rest of the argument is purely algebraic, 
similar to the proof of j5] — > fTU|) . 

The key idea in the proof of ([59| is to exploit the uniform boundedness of $ 
on C \ T (as in the proof of Lemma 13. 3p to make up for the lack of "Hardy space 
continuation" , plus a contour integration argument. 

It suffices to show that for every A € V, 

(C+h)(\)$ + (X) = C+({C+h)* + )(\) 

(Cf fc)(A)*_(A) = -Cr((C£h)$-)(\) 

We'll only show this for the + sign and for A S T2 U T3, the other is similar. 

Let z € {n — a < Arg(z) < n}. For any 7^ positively-oriented triangular contour 
with two edges of length R lying on T2 and T3 (so one vertex is 0), we have 

(CH0W . W - f f (*M1,a 
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for R sufficiently large. Using Carleson's measure theorem and boundedness of it 
is not hard to see that the contribution of the other edge of 7^ (which we'll denote 
by 7b) vanishes as R — > 00 

Consequently 

(C r /»)(z)*(*) = ^ / = C r (ir 2 ur 3 (C r + / l )<i> + )(z) 

2tt« Jr 2 ur 3 A - z V / 

Similar argument for similar contours for other connected components of the posi- 
tive part of C\T) (which will give (on the left-hand side) instead of (Crh)(z)Q(z) 
because z is outside these components), we obtain 

(C r h)(z)$(z)=C r ((C+h)$+)(z) 

so taking nontangential limit to T we get 

(C+h)(\)*+(\) = C+((C+h)<P + )(X) 



for A G T 2 U T 3 , as desired. 

(ii) Similar to (i). □ 

Remarks: The given proof of (TTT1) uses the fact that $± are bounded on R and 
at the same time belong to some Hardy spaces of the upper/lower half planes (so 
that we can use (J7J). This argument can certainly be reapplied to show ([55]) in 
part (i), because in the pre-model case 4> decays uniformly (allowing us to write 
<f> = Cr{&+ — <£>_)). However, we avoided that approach because it is not reusable 
for the deformation of the model weights of part (ii), where strong decay is not 
available for the corresponding $. 

Alternative proof using equivalence of families of RHPs. We largely follows |10j. As 
before, we only prove (i), the argument for (ii) is similar. For any pair of weights 
W~, W + G L°°(T), the invertibility of i - Cw on L P (T) is connected to the unique 
solvability of the following family of RHP: 

m+ = m_(J -W-y 1 ^ + W+) 
m±-f G H P (T±) 

which is parameterized by / G L P (T). Here H p (T±) denotes the complex Hardy 
space Hp of the positive/negative component of C \ T. Essentially, for a fixed /, 
the above RHP can be reduced to the functional equation 

n = f + c w( m, fi g L p (r) 

(m± can be computed from fi by m± = f + C^(fi(W~ + W + )) = ± W ± ).) 
Thus, the unique solvability of the above family is equivalent to the invertibility of 
1 — Cw- 

To show that the invertibility of 1 — C Wr and 1 — C w are equivalent, we'll ex- 
ploit analyticity of both pairs of weights to conjugate back and forth between the 
respective RHP families. Recall the relations 

/ + = {I + w+)- 1 (I + w+) 

/ + = (I — Wp) _1 (J — w~) 
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The extensions of our jump matrices are therefore connected via the conjugation: 

j = (/ + $_)- 1 jr(/ + $+) 

For any / e L P (T), we'll transform the RHP normalized to / in the wr family, 

m + =m_j r , m± e / + H P (T±), 

to a RHP in the w family by defining 

n+ = m+(J + <&+), n_ = m_ (/ + $-) 

(the inverse transformation can be done similarly). Clearly n± satisfies the jump 
relation for w, what's left to show is the existence of some F G L P (T) such that 
n± G F + H p (r±), and F depends only on f and the jump matrices. The last 
requirement is crucial, it basically ensures that F remains a coefficient (and n is 
the variable) of the new RHP; otherwise for any two functions f± £ L P (T) we can 
always write 

(60) f± = (c+(/_) - (/+)) + C±(/ + - /_) 

To find F, observe that we can write m± = f + C^h for h := mr+ — rriv- G 
L p (r) (here the completeness of T is important, it allows for cancelations C£ Cp = 
Cp C+ = 0). Now, 

n+ = f(I + $+) + C+h + (C+/i)$+ 
n_ = /(/ + + Cf ft + (Cp 

Now (JSSJ) implies that (Cp /i)4>± G ff p (r ± ). Applying ([50) to two /i-independent 
functions f$±, we can take 

F = f + C+(/$_) - Cp (/$+) = (1 - 

completing our transformation. 

Now, as discussed above, the conjugation on T from w —> wr doesn't alter /i, 
just like the second perturbation scheme. We therefore obtain 

(1 - C WT )-\f ={l- C w )- l F = (1 - C^-ifl - C«)/ 

Similarly, the reverse direction gives 

(i - c^)- 1 / = (i - c^-^i - c 9 ~i)f 

and again, these identities complete the proof. □ 

II. Reduction on T. In this section, the deformation of the pre-model weights 
will be denoted by wr and the deformation of the model weights will be denoted 
by wr- 

f (o,o), onr ur 3 ; f (0,0), onr ur 3 ; 

= < (0,-wti), onriUr 2 ; (wr,w+) = I (0,-w+), onTiUTs; 

[ (-wm,o), onr 4 ur 5 . [(-wm,o), onr 4 ur 5 . 

For future reference, we'll explicitly compute the potentials ur,«r associated with 
the RHP w r : 

(4) ^0 = "^Off/ r MrW(i2fW + ^W)cfa 
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here Off(A) stands for the off-diagonal part of any 2x2 matrix A. Similarly, 



vr(t) Ur 0J = -^~ os I ^( z )( w r( z ) + w r( z )) dz 



T 

(notice that the T-deformed model weights have strong decay, unlike the original 
model weights on R.) Assuming 

||(1 - C Wr ) _1 || L 2 (r) ^ L 2 (r) < 1, 

we can easily show following analogue of the direct perturbation estimate (|T4]) : 

/7r (w^r + ) - / ^r(u>r+ w r) 
it Jt 

< \\w r - wr||2||Ci„ r J||2 + 1 1 /Ir - Mr||2||w r || 2 + ||w r - w r \\i 
here |j/Ir — /irlb can be similarly controlled by 

||«5r - wv\\z + \\wr ~ w r Uglier — I\\ . for any 2 < q < oo 

1— 2 

Note that the above norms are taken on T. Consequently, for the reduction involving 
leading asymptotics of u(t),v(t) from the deformed pre-model RHP to the deformed 
model RHP, it suffices to show: 

Proposition 9.5 (Weight estimates on L). (i) If p and q have one L 2 derivative 
near then for 1 < p < oo and large t: 

\\wr -Wt\\lp(D ~ t~^ + p )h ^ 

~ — 1 

||«>r|Up(r) + ||^r||if(r) < t p< fc+1 > 
(ii) Ifp and q have two L 2 derivatives near then the first estimate can be improved 
tot- (1+ ^^+\ 

Remarks: Let p — oo, this proposition implies that for large t the resolvent norms 
||(1 — Ctup) -1 !! and ||(1 — C,jj r ) ~ 1 1 1 comparable on any L q . 

Proof, (i) We'll estimate \\wr — wr||i,p(r 1 ) below. Contribution of other Ti can be 
estimated similarly. 

Let 9 be the common phase of the weights on Li (9 = ±8). Then lm(9) > 
on Ti thanks to decaying property of our weights. It suffices to show that on Li: 

(5+ 2 [q] - 5v 2 q(0)y ue < r^TTTe- ct Imfl 

for some absolute constant c > (taking L p (Ti) and using a change of variable, 
we'll get the desired estimate). 

Since \5+ 2 (z) — 8^ 2 {z)\ < \z\? on Li (by our previous study of 5), and 5~ 2 ,5q 2 
are bounded, we have 

(s?[q]-6£q{Q))e ite = 5+ 2 (lq}-q(O)y te + (5+ 2 -5^)q(oy t0 

< (\z\ N +\z\*)e- tlm0 W 

Let u :— \z\, then on Ti we can write thxi9(z) as C(fc)it fe+1 for some C(k) > 0. 
Consequently, we can control the above right-hand side by 

fW) ((u N + u i) e -K(fc)«* +1 ) e -i* 
< t -2wm e -i tIla9 ( z ) t as desired. 



STATIONARY PHASE METHOD FOR OSCILLATORY RHPS 



67 



The estimation of ||n>r||p and ||u>r||p can be done similarly. Using boundedncss 
of So and S on C, similarly it comes down to showing 

l|e- tC|a;|fe+1 ||L,( R )<t-^, C>0, 
which follows by a change of variable. 

(ii) Similar argument, note that we can estimate \5{z) — 8q{z)\ by |z| 1_e on IV □ 

To sum up, regarding resolvent bounds (i.e bounds on (1 — C lu ) _1 ) we completely 
reduce our localized RHP to the model RHP with weights wj^ , while for asymptotics 
of u, v we can reduce it to the T-deformed model RHP w r . 

Reduction of the a priori estimate H43\ ) on T. Notice that our deformation intuitively 
pushes every i£l with |x| > 1 away from the supports of the deformed weights, 
which lie in T \ R. This separation allows us to reapply some ideas from the argu- 
ment in Section [8. 3. II 

The situation is indeed simpler in our case. Observe that the deformed (model/pre- 
model) weights vanishes on R, so for any id s.t. \x\ > 1 we can write 

Hr{x) = I + C r \ R (p r (w^ + Wp))(i) 

Jlr(x) = I + C r \jg,(Jir(w£ + «5p))(i) 

Assuming boundedness of ||(1 — C Wr )~ 1 \\L2(-r)^L 2 {r) (which will be shown in the 
next section), Proposition 19.51 easily implies 

\W - /|U=(d < rwi) 

Now, for ifl such that > 1, we can write Jj-r( x ) ~~ ^t{x) as 
Ct\r{V-tw t ){x) - C r \ R (fj, r w r )(x) 
= Cr\ R ((/ir - /j,t)wt){x) + C r \ R ((^ r - I)(wr - wr))(x) + C r \ R (uT - w r )(x) 

< |~| (llMr - A*r|U 2 (r)||wr||L2(r) + II Mr - I\\L 2 (r)\\w r - w r \\mr) + \\wr - w r \\ L ^( 

_ 3 

< t 2 ( fc +!) , as desired. 



10. Model Riemann-Hilbert problems 

In previous sections, we reduced our oscillatory RHP to the P-deformation of N 
model RHPs, one for each stationary point. Consider one such model RHP. Modulo 
symmetry, we can assume the stationary point is of order k. In this section and 
in the following weights, p and q are constants, they are the values of the original 
p,q at the given stationary point: 




if x 6 D 



(w M ,w+ 



Here is of the form a + bx k+1 with a, b € R and 



if x 6 D_ 



{xeR: G'(x) > 0}, D- = {x e R : 6'{x) < 0}, 
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S(z) := exp(iw + /3(z)) is defined on C\l with boundary values 6±(x) on R. Here, 

!0, if is an exterior point of £>_ ; 

iei/ In [eA] , if is an endpoint of D_ ; 

— 7ri/sgn[Im(A)], if is an interior point of D-. 

{0, if k is even; 

sgn (6), if k is odd. 

v = -— \n(l+pq) 

and w is a real number, which we can assumed by incorporating it into p, q. The 
deformed weights were naturally chosen to be 



(w r (z), Wjt(z)) = < 



(0,0), zer ur 3 ; 

(0, analytic cont. of — lujflg,), z G T\] 
(0, analytic cont. of — u>^lg_), z G 



(analytic cont. of u^1r_, 0), z G 
(analytic cont. of «;^-1r + , 0), z £ I5. 

These can be made explicit by considering the parity of k and the sign of b. 

Our task is to show the boundedness of (1 — C^)^ 1 in L 2 (T), the a priori 
estimate for /i r := I + (1 — C wr )~ 1 C WT I on a set P C K bounded away from Xj, 
and verify the desired asymptotics for those potentials recovered from Mr- Recall 
that the boundedness of (1 — C^) -1 on L 2 (T) is equivalent to the boundedness of 
(1 — C WM )~ 1 on L 2 (R) as proved in the last section (Proposition I9.4[) . 

Recall that 1 + pq > 0. We'll divide our analysis into the following cases: 

1. Degenerate (or linear) case, when pq = 0. 

2. Defocusing case, when — 1 < pq < (which corresponds to a defocusing 
system [33]: q = —p). 

3. Focusing case, when pq > 0. 

10.1. Degenerate and defocusing cases. It is probably not too hard to check 
the boundedness of (1 — C^,,)" 1 in the degenerate case pq — 0, however we'll 
prove the unique solvability of our model RHP in these cases using the following 
proposition: 

Proposition 10.1 (Resolvent bound). If — 1 < pq < 1 then 

|(1 - C WM ) 1 ||l 2 (R)^L 2 (R) ^p : q 1 

Proof. We call a pair (p, q) "good" if the above conclusion is true. Then the propo- 
sition is a consequence of the following two lemmas: 

Lemma 10.1. If (p,q) is good then so is (cp,c~ 1 q) for any c > 0. 

Lemma 10.2. If\p\, \q\ < 1 then (p,q) is good. 

First, it is not hard to sec that (w^, w^) can be written as: 



(61) 
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where A(t, x) = i(tQ(x) — 2ev\n\x\j. Note that A is purely imaginary on the 
real line, and depends on the product pq rather than individual p, q, so it remains 
invariant under the change (p, q) — > {cp, c~ 1 q). 

The first lemma now follows from three simple observations: 

(i) Vc > 0, the equality / = (1 — C WM )~ 1 F is equivalent to 

C CT3/ V C "- /2 =(l-^3/W-/ 2 )"(^ /2 ^ /2 ), -3 -(J 

(ii) The pair of weights c' 73 ^ 2 wmc^' 73 ^ 2 can be obtained from the pair of weight 
u>m if we replace (p, q) by (cp,c~ 1 q). 

(hi) For any g € L p , the LP norms of g and c" 3 ^ 2 gc~ a3 l 2 are comparable up to 
a harmless constant depending on c. 

Below we'll show the second lemma. Observe that the model jump matrix can 
be rewritten as 

'l+pq pe~ lte " 
qe it& x 

Using analyticity and uniform boundedness of S^ 1 , we can use a similar argument 
as before (Proposition 19. 4|) to deform the RHP associated with Jm to the RHP 
associated with the following jump matrix: 



Jm 



J 



1 + pq pe lte 
qe ite x 



More precisely, we can show that ||(1 — C WM ) 1 ||l 2 (r)-s-l 2 (R) is comparable to ||(1 — 
CV m ) _1 ||l 2 (K)-s-l 2 (R)j where the L°° pair of weights Wm is defined by 

I + W+ = (I + w+)5-°° , I - W+ =(I- w M )SZ as 

We can also think of W M as being obtained from a factorization of J, and recall that 
the invertibility and the inverse norm of (1 — Cw M ) doesn't depend on the particular 
factorization (as long as every factors and their inverses remain bounded). So it 
suffices to show the desired bound (on the resolvent norm) for one factorization 
only, which we'll take to be 

i.e. W- = (1 "^'1 , W* - I 



v J ' " ^e lte y 

Now assume that \p\, \q\ < 1. For any L 2 matrix-valued function / = (/y)f we 
can explicitly compute Cwf as 

(C-(f 12 qe ite ) C+(f llP e- lte ) 
Cwf - {C-(h 2 qe^) C + (hipe- ltB ) 

Consequently, using the fact that \\C± ||l 2 (r)->.l 2 (r) = 1, we have 

||CV/|| 2 <max(b|,M)||/|| 2 

so 1 — Ctjv is invertible on L 2 (R), with 

||(1 - Ci4/) _1 ||l 2 (R)^L 2 (K1 < 



1 - max(|p|, \q\) 

as desired. We note that if \pq\ are sufficiently small then L p boundedness is also 
true. □ 
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Remarks: When fc = 1 it is possible to show the above resolvent bounds on L p , 
2 < p < oo using a suitable hypergeometric function that solves a timeless model 
RHP equivalent to the above RHP. This is the approach taken in [TU] and [21] ■ 

10.1.1. Computation of model potentials. Assuming unique solvability of the model 
cases, we'll show the following result: 

Proposition 10.2. The model potentials recovered from M-p are of the form 

( o pur 



u r (t) 
v r (t) 



qVt~ w exp(2A) 

1 . ievlnt 



"*+t exp(-2A) 




with A := — ita 
2 



1 



for absolute constants U, V depending on pq, b, and k such that U = — V . 

Proof. To compute the potentials recovered from Mr, we'll make the following 
change of variable: 

i 

M r ^ N r (Q = cxp (Aer 3 )M r (z) exp ( - Acr 3 ) 

here A := ^ita+ * is a i-dependent constant. We'll abuse notation and use the 
same symbol T for the respective contour of £. Also, D± still denote the respective 
sets for £. 

Then Nr solves the L 2 -normalized RHP on V whose weights are analytic con- 
tinuation to T of the pair (6^ f (£), i>Jf(£)) defined as follows: 



(62) 



where = ib^ k+1 — 2ieu\w\^\. The analytic continuation rule is exactly as 

before, and the deformed jump matrix for £ is 




(63) (&f (£),&+(£)) = I 



(0, 0), 

(0, analytic cont. of — 6| / 1r + ), 
(0, analytic cont. of — 6^1r_), 



eer ur 3 ; 
(er i; 

?er 4 . 



(analytic cont. of 
(analytic cont. of 

The analytic continuation of A(£) from (the respective side of) K to T can be 
computed explicitly, in particular, we'll have A(£) — —A(£) for any £ G V. 

The model RHP associated with N-p is still normalized in L 2 (T) sense with 
L? n L°° weights, furthermore it is timeless i.e. the weights are independent of 
t. This RHP is clearly uniquely solvable, thanks to the unique solvability of Mr. 
Indeed, the O(l) bound on ||(1 — C6 r ) _1 |j can be seen as a particular case of the 
same bound on ||||(1 — C«; M ) || when t = 1. Consequently, we can recover the 
following "potentials" from Np by sending £ — > oo nontangentially in C \ T: 







NT 



Unt 




= lim#V r (0 = --^-Off 
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Note that UnFiVnt are constants depending only on p,q,k,b (it is not hard to 
see that these constants are independent of the contour T, as long as the deformed 
angle is small enough - anyway our chosen angle depends on k). 
Consequently 

( m Ur n ] ) = lim ^rW 

\vr{t) / 2^00 

= t^fc+r exp ( - Aer 3 ) lim £iV r (£) exp (Aa 3 ) 

£—►00 

= ^ exp ( - A. 3 ) C/ - r )exp(A. 3 ) 

/ U Nr t^T^r cxp(-2A) 

\v Nr t~T^ exp(2A) 

To finish up, we'll show that there exist constants U, V such that 

U N r = pU, V N r = qV, 
and U, V depend only on pq, k, b and are anti-complex conjugates: 

U = -V 

The existence of U, V is indeed natural from and the strict triangular structures 
of the weights (where p is always part of the top right and q is always part of the 
bottom left entry). 

To see that U,V depend on pq (instead of individuals p and g), observe that 
Vc G C \ {0}, the transformation (p, q) M> (c 2 p,c~ 2 q) corresponds to the following 
changes: 

b r i->- c a3 brc~ a3 , ^ Nr H> d 73 ^ NT c~' 73 

which imply that the diagonal entries of /jjvr are unaffected, so it leaves U, V 
invariant. 

Below, we'll show that U, V are anti-complex conjugates. By 

V = ^-J^ T e A ^g 2 {z)dz 

where g%, g 2 take values in {—1, 0, 1} and they depend on D-, D + : 

• If Z?_ = M_ and D + = W.+ then g\ = lr 2 + lr 5 , and g 2 = lri + lr 4 - 

• If D- = and D + = M_ U R + then g 2 = l Tl + lr 5 and g 1 = l r , 2 + lr 4 
The remaining cases are similar. By considering all cases it is not hard to see that 

= 52 (0- Thus, it suffices to show that 

For simplicity, we'll suppress the subscript NT in /Lijvr in our computation. Since 
br = on ro ur 3 , the operator Cb r is technically operating on functions supported 
on T° :=r\M. 

For convenience of notation (which will be explained shortly), we'll orient T° by 
orienting \Ti\, [T^ outwards and | | inwards. Then, 

c br f = c br j v/GL p (r) 
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here C& is the Beals-Coifman operators on T° with the following pair of weights: 



(&Fo(0,#o(0) = 



{(analytic cont. of 6^1r + , 0), £ € Ti; 

(analytic cont. of 6^1e_, 0), £ G T 2 ; 

(analytic cont. of — 6^1r_, 0), £ G T 4 . 

(analytic cont. of - b M l R+ , 0), f £ T 5 . 



As can be seen, 6p~ = 0, so it is more notationally convenient to work with Cb v0 ■ 
This is the main advantage of &p° over &r- 

Let g : r° — >• {±1} defined by <? = lr^u^ — lr 4 ur 5 . Unraveling the equality 
/i = J + Cb r ^ = J + Cb r0 /x we get 

fi 11 = l + qT 2 p 12 , n 12 =pT l ^ 11 

,,22 _ -I , 21 ,,21 _ _7i ,,22 

where Ti , T 2 denote the following operators 

Tif = C+ (fe- A 9l g) , T 2 f = C+ (fe A g 2 g) 

From here it is not hard to see that xt 11 = 1 + pqT 2 Tip n and xt 22 = 1 + pqTiT 2 ^i 22 . 
For any / G L P (T) (1 < p < oo) we denote by / the function f(z) = f(z). We'll 
show that 

(65) (T^\f)(z) = (T 2 Tj)(z) 

Since pg G R, this equality and uniqueness of jj, automatically implies /x 11 = /i 22 , 
as desired. To see (I55|) . notice that our previous discussion gives A = —A, while 
clearly g = —5. Consequently, together with the previous observation g\ — g 2 we 
have 

T\f(z) = -T 2 f(z), and Trf[z) = -Tj{z) 
for any z G T, which easily imply (|65[) . □ 

10.1.2. ^4 priori estimate for pr- We'll show the following estimate: 
Proposition 10.3. For any i£l such that \x\ > 1, me /iawe 

A*r(x,t)=/ + * U0 M 

Proof. We'll make the same change of variable as in the proof of Proposition 110.21 

Z \-> £ = t Z 

M r !-)• N r (0 = exp (cr 3 A) Af r (z) exp ( - a 3 A) 

Under this change of variable the weights will become &p which are defined by ([53"|) . 
As before, we use the same symbol T for the respective contour of £ (which is really 
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a dilation of T z ). Since bf = on R = T U T 3 , V£o G M \ {0} we have 

Mwr(Co) = / + — ? — t ^ 

^772 Jr\R 4 - 4o 

7 | 1 /r\R^r(Q6r(Qdg | 1 f ^(OMO^ 
27ri -Co Jr\R £o(£-fo) 

/ o l 



Vjvr 
. Co 



+ 0(7^) 



o ; v iCoi 2 - 

notice that, thanks to boundedness of (1 — Cf,,,)" 1 and strong decay of 6r, 
||Miv-r(0 & r(0^lli|(r) ^ I^Jvr - 7||L2(r)||£&r|U 2 (r) + ||£&r|Ui(r) < 1 
and the trivial inequality 



Finally, recover x from £o we obtain the desired estimate (notice that A is purely 
imaginary so multiplication by e ±2A won't destroy our error bound). □ 



10.2. Focusing cases. As mentioned above, the quadratic (focusing) case when 
k = 1 can be studied using parabolic cylinder functions, see for instance |31j . In 
this section, we'll discuss the cubic (focusing) case, i.e. when k — 2 and pq > 0. 

The main idea is to look at our deformed model RHP as an inverse monodromy 
problem (in the sense of Birkhoff [3 ) with one irregular singularity at oo of Poincare 
rank fe+1 [22]). Essentially, the solvability of this inverse problem can be understood 
by studying the location of the poles of suitable deformation equations, which exist 
if k>2 [13]. 

In the cubic case (i.e. k = 2), it is classical that these deformation equations form 
a system of two Painleve II equations, and under the imaginary assumption pq > 
they can be shown to be free of real poles. This approach and the corresponding 
result go back to Its-Novokshenov [21] (see also [4] for a more direct perspective) . 

To see how to formulate our deformed model RHP as an inverse monodromy 
problem, we first repeat the same variable change z i— > £ to obtain a timeless 
normalized RHP (N r ,jr) with jump matrix j' r (f) = (7 - 6f (f)) -1 ^ + &r(0)) 
where br are defined by (|63p . In other words, 



analytic cont. of I — 
jr = \ analytic cont. of 7 — b~t[ 1r_ , 
analytic cont. of 7 — 6^1k_, 
analytic cont. of 7 — 6^1r + , 



?er„ur 3 ; 

£Gr 2 ; 
£Gr 4 ; 
£Gr 5 . 



with bu defined by |62|) . For simplicity of notation, let c = |. 
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If c > 0, using (|62|) the jump matrix jr can be written explicitly as 



Jr 



I, 




—pe 2lc f 1 
1 



eer ur 3 ; 
eeriur 2 ; 

£er 4 ur 5 . 



while if c < the picture if essentially reflected across R and can be handled 
similarly. So without loss of generality we'll assume that c > 0. We can think of 
our model RHP as a search for a 2 x 2 matrix valued function 

-.= N r (0e~ icSh+1(T3 

analytic in the six sectors formed by T with asymptotics e~ lc ^ + <73 , so that we can 
travel from one sector to the other using connection matrix e %c ^ <T3 jr (£) e_ " 3 ■ 
If we fix the travel direction to be counter clockwise (instead of from the + side 
to the — of r as we usually do), then we have six connection matrices defined on 
To, • ■ • , T 5 : 

1 0W 1 0\ (I -p\ (I ^ 
y q l) \-q 1) ' \Q 1 j ' 1,0 1, 

(observe that Ti and are oriented inwards while I^ and T4 are oriented out- 
wards). This suggests an inverse monodromy problem in the sense of Birkhoff [3J, 
although some adaptations will be required. Let fifc denote the sector: 



I, 



(k 



- 2)77 fc77 

- < argz < — } 



S21 — 



I = 1,2,3 



for each 1 < k < 7. Consider the following set of Stokes matrices 

'1 s 2l \ f 1 X 

,0 1 J> S *-i- \ S2l _ x 1, 

satisfying the cyclic relation S1S2 ■ ■ ■ Se — I (in general the right-hand side is 
e -2ni~fa 3 £ or a m0 nodromy exponent 7, but in our case 7 will be 0). 

Theorem 10.3 (Its-Novokshenov) . For all but a discrete subset of x G R, there 
exists uniquely a family of (2 x 2 matrix valued) functions ^k(;x) holomorphic in 
C, such that has the following asymptotics inside flk 

mi{x) m2(x) 



* fc (£,aO = (I- 



)e 



-i(c£ +x£)<t 3 



t e 

and are related via ^l>k+i(£,) — ^k(£,)Sk in &k H Qk+i for any 1 < k < 7. 
Furthermore, if u(x) :=2m\ 2 (x) and v(x) :=2m 21 (x) thenu,v satisfy 

= 2u 2 v+^u 



2v 2 u- 



3c 



and the above discrete subset o/R is exactly the set of real poles ofu(x). 

Remarks: 1. The above asymptotics should be understood as follows: Given any 
strict subsector Q,' k C fife and any sufficiently small neighborhood in C of x, we 
have 

7711(2;) 77l„(2;). 



* fc (£,x)e j(c « (I + 



-) = o(t 
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More over, this asymptotics is differentiable, i.e. similar bounds are true for the 
derivatives of ^k- 

2. The above system of two Painleve equations is called the deformation equation 
for the respective isomonodromy problem. 

Sketch of proof [JJ [50] . For convenience, we sketch the proof of the above theorem, 
more details can be found in the references. First, using a classical theorem of 
Sibuya |29) . we can always find a weak solution to the above inverse problem, i.e. 
a family $^ satisfying the required properties in a neighborhood of infinity. Fix k. 
Using Birkhoff-Grothendieck's theorem, we can factorize 

**(0 = r it - 1 (0r fcO1, *fc(0 

where is a nonnegative integer depending on x, and Tk,^k & re holomorphic 
invertible functions on C \ {0} and C respectively. By a suitable left multiplication 
and an appropriate normalization, we can assume that ^k has an asymptotics of 
the form 

(66) x) = (I+ + + . . . ) e -Kd s +^-in k taflo* 

? ? 

inside fi^. Here, as a result of the left multiplication, the 12-entry of rrij will be 
zero for every 1 < j < Ink- If 2rik > fco + 1 = 3 it is possible to show that 
the formal monodromy exponent nk associated with the above ^k is exactly by 
using a recursive formula of Jimbo, Miwa, and Ueno [22j . giving a contradiction. 
Consequently, = or 1. 

On the other hand, using a generalized version of the Birkhoff-Grothendieck the- 
orem (initially proved by Malgrange [25] , elementary proof given later by Bolibruch 
[4]) we can show that rik(x) are the same for all but a discrete set of x € R. This 
means for these x's the monodromy of the associated ODE are the same, hence 
using a compatibility condition we can obtain the following deformation system 

U xx = 2u 2 w + |f u 
v xx = 2v 2 u + |f v 

with u(x) :— 2m\ 2 (x) and v(x) := 2m 21 (x). Using a suitable Schlesinger transfor- 
mation and the Painleve property of the above system, we can show that the set 
of x where nk(x) = 1 are exactly the poles of the above system. Now, for those x 
with nk{x) = 0, it is not hard to see that ^fc are exactly what we need. □ 

To apply this theorem to our case, we first notice that by a symmetry observation 
as in the proof of Proposition [TUT] we can assume that p = q. Let ^k's satisfy the 
conclusion of the above theorem for the following values of Si , . . . , S$ : 

C U U M U 

Recall the second Pauli matrix 02 = f ? . Using p — q, we observe that 

o~2SkO-2 — Sj} k V 1 < k < 6 

Thus, if we define * fe (£) by $ fe (0 := o- 2 *7-fc(0c r 2 VI < k < 6 and $7 = *i, then 
^k also satisfies the conclusion of the theorem. Consequently, = ^k, hence 
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m\{x) = <72mi(a;)cr2) which implies 

u(x) — —v(x) for any 
Thus, u solves the following imaginary Painleve II equation: 

Ax 

u xx = -2\u\ 2 u + —u 
3c 

By looking at the Laurent series of u, it is clear from the above equation that u 
doesn't have any pole on K. Consequently, the above inverse monodromy problem 
is solvable for any real x, in particular for x = (which is what we are interested 
in). Since S% = S$ = I, the asymptotics (|66l) is true on C+ for 'J 2 = ^3, and on 
C_ for * 5 = * 6 . Now, define iV r on C \ T by 

*i(0, for £ e r 50 u r oi; 

*a(0, foreer 12 ; 

* 4 (0. for e £ r 23 u r 34 ; 

It is clear that Np satisfies the required jump relation on T with jump matrix jp, 
furthermore using strong decay of N-p and contour integration, we can easily show 

(67) N r (0 = I + (C r g)(0, for £ G C \ T 

here g(x) := Nr+(x) — Ny-(x). Clearly g G L P (T) for any 1 < p < 00. Let unf be 
defined on T by 

Mat = N r+ (I + b+) _1 = JV r - (I - b-) _1 
From (|57)l . it is not hard to see that //jvr satisfies: 

Furthermore, /Ltjvr S I + L P (F) for any 1 < p < 00. This is because fiNr is 
asymptotically / + 0( j|j ) and at the same time continuous on T. We can recover 
A*rOM) from hnt(0 by 

Hr(x,t) = (J-Nritw x) 

Notice that /ir is invertible as a matrix (indeed its determinant is 1) and fir, Mr 1 £ 
L°°(r). The invertibility (and boundedness of the inverse) of 1 — C Wr on L^(T) now 
follow from the T analogue of Proposition [3TT1 which can be summarized briefly as: 

Corollary 10.4. For any 1 < p < 00, 1 — C Wv is invertible on LP(T) with bounded 
norm: 

II 1 - C™ r IUi(r)_>L£(r) < p 1 

In particular, when t = 1 we have ||1 — C& r ||,LP(r)->L p (r) 1- The rest of the 
computations (i.e. computing asymptotics of u(t), v(t) and proving the a priori 
estimates for fi) can be done exactly as in Section 110. 1.11 and Section 110.1.21 
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